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We review the field of cavity optomechanics, which explores the interaction between electromag- 
netic radiation and nano- or micromechanical motion. This review covers the basics of optical 
cavities and mechanical resonators, their mutual optomechanical interaction mediated by the ra- 
diation pressure force, the large variety of experimental systems which exhibit this interaction, 
optical measurements of mechanical motion, dynamical backaction amplification and cooling, 
nonlinear dynamics, multimode optomechanics, and proposals for future cavity quantum optome- 
chanics experiments. In addition, we describe the perspectives for fundamental quantum physics 
and for possible applications of optomechanical devices. 



Contents 



I. Introduction [2] 

II. Optical cavities and mechanical resonators 0] 

A. Optical resonators [4] 

1. Basic properties [4] 

2. Input-output formalism for an optical cavity [5] 

B. Mechanical resonators [H 

1. Mechanical normal modes [71 

2. Mechanical dissipation [8l 

3. Susceptibility, noise spectra and fluctuation 
dissipation theorem [9] 

III. Principles of optomechanical coupling 1101 

A. The radiation pressure force and optomechanical 
coupling IIQI 

B. Hamiltonian formulation 1111 

C. Optomechanical equations of motion 1131 

IV. Experimental realizations and optomechanical 
parameters ^] 

A. Optomechanical parameters 1141 

B. Suspended mirrors 1161 

C. Optical microresonators 1181 

D. Waveguides and photonic crystal cavities 1181 

E. Suspended and levitated nano-objects 1191 

F. Microwave resonators 1191 

G. Ultracold atoms |20] 

V. Basic consequences of the optomechanical 

interaction |20] 

A. Static phenomena: Optical potential and bistability |20] 

B. Dynamical backaction |^ 

1. Optical spring effect |22] 

2. Optomechanical damping rate 1231 



*Electronic address: markus. aspelmeyer@univie.ac. at Electronic 

address : opt omechanicsr mp@gmail . com 

^Electronic address: tobias.kippenberg@epfl.ch 

■1- Electronic address: Florian.Marquardt@physik.uni-erlangen.de 



VI. Quantum optical measurements of mechanical 



motion 1251 

A. Parametric displacement sensing and the standard 
quantum limit (SQL) 1251 

1. Introduction and qualitative discussion 1251 

2. The standard quantum limit 1261 

B. Optical QND measurements 1281 

1. Single quadrature measurements [29] 

2. Mechanical Fock state detection [29] 

3. Optical feedback cooling (cold damping) [30] 

VII. Optomechanical cooling [31] 

A. Quantum theory of radiation pressure cooling [32] 

B. Strong coupling regime 1361 

1. Optomechanical normal- mode splitting 1361 

2. Optomechanically induced transparency 1391 

3. Optomechanically induced amplification and 
electromagnetically induced absorption 1401 

VIII. Classical Nonlinear Dynamics 1411 

A. Parametric instability and attract or diagram 1411 

B. Chaotic dynamics 03] 

IX. Multimode optomechanics 1431 

X. Quantum Optomechanics 1451 

A. Light-assisted coherent manipulation of mechanics 1461 

B. Mechanics-assisted readout and manipulation of light l46l 

C. Optomechanical entanglement 27] 

D. Quantum hybrid systems 09] 

E. Quantum protocols 1501 

F. Nonlinear quantum optomechanics 1521 

XI. Future Perspectives 1551 

A. Foundational aspects [55] 

B. Applications [56] 

XII. Acknowledgements [53 

XIII. Appendix: Experimental Challenges 1571 

A. Influence of classical excess laser phase noise on laser 

cooling [57] 



B. Influence of cavity frequency noise on laser cooling [58] 



2 



C. Influence of classical excess laser noise on sideband 

thermometry 1581 



References 



I. INTRODUCTION 



EH 



Light carries momentum which gives rise to radiation 
pressure forces. These forces were already postulated in 
the 17*^ century by Kepler, who noted that the dust tails 
of c omets point aw ay from the sun during a comet tran- 
sit (Kepler 1619). The first unambiguous experimental 



demonstrations of the radiation pressure force predicted 
by Maxwell were performed using a light mill configura- 
tion (iLebedewl 119011 iNichols and Hulll 119011). A careful 



analysis of these experiments was required to distinguish 
the phenomenon from thermal effects that had domi- 
nated earlier observations. As early as 1909, Einstein 
derived the statistics of the radiation pressure force fiuc- 



tuations acting on a moveable mirror (Einstein 1909), 



including the frictional effects of the radiation force, and 
this analysis allowed him to reveal the dual wave-particle 
nature of blackbody radiation. In pioneering experi- 
ments, both the linear and angular momentum transfer of 
photons to atoms and macros copic objects were demon - 
strated by Frisch ( |Frisch[ 11933] and by Beth ( |Beth[|1936D , 
respectively. 

In the 1970s Arthur Ashkin demonstrated that focused 
lasers beams can be used to trap and control dielectric 



particles, which also included feedback cooling (Ashkin 



1978 2006). The non-conservative nature of the radia- 



tion pressure force and the resulting possibility to use it 
for cooling atomic motion was first pointed out by Han- 
sch and Schawlow and by Dehmelt and Winelan d (Han- 
sch and Schawlow 1975 Wineland and Dehmelt [ |1975| T 
Laser cooling was subsequently realized experimentally 
in the 1980s and has beco me since then an extraordi- 
narily important technique (Stenholm 1986). It has, for 



example, allowed cooling of ions to their motional ground 
state and it is the underlying resource for ultracold atom 
experiments. Many appl ications have be en enabled by 
laser cooling ( Metcalf and van der Straten |1999j , includ- 
ing optical atomic clocks, precision measurements of the 
gravitational field, and systematic studies of qu antum 
many-body physics in trapped clouds of atoms (Bloch 
and Zwergerl|2008[). 



The role of radiation pressure and its ability to pro- 
vide cooling for larger objects was already investigated 
earlier by Braginsky in the context of interferometers. 
Braginsky considered the dynamical influence of radia- 
tion pressure on a harmonically suspended end-mirror 
of a cavity. His analysis revealed that the retarded 
nature of the force, due to the finite cavity lifetime, 
provides either damping or anti-damping of mechani- 
cal motion, two effects that he was able to demon- 
strate in pioneering experiments using a microwave cav- 
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Figure 1 Schematic of a generic optomechanical system, both 
in the optical domain (top), with a laser-driven optical cav- 
ity and a vibrating end mirror, as well as in the microwave 
domain (bottom), with a vibrating capacitor. Here we have 
depicted a microwave drive entering along a transmission line 
that is inductively coupled to the LC circuit representing the 
microwave resonator. 



observed in microwave-coup l ed fe^ -scale mechanical res- 
onators (Cuthbertson et al. 1996). Independently, sim- 



ilar physics was explored theoretically for solid-state vi- 



brations (Dykman 



1978 ). In the optical domain, the first 



cavity optomechanical experiment (|Dorsel et al 1983) 



demonstrated bistability of the radiation pressure force 
acting on a macroscopic end-mirror. 

Braginsky also addressed the fundamental conse- 
quences of the quantum fluctuations of radiation pres- 
sure and demonstrated that they impose a limit on how 
accurately the pos ition of a free test mass (e.g. a mirror) 



can be measured (Braginsky and Khalili 1995[ Bragin- 



[sky and Manukin 1977). A detailed analysis by Caves 
clarified the role of this ponderomotive quantum noise in 
interferometers (Caves 1980). These works established 



the standard quantum limit for continuous position de- 
tection, which is essential for gravitational wave detectors 
such as LIGO or VIRGO. 

During the 1990s, several aspects of quantum cavity 
optomechanical systems started to be exp lored theoret 
ically. These include squeezing of light (Fabre et al. 



1994 Mancini and Tombesi 1994) and quantum non- 



it}^ Braginsky and ManukinJ |1967| [Braginsky et al. et al. 1994 Pinard et al. 1995), which exploit the effec 



demo l ition (QND) detect ioiTof t he light intensity ([Jacobs 



1970). In later experiments, these phenomena were also tive Kerr nonlinearity generated by the optomechanical 
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interaction. It was also pointed out that for extremely 
strong optomechanical coupling the resulting quantum 
nonlinearities could give rise to nonclassical and entan- 
gled states of the light field an d the mechanics ("Bose| 
let aLjl997 Mancini et al. 1997). Furthermore, feedback 



([Thompson et al] '2008) and nanorods (iFavero et 



cooling by radiation pressure was suggested (Mancini 



et al. 1998). Around the same time, in a parallel de- 
velopment, cavity-assisted laser-cooling was proposed as 
a method to cool the motion of atoms and molecules that 



et all 



2009[ ) insid e Fabry-Perot resonators, whispering gallery 

2QQ9t 

and microspheres ( Ma et aL\ 2QQ7| Park and Wang , 20091 
Tomes and Carmonj " 2009 ) , photonic crystals ( Eichenfield] 
2009a|bD, and eva nescently coupled nanobeams 



fAnetsberger e t al] |200 9y In addition, cavity optome- 
chanics has been demonstrated for the mechanical ex- 



1998 Vuletic and Chu 2000) 



lack closed internal transitions ( Hechenblaikner et al Murch et al 2008j). Optomechanical interactions are 
^ also present in optical waveguides - as first studied and 



citations of cold at om clouds ( jBrennecke et al 



2008 



On the experimental side, optical feedback cooling 
based on the radiation pressure force was first demon- 



strated in (Cohadon et al 1999) for the vibrational 



modes of a macroscopic end-mirror. This approach was 



later taken to much lower temperatures (Kleckner and 



Bouwmeester 2006 Poggio et al 2007). At the same 



time, there was a trend to miniaturize the mechanical 
element. For example, the thermal motion of a mm- 
scale mirror was monitored in a cryogenic optical cavity 
( [Titt onen et al 1999). Producing high-quality optical 
Fabry-Perot cavities below that scale, however, turned 
out to be very challenging. In spite of this, it was still 
possible to observe optomechanical effects of retarded ra- 
diation forces in microscale setups where the forces were 
of photothermal origin, effectively replacing the cavity 
lifetime with a thermal time constant. Examples in- 



clude demo nstration of the opt ical spring effect ( Vogel 
et al. 2003), feedback damping (Mertz et a/. ,1993*), self- 



induced oscillations (Hohberger and Karrai^ 2004^ 



Zala- 



lutdinov et al. 2001), and cavity cooling due to the dy- 
namical backaction of retarded photothermal light forces 
(Hohberger Metzger and Karrai 2004). 



Yet, for future applications in quantum coherent op- 
tomechanics it is highly desirable to be able to exploit the 
non-dissipative radiation pressure force. Both the advent 
of optical microcavities and of advanced nanofabrication 
techniques eventually allowed to enter this regime. In 
2005 it was discovered that optical microtoroid resonators 
with their high optical finesse at the same time contain 
mechanical modes and thus are able to display optome- 
chanical effects, in particular radiation-pressure induced 
self-oscillations (Carmon et al.\ |2005[ Kippenberg et al.\ 



2005| [Rokhsari et al. 2005)(i.e. the effect Braginsky 
termed "parametric instability"^). In 2006 three differ- 
ent teams demonstrated radiation-pressure cavity cool- 



ing, for suspended micromirrors (Arcizet et al. 2006a 



Gigan et al. 2006 ) and for microtoroids ( Schliesser et al. 



2006p . Since then, cavity optomechanics has advancec 



rapidly and optomechanical coupling has been reported 
in numerous novel systems. These include membranes 



^ Braginsky called the process of dynamical backaction ampli- 
fication and the concomitant self-induced coherent oscillations 
"parametric oscillatory instability", as this effect is undesirable 
in gravitational wave interferometers which were the basis of his 
analysis. 



observed in the context of squeezing, where the confined 
mechanical modes of fibers lead to Guided Acoustic Wave 
scattering (jShelby et al. 1985). Nowadays there are a 
number of systems where such optomechanical interac- 
tions are explored in the absence of a cavity, such as 
waveguides in photonic circuits or photonic cr ystal fi- 
bres, see e.g. (Kan g et a/.| |2009[ |Li et a/.||2008D . These 
setups lie somewhat outside the scope of the concepts 
presented in this review, but we emphasize that they are 
very promising for applications due to their large band- 
width. 

Optomechanical coupling has also been realized using 
microfabricated superconducting resonators, by embed- 
ding a nanomechanical beam inside a superconducting 



transmission line microwave cavity (Regal et al. 2008) 



or by incorporating a fiexible aluminum membrane into a 



lumped element superconducting resonator ( Teufel et al. 



2011bj). In these systems the mechanical motion capaci- 
tively couples to the microwave cavity. This approach 
ties cavity optomechanics to an independent develop- 
ment that has been garnering momentum since the late 
1990s, which is concerned with measuring and controlling 
the motion of nano- and micromechanical oscillators us- 
ing electrical and other non-optical coupling techniques. 
Examples include coupling of mechanical oscillators to 



single electron transistor s ( Cleland et al. 2002 LaHaye 
et a/.[ 2Q04a[ Naik et al. 2006 D or a quantum point con - 



tact ([Cleland et a/. [ |2002| [Flowers- Jacobs et aT\\2mi ). 



Besides a wealth of possible applications for such de- 



vices in sensitive detection (Cleland and Roukes 1998 



LaHaye et al. 2004b Rugar et al. 2004), these meth- 



ods provide the possiblility of realizing mechanical quan- 



Schwab and Roukes 2005) by direct interaction with 


two- level quantum systems ( Arcizet et al. 


2011 Cleland 


and Gellerl 2004| Kolkowitz et al.\ [2012| 


jaHaye et al.\ 


2009||O'Connell et a/.l|2010||Rugar et al. 


20041 Wilson-j 


Rae et al.\ 2004[). For recent comprehensive general re- 



views of nanomechanical systems (in particul ar electro- 



mechanical devices), we refer the reader to ([Blencowe 



2005 



[Greenberg eTalj [2QT2| [Foot and van der Zant 



2012). 

It should be noted that in atomic systems quantum 
coherent control of mechanical motion is state of the art 
since early pioneering experiments with trapped ions 



for reviews see ([Blatt and Wineland[ [2008[ Jost et al. 
2009 Leibfried et al.\ 2003). In fact, quantum informa- 
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tion processing in these systems relies on using the quan- 
tum states of moti on to mediate interaction s between dis- 
tant atomic spins (Cirac and Zoller 19951). In contrast, 



the fabricated nano- and micromechanical structures that 
form the subject of this review will extend this level of 
control to a different realm, of objects with large masses 
and of devices with a great flexibility in design and the 
possibility to integrate them in on-chip architectures. 

There are several different motivations that drive the 
rapidly growing interest into cavity optomechanics. On 
the one side, there is the highly sensitive optical detec- 
tion of small forces, displacements, masses, and accelera- 
tions. On the other hand, cavity quantum optomechan- 
ics promises to manipulate and detect mechanical motion 
in the quantum regime using light, creating nonclassical 
states of light and mechanical motion. These tools will 
form the basis for applications in quantum information 
processing, where optomechanical devices could serve as 
coherent light-matter interfaces, for example to intercon- 
vert information stored in solid-state qubits into flying 
photonic qubits. Another example is the ability to build 
hybrid quantum devices that combine otherwise incom- 
patible degrees of freedoms of different physical systems. 
At the same time, it offers a route towards fundamental 
tests of quantum mechanics in an hitherto unaccessible 
parameter regime of size and mass. 



A number of reviews (Aspelmeyer et al.\ 2010| 


2012 


Cole and Aspelmeyer |2012t |Favero and 


Karrai| 


2009 


Genes et a/. , 2009a; Kippenberg and Vahala, 2007, 20081 


Marquardt and Girvin^ ,2009^ Meystre^ .2012^ ^Schliesser] 


and Kippenberg |2010) and brief commentary 


)apers 


(ICleland; "2009"; "Cole and Aspelmeyer 2011 


Karrai 


2006 



published during the past few years, and the topic 
has also been discussed a s part of a larger reviews on 
nanomechanical systems (Greenberg et al. 2012| |Poot| 
[and van der Zant[ [2012). Here we aim for a compre- 
hensive treatment that incorporates the most recent ad- 
vances and points the way towards future challenges. 

This review is organized follows: We first discuss op- 
tical cavities, mechanical resonators, the basic optome- 
chanical interaction between them and the large range of 
experimental setups and parameters that are now avail- 
able. We then go on to derive the basic consequences 
of the interaction (such as optomechanical damping and 
the optical spring effect), describe various measurement 
schemes, and present the quantum theory of optomechan- 
ical cooling. After studying nonlinear effects in the clas- 
sical regime, we address multimode setups and the wide 
field of proposed applications in the quantum domain, 
before concluding with an outlook. 



II. OPTICAL CAVITIES AND MECHANICAL 
RESONATORS 

In this section we recall the basic aspects of optical cav- 
ities and of mechanical resonators, as needed to describe 



cavity optomechanical systems. Much more about these 
topics can b e found in standard textb ooks on quantum 
optics, e.g. (Walls and Milburn 1994), and on nanome- 



(Walls and Milburn 


1994) 


bems (Cleland 2003 


). 



A. Optical resonators 

Optical resonators can be realized experimentally in 
a multitude of forms of which several types will be dis- 
cussed later in the review. Here we give a unifying de- 
scription of the optical properties and provide the math- 
ematical description of a cavity that is pumped with a 
single monochromatic laser source. 



1. Basic properties 

We first consider the classical response of a simple 
Fabry-Perot resonator, which will allow to introduce the 
relevant parameters to characterize an optical cavity. A 
Fabry-Perot resonator or etalon consisting of two highly 
reflective mirrors, separated by a distance L, contains a 
series of resonances which are given by the angular fre- 
quency ujm ~ m-TTj^. Here m is the integer mode number. 
The separation of two longitudinal resonances is denoted 
as the free spectral range (FSR) of the cavity: 



FSR 



C 



(1) 



In the following, we will almost always focus on a single 
optical mode, whose frequency we will denote oJcav 

Both the finite mirror transparencies and the internal 
absorption or scattering out of the cavity lead to a finite 
photon (intensity) cavity decay rate^ hz. 

A further useful quantity is the optical finesse, J^, 
which gives the average number of round-trips before a 
photon leaves the cavity: 



J-. 



(2) 



The optical finesse is a useful parameter as it gives the 
enhancement of the circulating power over the power that 
is coupled into the resonator. Alternatively, we can in- 
troduce the quality factor of the optical resonator. 



Q. 



opt 



(3) 



where r = is the photon lifetime. Note that the qual- 
ity factor is also used to characterize the damping rate of 
mechanical resonators (see below). Generally speaking, 
the cavity decay rate tz can have two contributions, one 



^ In this review we shall use k for the photon (energy) decay rate, 
such that the amplitude decay rate is given by hi/2. In some 
papers the latter is denoted as n. 
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Symbol 


Meaning 




Mechanical frequency 


Vm 


Mechanical damping rate 


Qm 


Mechanical quality factor, Qm = ^m/^m 


A 


Laser detuning from the cavity resonance, A = cjl — ojcav 


K 


Overall cavity intensity decay rate, from input coupling and intrinsic losses, k = + kq. 


90 


Optomechanical single-photon coupling strength, in Hint = —hgoa^a(b-\- ¥) 


9 


Light-enhanced optomechanical coupling for the linearized regime, g = go a/ ricav 


G 


Optical frequency shift per displacement, go = Gxzpf 


XZPF 


Mechanical zero-point fluctuation amplitude, xzpf = y^h/2mef£^in 


a 


Photon annihilation operator, with a the circulating photon number 


b 


Phonon annihilation operator, with ¥b the phonon number 


n 


Average number of phonons stored in the mechanical resonator, n = (b^b^ 


nth 


Average phonon number in thermal equilibrium, nth = (^g^^m/fcs^ _ 


ncav 


Photon number circulating inside the cavity, ncav = ^a^a^ 




Optical susceptibility of the cavity, Xaa{(^) = — ^(cj + A)]~^ 




Mechanical susceptibility, Xxx{(^) = niefi{^m — cj^) — imeffFmCJ 




Quantum noise spectrum, Sxx{(jo) = f dte^^^ (x(t)x(O)) (Sec. ILB.3) 




(Symmetrized) mechanical zero-point fluctuations, Sxx^ {(^) = h\lmxxx((^)\ 




^xt^{^) ^ ^xx^i^)'- standard quantum limit result for added noise in displacement measurement (Sec. VLA) 


Topt 


Optomechanical damping rate (Sec. V.B.2): max. Aficavgo/i^ ^ ^ 




Optical spring (mechanical frequency shift. Sec. V.B.I): 2ncavgo^ / ^ for n <^ Qm, A 


n-niin 


Minimum reachable phonon number in laser-cooling, nmin = {K,/4Qrn)'^ for k, <^ 



Table I List of most important symbols and some formulas used in this review 



from losses that are associated with the (useful) input 
(and output) coupling and a second contribution from 
the internal losses. It is useful to differentiate these two 
contributions. For the case of a high-Q cavity, the total 
cavity loss rate can be written as the sum of the individ- 
ual contributions: 



Here, K:ex refers to the loss rate associated with the 
input coupling, and tvo refers to the remaining loss rate. 
For example, in the case of a waveguide coupled to a mi- 
crotoroidal or microsphere resonator, f<iex is the loss rate 
associated with the waveguide-resonator interface and 
describes the light absorption inside the resonator. For 
the case of a Fabry-Perot cavity, f<iex is the loss rate at 
the input cavity mirror and a^o summarizes the loss rate 
inside the cavity, including transmission losses at the sec- 
ond cavity mirror as well as all scattering and absorption 
losses behind the first mirror. Note that by splitting the 
total decay rate into these two contributions, we are as- 
suming that the photons going into the hzo decay channel 
will not be recorded. More generally, one could distin- 
guish between more decay channels (e.g. input mirror, 
output mirror, absorption). 



2. Input-output formalism for an optical cavity 



A quantum mechanical description of a cavity that is 
coupled to the outside electromagnetic environment can 
be given either via master equations (if only the inter- 
nal dynamics is of interest) or via a framework known 
as input-output theory, if one also wants to access the 
light field being emitted by (or reflected from) the cav- 
ity. Input-output theory allows us to directly model the 
quantum fluctuations injected from any coupling port 
(such as the input mirror) into the cavity. In addition, 
it takes into account any coherent laser drive that may 
be present. For more deta ils beyond the brief discussion 
provided below, see e.g. (Clerk et al 2010a Gardiner 



and Zoller 



2004|. 



Input-output theory is formulated on the level of 
Heisenberg equations of motion, describing the time- 
evolution of the field amplitude a inside the cavity. One 
finds that the amplitude a experiences decay at a rate 
K:/2. At the same time, its fluctuations are constantly 
replenished via the quantum noise entering through the 
various ports of the cavity. In the present case, we dis- 
tinguish between the channels associated with the input 
coupling (decay rate A^ex) and the other loss processes 
(overall decay rate /^o, including loss through the second 
mirror). The equation of motion reads: 



6 



(a) 



Optical / Microwave Waveguide 

— — r 



Uni directional 
Optical mode 



(b) 




n 



(0 
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Bi-directional 
Optical nnode 
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Figure 2 Comparison of uni- and bi-directionally coupled cav- 
ities and the notion of reflected and transmitted field am- 
plitude, (a) denotes a waveguide- coupled unidirectional res- 
onator (e.g. WGM cavity coupled to waveguide), (b) denotes 
a single-sided cavity coupled in reflection (or a double-sided 
cavity where the transmission signal is disregarded), (c) de- 
notes a waveguide coupled bi-directional cavity, which can 
decay in both forward and backward propagating waveguide 
modes (i.e. tapered fiber coupled photonic crystal mode, or 
a waveguide-coupled quarter- wave stripline resonators), (d) 
denotes the coupling to a double-sided Fabry Perot resonator, 
in which both transmitted and reflected field are measured. 



(4) 



In the classical case, a would be replaced by a prop- 
erly normalized complex amplitude of the electric field of 
the cavity mode under consideration. Indeed, the clas- 
sical version of this equation (and the following ones) 
can be obtained by simply taking the average, such that 
a ^ {a). We have chosen a frame rotating with the laser 
frequency cjl, i-e. a°^^^ = e 
duced the laser detuning A = 



I, here havc iutro- 
ujl — u jca,v with respect to 



lar equation can also be written down for the mechanical 
oscillator in order to describe its dissipation and the as- 
sociated noise force, comprising quantum and thermal 
contributions (see Sec. |IILC ). 

The input field ain{t) should be thought of as a stochas- 
tic quantum field. In the simplest case, it represents the 
fluctuating vacuum electric field coupling to the cavity 
at time t, plus a coherent laser drive. However, the same 
formalism can also be used to describe squeezed states 
and other more complex field states. The field is normal- 
ized in such a way that 



p = hujc 



is the input power launched into the cavity, i. e. y^jn^: 

is the rate of photons arriving at the cavity. The same 
kind of description holds for the "unwanted" channel as- 
sociated with /in- 

According to the input-output theory of open quantum 
systems, the field that is reflected from the Fabry Perot 
resonator (or coupled back into the coupling waveguide) 
is given by: 



^out — 



(5) 



the cavity mode (see also Sec. III.B). Note that a simi- 



Note that this input-output relation describes correctly 
the field reflected from the input mirror of a Fabry- 
Perot resonator. The above equation describes also the 
transmitted pump field of an evanescently coupled uni- 
directional waveguide resonator system, such as a whis- 
pering gallery m ode resonator coupled to a waveguide 
(Cai et a/.' "2000). In this case the above expression would 
yield the transmitted pump field. 

We still have to consider the case of a two-sided cavity, 
e.g. a two-sided Fabry Perot cavity. Other examples in 
this review include a waveguide coupled to superconduct- 
ing stripline cavities or fiber-taper coupled photonic crys- 
tal defect cavities. In these cases there are both trans- 
mitted and reflected fields. In all of these cases there are 
two options for the description. If the field transmitted 
through the second mirror is not of interest to the anal- 
ysis, one may lump the effects of that mirror into the 
decay rate a^o, which now represents both internal losses 
and output coupling through the second mirror. If, how- 
ever, the field is important, it should be represented by 

/ (2) (2) 

an additional term of the type \/K,ex d\J in Eq. (4). Then 
an equation analogous to Eq. ^ will hold for the output 
field a^^^^ at that second mirror. 

In the following, we will not be concerned with noise 
properties, but focus instead on classical average quan- 
tities (for a single-sided cavity), taking the average of 
Eqs. and ([5|. 

We can solve the equation Q first for the steady-state 
amplitude in the presence of a monochromatic laser drive 
whose amplitude is given by (ain). Noting that (/in) = 0, 
we obtain: 
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(a) = 



iA 



(6) 



The expression linking the input field to the intracavity 
field will be referred to as the optical susceptibility, 



Xaa{oj) 



1 



-i{uj + A) + 1^/2 

Thus, the steady-state cavity population ncav = (a^a), 
i.e the average number of photons circulating inside the 
cavity, is given by: 

-ca. = |(a)r = ^^^^£^ (7) 

were P is the input power launched into the cavity. The 
reflection or transmission amplitude (for the case of a 
Fabry-Perot cavity or a waveguide-coupled resonator, re- 
spectively) can be calculated by inserting Eq. (|6| into 
Eq. (|5|. Using the symbol TZ for the reflection amplitude 
in the sense of figure [2] case (b), we obtain: 



(^out) _ (^0 - /^ex)/2 - iA 
{din) {i^o + /^ex)/2 - iA 



(8) 



The square \7l\ of this amplitude gives the probability of 
reflection from the cavity (for Fabry-Perot) or transmis- 
sion in the case of a uni-directional waveguide resonator 
system. From this expression, several regimes can be dif- 
ferentiated. If the external coupling /^ex dominates the 
cavity losses {hCq^ ^ ^ hcq), the cavity is called "over- 
coupled". In that case \7lf ~ 1 and the pump photons 
emerge from the cavity without having been absorbed or 
lost via the second mirror (a property that is important 
as discussed below in the context of quantum limited de- 
tection). The case where a^o = i^ex refers to the situation 
of "critical coupling". In this case, 1Z{A = 0) = on res- 
onance. This implies the input power is either fully dissi- 
pated within the resonator or fully transmitted through 
the second mirror (in the case of a Fabry-Perot cavity 
with hZQ denoting the decay through the second mirror). 
The situation tZex is referred to as "undercoupling" 

and is associated with cavity losses dominated by intrin- 
sic losses. For many experiments this coupling condition 
is not advantageous, as it leads to an effective loss of 
information. 

The physical meaning of reflection (or transmission) 
depends sensitively on the experimental realization under 
consideration. One can distinguish four scenarios, which 
are outlined in the figure [2] 

B. Mechanical resonators 

1. Mechanical normal modes 

The vibrational modes of any object can be calculated 
by solving the equations of the linear theory of elasticity 



under the appropriate boun dary conditions that are de- 
termined by the geometry^ (Cleland 2003). This eigen- 
value problem yields a set of normal modes and cor- 
responding eigenfrequencies r^(n)- The mechanical dis- 
placement patterns associated with mechanical motion 
are given by the strain field Un{r)^ where n designates 
the normal mode. 

For the purposes of this review, we will mostly focus 
on a single normal mode of vibration of frequency Q.^ 
(where 'm' stands for 'mechanical'), assuming that the 
mode spectrum is sufficiently sparse such that there is no 
spectral overlap with other mechanical modes. The loss 
of mechanical energy is described by the (energy) damp- 
ing rate F^, which is related to the mechanical quality 
factor^ by Qm = ^m/^m- If one is interested in the equa- 
tion of motion for the global amplitude x{t) of the mo- 
tion, one can utilize a suitably normalized (see below) 
dimensionless mode function u{f^t)^ such that the dis- 
placement field would be u{f^t) = x{t) • u{r). Then the 
temporal evolution of x{t) can be described by the canon- 
ical simple equation of motion of a harmonic oscillator of 
effective mass meffi 



^eff- 



dx(t) 

- rrief^Tm—^ + mef^^l,x{t) 



Fe,t{t) (9) 



Here Fext{t) denotes the sum of all forces that are acting 
on the mechanical oscillator. In the absence of any exter- 
nal forces, it is given by the thermal Langevin force (see 
Sec. II. B. 3). In the above equation the (energy) damping 
rate F^ has been assumed to be frequency independ ent. 
Deviations of this model are treated for example in ( |Saul-| 



son 



1990). 



sary at this point (Cleland 2003 



A brief remark about the effective mass meff is neces - 

1999]). 



Pinard et al. 



The normalization that has been chosen for the mode 
function u{f) affects the normalization of x{t). However, 
it will always be true that the potential energy is given by 
meff^m (^^(^)) /2. This value can then be compared to 
the expression for the potential energy that arises from 
a calculation according to the theory of elasticity. De- 
manding them to be equal yields the correct value for 
the effective mass meff (which therefore is seen to de- 
pend on the normalization that was chosen for the mode 
function). Of course, for the simple case of a center-of- 
mass oscillation of a solid object, a natural definition of 
x{t) is the center of mass displacement in which case the 
effective mass will be on the order of the total mass of the 
object. A treatment of effective mass in optomechanical 



experiments is found in (Pinard et al. 1999[) 



Eq. (|9| can be solved easily, which is best done in fre- 
quency space. We introduce the Fourier transform via 



^ A powerful simulation approach in this context are finite element 

(FEM) simulations. 
^ In the context of mechanical dissipation often the loss tangent 

S^is quoted, its relation to the quality factor being Qm = j^- 
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x{oj) = f_^dte"^'x{t). Then Sx{co) = XxA<^)Fe.t{uj) 
defines the susceptibihty Xxx, connecting the external 
force to the response 5x of the coordinate: 



(10) 



The low frequency response is given by Xxx{^) = 
{meff^'^)~^ = 1/k where k is the spring constant^. 

The quantum mechanical treatment of the mechanical 
harmonic oscillator leads to the Hamiltonian 

Here the phonon creation (b^) and annihilation (b) oper- 
ators have been introduced, with 



X = xzpf(^ + b^), P = -imeffl^m^ZPF(^ - b^), 



where 



XZPF 



h 



the mechanical vacuum state. The position and momen- 
tum satisfy the commutator relation [xzpf^Pzpf] = 
The quantity b^b is the phonon number operator, whose 
average is denoted by n = (bW). In the following, we will 
typically not display explicitly the contribution ^hQm of 
the zero-point energy to the energy of the oscillator. 

We briefly discuss the effect of dissipation. If the me- 
chanical oscillator is coupled to a high temperature bath, 
the average phonon number will evolve according to the 
expression: 



dt 



(n) = -Tm{{n) - nth) 



For an oscillator which is initially in the ground state, 
{n){t = 0) = this implies a simple time dependence 
of the occupation according to {n){t) = n^^(l — e"^/^"^), 
where nth is the average phonon number of the environ- 
ment. Consequently, the rate at which the mechanical 
oscillator heats out of the ground state is given by: 



d 
dt 



{n{t = 0)) = hth ■ r. 



th 



^ To describe the response of a high Q oscillator near resonance 
uj r^mone can approximate Xxx by a Lorentzian, i.e. using 
Q'^ -ixP- = {Vim -u;)(r2m +w) ~ 2(^7rn - ^S)VLm yields Xxx(c^) = 



The latter is often referred to as the thermal decoher- 
ence rate, and given by the inverse time it takes for one 
quantum to enter from the environment. In the above 
expression the high temperature limit has been taken, 
i.e. hth ^ ksTiyath/^^m' This expression shows that to 
attain low decoherence a high mechanical Q factor and a 
low temperature bath are important. The change of pop- 
ulation of a certain Fock state can be described within 
the framework of the Master equation approach. This ap- 
proach allows to calculate the decoherence rate of other 
quantum states s uch as a Foc k state \ n). T he latter is 
given by (see e.g. (Gardiner and Zoller 2QQ4[ )): 



(n + l)nth^m + n{nth + l)^m 

revealing that higher Fock states exhibit a progressively 
higher rate of decoherence. 



2. Mechanical dissipation 

The loss of mechanical excitations, i.e. phonons, is 
quantified by the energy dissipation rate Tm = ^m/Qm- 
The origins of mechanical dissipation have been inten- 
sively studied over the last decades and comprehensive 



is the zero-point fluctuation amplitude of the mechan- 
ical oscillator, i.e. the spread of the coordinate in the 

ground-state: (0|x^|0) = x|pp, and where |0) denotes include: 



reviews are found for example in (Cleland 2QQ3[ Ekinci 



and Roukes 2005). The most relevant loss mechanisms 



• viscous damping, which is caused by interactions 
with the surrounding gas atoms or by compression 



of thin fluidic laye rs (Karabacak et all 2007 Ver- 
bridge et a/T| |2008[ [V^nola et a/.| |2"QQ6]); 



clamping losses, which are due to the radiation of 
elastic waves into the substrate through the sup- 



ports of the oscillator (Anetsberger et al 


2008 


IBindel and Govindjee||2005t Clark et al 


2005 


ICole 


et al\ |2011| Cross and Lifshitzl |2001| 


^ichenfield 


et a/.||2009bt|Jockel et aL\ 20111 lJudge et a/.l|2007 



Judge[ [2004, |Wangeram2000, .Wilson-Rae^ ,2008^ ; 



• fundamental anharmonic effects such as thermoe- 
lastic damping (TED) and phonon-phonon interac- 
tions (Duwel et oD, 2006 ; 'Klselev and lafrate 2008 



Lifshitz and Roukes, ,2000| |Zener[ 1 1938^ 



• materials-induced losses, which are caused by the 
relaxation of intrinsic or extrinsic defect states 



in the bulk or surface of the resonator (|Mohanty 


et aL\ 20021 |Southworth et aL\ |2009t 


Unterreith- 


meier et al] 


2010) IVenkatesan et al. 


|2010| |Ya- 


sumura et al. 


2000). Such losses have been suc- 



involving two level defect states, which are coupled 
to the strain via the deformation potential ([Ander- 



son et al.\ |1972t |Hunklinger et al.\ |1973t |Phillips| 


|1987t|Remus et al.\ 


20091 Seoanez et a/.l|2008t|Tiel-| 


biirger et al. 


1992 


). In the context of nano- and 
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Figure 3 Brownian motion (thermal fluctuations) of a 
nanomechanical resonator in the time-domain (schematic), 
with amplitude and phase fluctuating on a time scale set by 
the damping time F" ^ . 



micromechanical oscillators the two level fluctua- 



tor damping has been revisited ( Remus et al. 2009 



Seoanez et al. 2008). 



The various dissipation processes contribute indepen- 
dently to the overall mechanical losses and hence add 
up incoherently. The resulting mechanical quality fac- 
tor Q total is given by ^ = XI where i labels the 
different loss mechanisms. 

Another helpful quantity is the so-called "Q • /" prod- 
uct, which plays an important role in the phase noise 
performance of oscillators. In the context of optome- 
chanics, it quantifies the decoupling of the mechani- 
cal resonator from a thermal environment. Specifically, 
' fm ^ { ) J denotes the number of coher- 



ent oscillations in the presence of thermal decoherence, 
and evidently scales with Q ■ f. 



3. Susceptibility, noise spectra and fluctuation dissipation 
theorem 

If one measures the motion of a single harmonic oscil- 
lator in thermal equilibrium, one will observe a trajec- 
tory x{t) oscillating at the eigenfrequency ftrn- However, 
due to the influence of both mechanical damping and 
the fluctuating thermal Langevin force, these oscillations 
will have a randomly time- varying amplitude and phase. 
Both amplitude and phase change on the time scale given 
by the damping time F~^. Such real-time measurements 
have been p erformed in optomechanical systems (JHadjaij 
et fl/.[|1999[ ) (see Fig.[4|. 



In experiments, the mechanical motion is often not an- 
alyzed in real-time but instead as a noise spectrum in 
frequency space. This allows to easily separate the con- 
tributions from different normal modes. We briefly re- 
capitulate the relevant concepts. Given one particular 
realization of the trajectory x{t) obtained during a mea- 
surement time r, we define the gated Fourier transform 
over a finite time interval r: 



x{uj) 



(11) 



Averaging over independent experimental runs, we obtain 
the spectral density (|x(cc;)p). In the limit r ^ oo, the 




Figure 4 Brownian motion of a mechanical resonator ob- 
tained in an optomechanical setup. The two "quadratures" 
x{t) = Xi{t)cos{Qmt) + X2{t)sin{Qmt), are displayed in a 
frame rotating at the angular mechanical resonance frequency 
Qrn, such that the unperturbed undamped motion would cor- 
respond to a stationary single point sitting at the origin of 
the phase space. The fluctuations are a consequence of the 
thermal Brownian mo tion (f igure courtesy of A. Heidmann, 
see also (Hadjar et al. 1999)). 



Wiener-Khinchin theorem connects this to the Fourier 
transform Sxx{^) of the autocorrelation function, also 
called the noise power spectral density: 



Here we have defined: 



+ CO 



{x{t)x{0)) e'""* dt . 



(12) 



(13) 



The only assumption which has been made is that x{t) is 
a stationary random process. From Eqs. ( 12|13 ), we im- 
mediately obtain the important result that the area un- 
der the experimentally measured mechanical noise spec- 
trum yields the variance of the mechanical displacement. 



/ + 00 
-co 



(14) 
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Figure 5 Noise spectrum of a damped harmonic oscillator in 
thermal equilibrium (symmetric in cj i-^ — cj in the classical 
limit ksT > hw). 



Furthermore, in thermal equilibrium, the fluctuation- 
dissipation theorem (FDT) relates the noise to the dissi- 
pative part of the linear response. 



(15) 



where Xxx{^) denotes the mechanical susceptibility intro- 
duced above and we have treated the high-temperature 
(classical) case. For weak damping (F^ <C ^m)^ this 
gives rise to Lorentzian peaks of width F^ in the noise 
spectrum, located at cj = i^^rn (see Fig. [5|. Integra- 
tion of Sxx{^) according to Eq. (14) yields the variance. 



which for weak damping is set by the equipartition the- 
orem: (x^) = ksT /mQf^Vt^. 

In the quantum regime, the natural generalization of 
Eq. (13) contains the product of Heisenberg time-evolved 



operators, (x(t)x(O)), which do not commute. As a conse- 
quence, the spectrum Sxx{^) is asymmetric in frequency. 
The quantum FDT 



Sxx{^) 



2h 



1 



-huj/kBT 



lmxxx{^) 



(16) 



implies that Sxx{^) = for cj < at T = 0. Our discus- 
sion of dynamical backaction cooling will mention that 
this means the T = bath is not able to supply en- 
ergy, as there are no thermal excitations. In this review 
we will also consider the symmetrized noise spectrum, 
Sxx{^) = {Sxxi'^) + Sxx{—^)} /2. Fo r more on noise 
spectra, we refer to ( [Clerk etal. 2010a). 



III. PRINCIPLES OF OPTOMECHANICAL COUPLING 

A. The radiation pressure force and optomechanical 
coupling 

In our discussion the fundamental mechanism that cou- 
ples the properties of the cavity radiation field to the me- 



chanical motion is the momentum transfer of photons, i.e. 
radiation pressure. The simplest form of radiation pres- 
sure coupling is the momentum transfer due to reflection 
that occurs in a Fabry Perot cavity. A single photon 
transfers the momentum |Ap| = 2h/\ (A: photon wave- 
length). As a consequence the radiation pressure force is 
given by 



(F) = 2hk-^ '- = h-{d)a) 

Tn L 



Here Tc = 2L/c denotes the cavity round trip time. 
Therefore, hj^ describes the radiation pressure force 
caused by one intracavity photon. The parameter G = 
uj/L which appears in this expression describes also the 
change of cavity resonance frequency with position, i.e. 
the frequency pull parameter. In the next section, which 
introduces a Hamiltonian description of the interaction 
between a movable mirror and optical cavity, this relation 
will be derived in its full generality. 

More generally, the optomechanical coupling can arise 
for example by direct momentum transfer via reflec- 
tion (Fabry-Perot type cavities with a moveable end- 
mirror, microtoroids) , by coupling via a dispersive shift 
of the cavity frequency (membrane in the middle, levi- 
tated nano-objects trapped inside the cavity) or by opti- 
cal near-field effects (e.g. nano-objects in the evanescent 
field of a resonator or a waveguide just above a substrate). 

Various radiation pressure forces have been investi- 
gated in the pioneering work of Ashkin, who first demon- 
strated th at small dielectric particles can be trapped in 
laser light rt Ashkin , 2006). The relevant forces are gener- 



ally referred to as gradient (or dipole) forces, as the force 
arises from the gradient of the laser field intensity. The 
particle is attracted to the center of the Gaussian trap- 
ping laser beam. If E{r) denotes the laser electric field 
distribution, the time-averaged dielectric energy of the 
particle in the field is given hy U = —^p-E = —^a\E{f)\'^ 
(with a the polarizability) which correspondingly yields 
a force F = —VU. In addition to the gradient force, scat- 
tering forces occur for a traveling wave. These forces scale 



with 



i.e. the wavenumber of the electromagnetic ra- 
in contrast to the gradient forces. In addition 



diation 

there is also a contribution from the strain-optical ef- 
fect, i.e. the strain-dependent polarizability. The strain- 
optical coupling is the dominant c oupling mechanism 
in guided acoustic wave scattering ( Locke et al\ 1998 



Shelby et al |1985 ). Independent of the physical inter- 
pretation of the force, however, the optomechanical in- 
teraction in an optomechanical system can always be de- 
rived by considering the cavity resonance frequency shift 
as a function of displacement (i.e. the "dispersive" shift). 
This will be the basis for our Hamiltonian description 
adopted in the next section. 

It is important to note at this point that a significant 
difference between the trapping of particles in free space 
and micromechanical systems is the fact that the lat- 
ter are also subject to radiation forces based on thermal 
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effects. Absorption of light can heat a structure and de- 
form it, which corresponds to the action of a force (e.g. in 
an asymmetric, bimorph structure, including materials of 
different thermal expansion). These photothermal forces 
can in many ways lead to effects similar to retarded radi- 
ation pressure forces, with the thermal relaxation time of 
the structure replacing the cavity photon lifetime. How- 
ever, since such forces are based on absorption of light, 
they cannot form the basis for future fully coherent quan- 
tum optomechanical setups, since at least the coherence 
of the light field is thereby irretrievably lost. 



B. Hamiltonian formulation 

The starting point of all our subsequent discussions 
will be the Hamiltonian describing the coupled system 
of a radiation mode interacting with a vibrational mode 
(Fig.[T]). For brevity we will refer to the radiation field as 
"optical", even though the important case of microwave 
setups is included here as well. 

We will focus here on the simplest possible model sys- 
tem in cavity optomechanics, which has been used to 
successfully describe most of the experiments to date. In 
this model, we restrict our attention to one of the many 
optical modes, i.e. the one closest to resonance with the 
driving laser. Moreover, we also describe only one of the 
many mechanical normal modes. This is mostly arbi- 
trary, as the displacement frequency spectrum will show 
peaks at any of the mechanical resonances. Still, as long 
as the dynamics is linear with independently evolving 
normal modes, the model will provide a valid approxi- 
mation. In some cases, like sideband-resolved cooling, 
it may be possible to experimentally select a particu- 
lar mechanical mode by adjusting the laser detuning, 
whereas in other cases, like nonlinear dynamics, an ex- 
tended description involving several mechanical modes 
may become crucial. 

The uncoupled optical (c^cav) and mechanical {ftrn) 
mode are represented by two harmonic oscillators, which 
is typically an excellent approximation at the displace- 
ments generated in the experiments: 



Ha hur 



(17) 



In the case of a cavity with a movable end mirror the 
coupling of optical and mechanical mode is parametric, 
i.e. the cavity resonance frequency is modulated by the 
mechanical amplitude^: 



+ xdujcci^/dx + . . . 



Note that such a setup is also considered for discussions of the 
dynamical Casimir effect, where cavity photons are created by 
the mechanical modulation of the boundaries. In the optome- 
chanical scenarios considered here, however, the mechanical fre- 
quencies are too small for this effect to play a role. 



For most experimental realizations discussed in this re- 
view, it suffices to keep the linear term, where we de- 
fine the optical frequency shift per displacement as G = 
—dcuceiv/dx (but see Sec. VI.B.2| for another example). A 
more detailed derivation of the optomechanical Hamilto- 
nian can be found in an early paper (Law 1995). 



We mention in passing that other coupling mechanisms 
have been discussed. For example, the transparency of 
a moving Bragg mirror , and hence z^, can depend on its 
velocity ( [Karrai et al.\ 2QQ8| ). More generally, the dis- 
placement may couple to the cavity decay rate, yielding 

= i<i{x). This case (sometimes termed "dissipative cou- 
pling' ' ), whic h is of practical relevance in some setups (Li 
et a/.[|2QQ9c), can lead to novel physical effects, e.g. in 
cooling plste et a/. [[20091 ). 

For a simple cavity of length L, we have G = uJca,v/L. 
The sign reflects the fact that we take x > to indicate an 
increase in cavity length, leading to a decrease in cjcav(^) 
if G > 0. In general, expanding to leading order in the 
displacement, we have: 



^cav(^)^^^ ^ ^((^cav — Gx)o) CL . 



(18) 



Here x = xzpf(^ + ^^j^ as defined before. Thus, the 
interaction part of the Hamiltonian can be written 



where 



go GxzpF 



(19) 



(20) 



is the vacuum optomechanical coupling strength, ex- 
pressed as a frequency. It quantifies the interaction be- 
tween a single phonon and a single photon. We stress 
that, generally speaking, is more fundamental than G, 
since G is affected by the definition of the displacement 
that is to some extent arbitrary for more complic ated me - 
chanical normal modes (see the discussion in Sec. II.B.l). 
Therefore, in the following we will almost always refer 
to go- Further below, we will also mention ^, which is 
an often-used measure for the effective optomechanical 
coupling in the linearized regime. It will be enhanced 
compared to g^ by the amplitude of the photon field. 
The Hamiltonian reveals that the interaction of a mov- 
able mirror with the radiation field is fundamentally a 
nonlinear process, involving three operators (three wave 
mixing) . 

The radiation pressure force is simply the derivative of 
^int with respect to displacement: 



dx 



hGo) d = h- 



XzPF 



(21) 



The full Hamiltonian H will also include terms that 
describe dissipation (photon decay and mechanical fric- 
tion), fluctuations (influx of thermal phonons), and driv- 
ing by an external laser. These effects are formulated 
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most efficiently using the equations of motion and the 
input-output formaUsm (see Sec. 11. A. 2 , and also the next 
section). Here, we just remark that it is convenient to 
change the description of the optical mode by switching 
to a frame rotating at the laser frequency ul. Applying 
the unitary transformation U = ex.p{iujLa^ at) makes the 
driving terms time-independent^, and generates a new 
Hamiltonian H = UHoidU^ + ihdU /dt of the form 

H = -HAd^a + nnrnb^b - hgod^d{b -^P)^ ... , (22) 
where 



mechanical oscillator driven optical cavity 



rm(^th + 1) 




Figure 6 Schematic: Optomechanical (linearized) interaction 
between a driven optical mode and a mechanical resonator. 



(23) 



is the laser detuning introduced already in Sec. |II.A.2[ 
and where we have omitted (...) driving, decay, and fluc- 
tuation terms, which will be discussed below. Eq. (22) is 



the frequently used starting point in cavity optomechan- 
ics. 

We now introduce the so-called "linearized" approxi- 
mate description of cavity optomechanics. To this end, 
we split the cavity field into an average coherent ampli- 
tude (a) = a and a fiuctuating term: 



d = a^Sd (24) 
Then, the interaction part of the Hamiltonian 



-hgo{a + Sd^ia + Sd){b + b^) (25) 



may be expanded in powers of a. The first term, 
—hgo \df (b -\- b^), indicates the presence of an average 
radiation pressure iovce F = HG \a\. It may be omitted 
after implementing an appropriate shift of the displace- 
ment's origin by Sx = F /mQf^Q'^. The second term, of 
order |q;|^, is the one we keep: 



hgQ{a*5d + a5d^){b + b^) 



(26) 



The third term, —hgoSd^Sd, is omitted as being smaller 
by a factor \a\. Without loss of generality, we will now 
cav real-valued. Thus, the Hamiltonian 



assume a = ^/n, 
in the rotating frame reads 



where the quadratic interaction part 



(27) 



^iJf ^ = -%ov^(Mt + 5a){b + P) (28) 



is referred to as "linearized", since the resulting coupled 
equations of motion will be linear in this approximation. 
Note that the remaining terms in Eq. (27) no longer con- 
tain the laser driving, as that has already been taken care 
of by the shift implemented in Eq. (24). In the literature 
up to now, the combination 



(29) 



is often referred to as "the optomechanical coupling 
strength". Obviously, it depends on the laser intensity 
and is thus less fundamental than the single-photon cou- 
pling ^0 (obtained for ncav = !)• In the linearized regime 
described here, the optomechani cal system can be vi ewed 
in analogy to a linear amplifier ( [Hotter et al. 2012) that 
receives optical and mechanical input fields. 

The linearized description can be good even if the av- 
erage photon number circulating inside the cavity is not 
large. This is because the mechanical system may not be 
able to resolve the individual photons if the decay rate 
n is sufficiently large. The detailed conditions for the 
linearized approximation to be valid may depend on the 
questions that are asked. We will return to this ques- 
tion in the section on nonlinear quantum optomechanics 
(Sec. 



X.F) 



We briefiy note that ^ > is one neccessary condi- 
tion for the so-called "strong coupling" regime of cavity 
optomechanics, where the mechanical oscillator and the 
driven optical mode hybridize (Sec. VH.B). A much more 
challenging condition is to have g^ > i.e. the single- 
photon optomechanical coupling rate exceeding the cav- 
ity decay rate. In the latter regime, nonlinear quantum 
effects will become observable (see Sec. X.F). 



Depending on the detuning, three different regimes can 
be distinguished with respect to the interaction (28), 



especially in the sideband-resolved regime <C ^m-, 
which we assume in the remainder of this section). For 
A ^ — we have two harmonic oscillators of (nearly) 
equal frequency that can interchange quanta: the me- 
chanical oscillator and the driven cavity mode. Within 
the rotating-wave approximation (RWA) we thus can 
write the interaction as 



hg{6d^b^ 6dP) . 



(30) 
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This is the case relevant for cooling (transferring all ther- 
mal phonons into the cold photon mode; Sec. VILA) and 



for quantum state transfer between light and mechanics 
(Sec.[X|. In the quantum-optical domain, it is referred to 
as a "beam-splitter" interaction. 



For A ^ -\-^rm the dominant terms in RWA 



hg{5d)P + 5ah) 



(31) 



represent a "two-mode squeezing" interaction that lies 
at the heart of parametric amplification ( Clerk et al\ 
2010a). In the absence of dissipation, this would lead 



to an exponential growth of the energies stored both in 
the vibrational mode and the driven optical mode, with 
strong quantum correlations between the two. Thus, 
it may be used for efficiently entangling both modes 
(Sec. [x|. Focussing on the mechanical mode alone, the 
growth of energy can be inte rpreted as "anti-damping" 
or amplification (Sec. V.B.2[ ). If the intrinsic dissipa- 
tion is low enough, this behaviour may trigger a dy- 
namical instability that leads to self-induced mechanical 
oscillations. The resulting features will be discussed in 

Sec. Eml 

Finally, when A = 0, the interaction 



(32) 



means that the mechanical position x ex h ^ P leads 
to a phase shift of the light field, which is the situa- 
tion encountered in optomechanical displacement detec- 
tion (Sec. VI). In addition, this interaction Hamiltonian 



can be viewed as implementing QND detection of the op- 
tical amplitude quadrature 5a + 5a\ since that operator 
commutes with the full Hamiltonian in this case. 



C. Optomechanical equations of motion 

The mechanical motion induces a shift of the optical 
resonance frequency, which in turn results in a change 
of circulating light intensity and, therefore, of the radi- 
ation pressure force acting on the motion. This kind of 
feedback loop is known as optomechanical "backaction". 
The finite cavity decay rate n introduces some retarda- 
tion between the motion and the resulting changes of the 
force, hence the term "dynamical" backaction. 

A convenient starting point for the analytical treat- 
ment of the radiation-pressure dynamical backaction 
phenomena (Sec. V.B and Sec. |VII| ) is the input- 
output formalism. This formalism (briefiy introduced in 
Sec. II. A. 2) provides us with equations of motion for the 



cavity field amplitude a and, analogously, for the me- 
chanical amplitude h. These equations have the form of 
Quantum Langevin equations^, since both the light am- 



plitude and the mechanical motion are driven by noise 
terms that comprise the vacuum noise and any thermal 
noise entering the system: 



a — ^ ^'^^ ^ Gx)a - 



iQrn — ] i> + ig{)d)a 



^finit) (33) 



Tm^in (^) 



(34) 



Please see Sec. |II.A.2| for remarks on the input-output 
treatment and the optical decay rates k,^ k^q^^ k^q. With 
regard to the damping term — rm/2 for the mechanical 
dissipation, we note that this treatment is correct as long 
as ^ Tm- Otherwise the equations would have to be 
formulated on the level of the displacement f , with a 
damping force —mQf^T^x. 

The noise correlators associated with the input fiuctu- 
ations are given by: 

{dUt)ciUt'))=S{t-t') (35) 

{alit)&in{t')) = (36) 

{kn{t)bl{t')) = {nth + mt-t') (37) 

{bl{t)bUt'))=nthS{t-t') (38) 

Here we have assumed that the optical field has zero ther- 
mal occupation (/c^T/ficJcav ~ 0), which is an approxima- 
tion that is valid for optical fields at room temperature, 
although it may fail for the case of microwave fields, un- 
less the setup is cooled to sufficiently low temperatures. 
In contrast, the mechanical degree of freedom is typically 
coupled to a hot environment, with an average number 
of quanta given by nth ^ kBT/hQrn ^ 1- Together with 
these correlators, the quantum Langevin equations de- 
scribe the evolution of the optical cavity ffeld and the 
mechanical oscillator, including all ffuctuation effects. 

It is equally useful to give the classical, averaged ver- 
sion of these equations that will be valid for sufficiently 
large photon and phonon numbers, in the semiclassical 
limit. Then, we can write down the equations for the 
complex light amplitude a{t) = {a{t)) and the oscillator 
position x{t) = {x{t)): 



rrief^X = — meff^rn^ ~ ^effrm^ + \a\ 



(39) 
(40) 



Here we have neglected all ffuct nations, although these 
could be added to describe thermal and even, in a semi- 
classical approximation, quantum noise forces. The term 
ain represents the laser drive. Note that we have also cho- 
sen to write the mechanical equation of motion in terms 

of the displacement, where x = 2xzpFRe ^^6^^ This 



^ In the standard approximation adopted here, these equations are 



Markoffian, i.e. without memory and with delta-correlated noise. 
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becomes equivalent to the equation given above only for 
weak damping, Fm ^ ^^m- These fully nonlinear coupled 
differential equations will be the basis for our discussion 
of nonlinear phenomena, in particular the optomechani- 
cal parametric instability (also called "self- induced oscil- 
lations" or "mechanical lasing". Sec. VIII). 

(and likewise 



The equations of motion Eq. (33), (34 



their classical versions) are inherently nonlinear as they 
contain the product of the mechanical oscillator ampli- 
tude and the cavity field (first line) or the radiation pres- 
sure force (X d)a that is quadratic in photon operators 
(second line). While they can still be solved numerically 
in the classical case, for the quantum regime they are of 
purely formal use and in practice cannot be solved ex- 
actly, neither analytically nor numerically. However, in 
many situations that we will encounter it is permissible 
to linearize this set of equations around some steady- 
state solution: a = a^5a. Using a = ^/n^ and keeping 
only the linear terms, we find the following set of coupled 
linear equations of motion: 



(41) 
(42) 

b + ig {5a + 6a^) + 

(43) 

These correspond to what one would have obtained alter- 
natively by employing the "linearized" coupling Hamilto- 
nian of Eq. ( 28 ) and then applying input-output theory. 




Here we havefas is common practice) redefined the ori- 
gin of the mechanical oscillations to take into account the 
constant displacement hG\a\'^ / rrief^Q'^ that is induced by 
the average radiation pressure force. It is evident that 
now the mutual coupling terms between the optical and 
mechanical degrees of freedom are linear in the field op- 
erators, and that the strength is set by the field-enhanced 
coupling rate g = ^oV^cav • 

As shown in the later sections, these linearized equa- 
tions are able to fully describe several phenomena, in- 
cluding optomechanical cooling, amplification, and para- 
metric normal mode splitting (i.e. strong, coherent cou- 
pling). They can be solved analytically, which is best 
performed in the frequency domain (see Sec. V.B). 

For completeness, we display the linearized quantum 
equations in frequency space: 



-iuj5a[uj] = (^iA - Sa[uj] + ig (b[uj] + (b"^^ [ou]^ 



(44) 



—iujb[uj] = \^—iQrn ^ j H^] + W {6a[uj] + ((5a''')[cj]) 

+ \/r^^mM. (45) 
Here b[uj] = dte^^^h{t) is the Fourier transform of h. 



Note the important relation {^"^^ [^] = •> which 

has to be taken care of while solving the equations. 

It is equally useful to consider the linearized version of 
the classical equations of motion for the light amplitude 
and the displacement, Eqs. (39) and (40): 



5a = (^iA — 6a -\- iGax (46) 

meffX = — meff^^m^ ~ meffFi + hG{a*6a + a6a*) (47) 

Finally, we display them in frequency space, in the form 
that we will employ in Sec. |V.B| 

-iuj6a[uj] = (^iA - 6a[Ljj] + iGax[uj] (48) 

—mef[(jj'^x[ijj] = — meff^m^[^] -\- i(^Taef[T^x[ljj] (49) 
+ hG{a''Sa[uj] + Q;(fo*)[cj]) 

Again, note {Sa''){uj) = Sa{—ujy. 



IV. EXPERIMENTAL REALIZATIONS AND 
OPTOMECHANICAL PARAMETERS 

The increasing availability of high-quality optomechan- 
ical devices, i.e. high-Q mechanical resonators that are 
efficiently coupled to high-Q optical cavities, has been 
driving a plethora of experiments during the last years 
that are successfully demonstrating the working princi- 
ples of cavity optomechanics. We now discuss some of 
the most frequently used architectures. 



A. Optomechanical parameters 

The following table summarizes the relevant optome- 
chanical parameters for some typical current experimen- 
tal implementations. These are: the mechanical res- 
onator frequency Qm and mass m; the fundamental me- 
chanical (phonon) and optical (photon) dissipation rates 
F^ = flm/Qm and hc^ respectively; the "Q • /" product, 
which is a direct measure for the degree of decoupling 
from the thermal environment (specifically, Qm • / = 
Qm ' ^m/27r > kBT/h is the condition for neglecting 
thermal decoherence over one mechanical period); the 
sideband suppression factor k/^}^ that determines the 
ability to realize ground-state cooling (see Sec |VII[ ); and 
finally the bare optomechanical coupling rate go^ which 
corresponds to the cavity frequency shift upon excitation 
of a single phonon. 

Some parameter combinations are of particular rele- 
vance for optomechanical tasks. The following figures 
provide an overview of the state of the art in current 
experiments. The data is compiled from published ex- 
periments. These are labelled as follows: 



1 dCuthbertson et a/ ||1996D, 2 (|Massel et a/.||2QllD, 3 
dRegal et a/.||2QQ8D, 4 ([Rocheleau a/.[|2QlQp, 5 ([Teufel 
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Figure 7 A. gallery illustrating the variety of optomechanical devices, arranged according to mass. Pictures courtesy (from 
top left, down): N. Mavalvala, M. Aspelmeyer, A. Heidmann, D. Bouwmeester, J. Harris, P. Treutlein, T. J. Kippenberg, L 
Favero, M. Lipson, T. J. Kippenberg/E. Weig/J. Kotthaus, H. Tang, K. Vahala/T. Carmon, J. Teufel/K. Lehnert, I. Robert, 
O. Painter, O. Painter, I. Favero/E. Weig/K. Karrai, D. Stamper-Kurn 
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Reference 


Q,n/27r[Hz] 


m [kg] 


r^/27r[Hz] 


Q ■ f [Hz] 


/^/27r[Hz] 




^o/27r[Hz] 


(Murch et al.\ 


20081 


4.2 • 10^ 


1 . 10-22 


1 - 10^ 


1.7-10^ 


6.6 - 10^ 


15.7 


6-10^ 


dChan et a/.',' 


201l|) 


3.9 • 10^ 


3.1 • 10-^^ 


3.9 - 10^ 


3.9 - 10'^ 


5-10^ 


0.13 


9-10^ 


( Teufel et al , 


2011al 


1.1 • 10^ 


4.8-10-^^ 


32 


3.5- 10'2 


2-10^ 


0.02 


2- 102 


(Verhagen et al.\ 


20121 


7.8 • 10^ 


1.9-10-^2 


3.4 - 10^ 


1.8- 10'2 


7.1 - 10^ 


0.09 


3.4 - 10^ 


(Thompson et a/." 2008) 


1.3-10^ 


4-10-^^ 


0.12 


1.5-10'' 


5-10^ 


3.7 


5-10' 


(Kleckner et a/., 2011) 


9.7-10^ 


1.1 • 10-'° 


1.3-10-2 


9-10^ 


4.7 - 10^ 


55 


2.2-10' 


_[Groblacher et al 


^ 2009a) 


9.5-10^ 


1.4-10-'° 


1.4- 102 


6.3 - 10^ 


2-10^ 


0.22 


3.9 


(Arcizet et al. 


2006a|) 


8.14-10^ 


1.9- 10-'^ 


81 


8.1 - 10^ 


1 - 10^ 


1.3 


1.2 


{ Cuthbertson 


et al. 1996) 


318 


1.85 


2.5 - 10-^ 


4.1 - 10'° 


275 


0.9 


1.2-10-^ 



Table II Experimental parameters for a representative sampling of published cavity optomechanics experiments 
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Figure 8 The single-photon optomechanical coupling strength 
^0, vs. cavity decay rate At, for published experiments (see 
main text for references). A high value of go/n would be 
favourable for nonlinear quantum optomechanical experi- 
ments, working with single photons and phonons (Sec. X.F). 
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Figure 9 Mechanical quality factor Qm = ^^m/Fm vs me- 
chanical frequency Qrn = 27r/m for published experiments 
(see main text for references). The dashed lines represent 
Qm • fm = const, for various Q • / values. Note that full 
coherence over one mechanical period 1/fm is obtained for 
Qm • fm = kBT/h, i.e. Qm ' fm ^ ^ ' 10^^ Is a minimum 
requirement for room-temperature quantum optomechanics. 



et al. 2011a), 6 ( Chan et al.\ 2011), 7 (Gavartin et al. b. Suspended mirrors 



field et al.\ 



2009a ) , 14 ( Ding et al. 



2QT1I7 15 dPark and 
Wang! |2QQ9f7^ dSchliesser et a/."["|2QQ8bD , 17 ([Verha- 



gen a/.[|2QT2| , 18 ([Schliesser a/.[ |2009|), 19 ([Arcizet 



2006 



2006 



2011 



2Q11D, Spinkey et al.\ |2QlQ|7rhompsQn et al.\ |2QQ8[ 
9 ([Wilson et a/.[ [2QQ9D, 10 ([Jiang et fl/.| |2QQ9[ ), 11 ([Lin 
et a/.|[2QQ9|7^([Wiederhecker et fl/.|[2QQ9[ ), 1 3 ([Eichen- 



et al. 



et al. 



22 ( [Groblacher et al. 
24 ( Mow-Lowry et al. 
26 ([Grobla cher et a l\ 
(Murch 



20091S 



2008 



et g/ 1 [2QQ6b"|72Q (jFavero et a/.| [20071 ), 2 1 (jGigan et al. 



201 ip , 28 (M urch et al. 
2QQ8P , 30 (Goryachev et al. 



23 (Kleckner et al. 
25 (Kleckner et al. 



2QQ 9a[), 27 ( [Schleier-Smith 
|2QQ8[ ), 29 Prennecke 
20T2F, 31 dVerlot et a/. 



2009). Different symbols indicate the different optome- 



chanical implementations: suspended mirrors (^), opti- 
cal microresonators (♦), photonic crystal cavities (a), 
suspended nanoobjects (t), microwave resonators (•), 
cold atoms (►). 



An obvious way to realize optomechanical interactions 
in a cavity is to suspend one of the cavity's mirrors. The 
mechanical motion changes directly the cavity length and 
hence the frequency response of such a "rubber cavity"^ 

The first experimental implementations of this type 
were a Fabry- Perot cavity with moving mirrors and date 
back to the early attempts of laser interferometeric de- 



tection of gravitational waves (Abramovici et al. 1992). 



While there the purpose of suspending the macroscopic 
cavity mirrors is to achieve acoustic isolation, optome- 
chanical effects, in particular quantum mechanical radia- 



^ for the origin of this terminology see the Acknowledgements sec- 
tion of ( [Bose et a/.]|1997[ >. 
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Figure 10 Single-photon blockade parameter D — g^/ (Q^rn • i^) 
vs sideband resolution Qm/K, for published experiments (see 
main text for references). A single photon induces a cavity 
frequency shift AcJc — D • which results in a blockade effect 
for a subsequent photon for > 1, as discussed in Sec. |X.Fl 
The maximum two-photon suppression scales with [t^/Q^rn) 
and therefore sideband resolution (green shaded area) is an 
additional requirement for successful single-photon blockade. 




10 



10' 10' 10' 10' 

mechanical frequency Q / in [Hz] 



10' 



Figure 11 The single-photon cooperativity, C = 5'o/(^rm), 
vs. mechanical frequency. This quantity is important for as- 
pects such as the streng th of optomechanically induced trans- 
parency (Sec. VII. B. 2). Moreover, if glhca^^r / {i^T 

m'^th) is on 

the order of unity, the state transfer between light and me- 
chanics is faster than the mechanical decoherence rate. Con- 
tour lines indicate at which temperatures this would be true 
even for single photons (ricav = 1)- 



tion pressure fluctuations, eventually pose the fundamen- 



tal limit for its interferometric sensitivity (Caves 1980 



Unruh 1983[ ) (see Section VI. A for a detailed discussion) . 



At the same time this configuration allows to exploit 
cavity optomechanics for the center of mass motion of 
truly macroscopic test masses. Experiments of that type 



have to dat e resulted in the de monst ratio n of the optical 
bistability ( [Dorsel et"a^|1983D (Sec . |V.A[), of the optica l 



spring effect (jCorbitt et al. 
(Sec 



2007a 



Sheard et al. 



V.B.1|) and of optical cooling (Corbitt et al 



2004h 



2007b 



Mow-Lowry et a/. ^2008 ) (Secs. |V.B[|VII| ) with suspended 
mirrors on the gram scale and even of feedback cooling of 
suspended mirrors on the kilogram scale at the LIGO fa- 
cility (Abbott et al. 2009). A practical challenge of these 



experiments is their operation at very low mechanical fre- 
quencies (l]rn/27r <1 kHz), which requires sophisticated 
isolation against acoustic noise. For such experiments 
this is achieved by suspending the macroscopic mirrors 
over several stages. To minimize mechanical losses it has 
recently been sugge sted to levitate the macrosopic mirror 
in an optical trap (Singh et al. 2010). Finally, this kind 



of setup also allows to monitor and to optomechanically 
control internal mechanical modes of macroscopic mirrors 
( [Cohadon et al. 1999 Hadjar et al. 1999), which are a 
dominant source for unwanted cavity phase noise e.g. in 



gravitational wave detectors (Harry et al. 2011 2002), in 



cavity QED (BuckJrJ 2003 ) or in frequency st abilization 
of atomic optical clocks (Numata et al. 2004). 



Another possibility is to use highly reflecting mi- 
cromechanical devices as a Fabry-Perot end mirror. 



These systems include coated cantilevers (Arcizet et al. 
2006a[ |2QQ8l IHohberger Metzger and Karrai[ |2004| |Tit- 



tonen et al. 



1999) 



and micropillars (Verlot et 
201 ID, microm eter-sized mirror pads on top of cantilevers 
Groblacher et al. 



(IFavero et al. 



et al. 2006 



2009b Kleckner 



_ j2007 

201 ID, or micromechanically suspended op- 



tical coatings (Bohm et al. 2006| Cole et al.\ 2008| 2010) 



and photonic crystal slabs (Antoni et al. 2011 



Kemik- 



tarak et al. 2012). Efficient optomechanical coupling 
in this configuration requires the size of the mechanical 
structure to be much larger than the wavelength of the 
light: typical cavity lengths range from 10~^m to 10~^m, 
with an optical finesse up to 10^, which is generally lim- 
ited by losses due to the finite cavity mirror sizes. It 
has been pointed out that additional interference effects 
may be able to overcome this limit significantly (Kleck- 



ner et al. 20To|. Compared to the macroscopic mirrors 
discussed above, these micromechanical devices allow ac- 
cess to higher mechanical frequencies (up to some tens of 
MHz) and, in principle, to higher mechanical quality fac- 
tors. In particular, the possibility of exact geometric con- 
trol via microfabrication techniques allows to minimize 



mechanical losses due to clamping (Anetsberger et 
^OOSllCole et a/.l l20TTl iWilson-Rael |2008D. 



The accessible mass and frequency ranges in combina- 
tion with the restrictions on cavity length (L > A) and 
achievable cavity finesse set some practical limitations for 
this geometry. In particular, sideband-resolution (small 
n) and large optomechanical coupling (small L, lead- 
ing to a large n) impose conflicting conditions and need 
to be traded against each other. On the other hand, 
alo ng wit h the macroscopic microwave transducers (see 
Sec IV.F), this realization provides optomechanical con- 



trol over by far the largest range of mass and frequency. 
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C. Optical microresonators 



D. Waveguides and photonic crystal cavities 



A situation similar to the Fabry-Perot case occurs in 
optical microresonators, where light is guided in whisper- 
ing gallery modes along the rim of a circular resonator 
( jVaSalal |2QQ3D . There is a large number of different 
mechanical normal modes of vibration of these struc- 
tures. The resulting distortions of the structure directly 
modify the optical path length of the resonator, shift- 
ing its optical resonance frequency and hence generat- 
ing optomechanical coupling. The small size of microres- 
onators allows to achieve large coup ling rates go (Ding 

|2012D and to access me- 



et al 



2011 Verhagen et al 



chanical frequencies from several MHz up to some GHz. 
In essence, three different architectures can be distin- 
guished: (i) Microdisk resonators, which are the stan- 
dard resonator structure in planar photonic circuits and 
can be fabricated with high precision. Recent experi- 
ments have demonstrated large opto mechanical couplin g 
rates up to ^ 27r x 8 • 10^ Hz ( Ding et al. 2011 ). 
A fundamental limit in their performance is given by 
radiation losses at the sidewalls. Another limitation is 
due to internal materials losses, which could be improved 
by using single-crystalline materials. A first demonstra- 
tion in this direction, specifically optomechanical cou- 
pling to internal modes of a single-c rystalline CaF2 res - 
onator, has recently been reported (Hofer et al 2009). 



(ii) Microsphere resonators, which allow a larger opti- 



cal q uality (Park and Wang[ [2009 



Tomes and Carmonl 



2009| ); there, the mechanical quality is mainly limited 



by internal materials losses, in particular for the often 
used silica microspheres, (iii) Microtoroidal resonators, 
which are obtained from microdisk resonators by a ther- 
mal reflow process that melts the sidewalls into a toroidal 
topology. The generated smooth surface together with 
the microfabrication control provides a combination of 
high op tical (Armani et al. 2003") and high mechanical 
quality (Anetsberger et al. 2008). This has resulted in 



the first demonstration of radiation-pressure driven op- 



tomechanical parametric amp l ification (Carmon et al. 



2005 Kippenberg et al. 2005 Rokhsari et al.\ 2005) as 



well as of sideband-resolved operation (Schliesser et al. 



[2008b) . Recently, hybrid toroid devices have been devel- 
oped that combine optomechanical and electromechani- 



cal actuation (Lee et al. 2010) 



The practical benefits of these geometries are the avail- 
ability of large optical qualities in combination with the 
resolved sideband regime < essentially owed to the 
fact that the mechanical frequencies range from 10 MHz 
to several GHz. Possible limitations arise from the ne- 
cessity to propagate light inside a solid-state medium, 
which increases the sensitivity to optical absorption and 
thermorefractive cavity noise. 



On-chip waveguides and photonic crystal cavities offer 
a different implementation architecture. 

Photonic crystals are formed by a periodic modulation 
of the index of refraction of some material (typically sil- 
icon), which leads to the formation of optical bands, in 
analogy to the electronic bands for electron waves prop- 
agating in a crystal lattice. Light cannot propagate in 
the band-gaps. Thus, when artificial defects are intro- 
duced into th e periodic patt ern, localized electromagnetic 
field modes (Vahala 2004) can form that do not decay 
into the continuum inside the structure. These struc- 
tures are called photonic crystal cavities. To obtain an 
optomechanical device, in-plane photonic crystal cavities 
are underetched to form nanomechanical beams. The 
mechanical motion results in modulations of the cavity 
boundaries and stresses in the material, both of which 
contribute to the optomechanical coupling between the 
cavity photons and the mechanical modes of t he struc- 
ture. This has been demonstrate d both for ID ([Eichen- 
field et al.\ [2009a|b|) and for 2D ( [Gavartin et al.\ [2011 



Safavi-Naeini et al. 



2010) photonic crystal cavities . The 



small cavity dimensions in combination with the small 
mass of the localized mechanical mode result in an op- 
tomechanical coupling strength that is much larger than 
in regular Fabry-Perot approaches, with current experi- 
ments achieving ^o/27r ^ MHz. The available mechani- 
cal frequencies can range from several tens of MHz up to 
several GHz, which significantly reduces the environment 
thermal occupation hth = kBT/hVt^. The idea of creat- 
ing bandgaps by inducing periodic boundary conditions 
can be extended to the modes of the mechanical beam. 
Introducing a surrounding periodic structure matched to 
the phonon wavelength ("phonon shield") results in a ID 
photonic crystal cavity with co-localized photonic and 
phononic mod es with a significan tly i ncreased mechani- 
cal quality Qy^([Chan et al.\ 2011 ); see (Eichenfield et al. 



2009c Safavi-Naeini and Painter 2010) for more on the 



design of ID and 2D opto mechanical crystals (Safavi- 

It is also possible to in- 



Naeini and Painter 



2010) 



tegrate two-level quantum systems inside the photonic 
crystal nanobeam, e.g. by fabricating photonic crystal 



cavity nanobeams out of diamond (Riedrich-Moller et 
2012) or out of GaAs ( Rundquist and Majumdar{ " 2011[ ) , 



which can include artificial qubits formed by for example 
nitrogen vacany (NV) centers or by quantum dots, re- 
spectively. Another possibility is to have hybrid devices 



with both optical and electrical actuation (Winger et al. 
20TT1). 



Currently, because of the large available coupling rates 
go/ this approach may allow to enter the regime of non- 
linear photon-phonon interactions (see Sec. X.F). More- 
over, the large mechanical frequencies in the GHz range 
could allow for low-temperature operation in a regime 
where the average phonon number drops below one even 
without additional laser-cooling. This would be highly 
beneficial for quantum applications. Finally, the in-plane 
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architecture is immediately compatible with the archi- 
tectures of integrated (silicon) photonics and provides a 
direct route to larger-scale optomechanical arrays, which 
is interesting in the context of classical and quantum in- 
formation processing, and for the study of collective dy- 



namics (Sec. IX). 



It should be noted that optomechanical forces can be- 
come strong even in the absence of a cavity, for structures 
with waveguides running close to a substrate or close to 
each other. This approach (while somewhat outside the 
domain of the concepts covered in the present review) 
could be very fruitful for applications, since it does away 
with the bandwidth restrictions generated by a cavity 



(Bagheri et al. 


2011 Li et al 


2009a|b 2008 Pernice 


et al 2008 


2009). In another equally promising devel- 



opment, the GHz acoustic vibrations of photonic crystal 
fibres are being e xcited and controlled via optomechan- 

2012alb» "Dainese et al 



ical interactions (Butsch et al 



20061 [Kang et a/., ,2010| |2008| |201in2009t ,Wiederhecker 
et a/.l 120081). 



E. Suspended and levitated nano-objects 

This class of cavity optomechanics implementations 
uses a rigid optical cavity that contains a mechanical el- 
ement either inside the cavity or in the near field of the 
cavity. It allows in particular the efficient optomechan- 
ical coupling to sub- wavelength size mechanical objects, 
which has been demonstrated for systems such as high- 
quality mechanical membranes made of high-stress SiN 
fSankey g/.l 1 2010 1 [Thompson et aI[|2008D, stochiomet - 
ric SiN (I Wilson et fl/.||2009D or AlGaAs (|Liu et flr||201lD , 



and for carbon nanowires (Favero et al\ 2009), which 



have been suspended inside state-of-the-art Fabry-Perot 
cavities. The embedded n ano-objects modify the cavity 
field either via dispersion ( Thompson et al\ 2008 ) or via 
dissipation, as suggested in fXuereb et al [2011 ). 



An alternative approach to Fabry-Perot resonators is 
to exploit near-field effects close to the surface of optical 
microresonators, where the evanescent optical field allows 
dispersive coupling to other structures. In essence, the 
mechanical motion modulates the distance d between the 
interfaces. Due to the near-field character the optome- 
chanical coupling strength scales exponentially with d 
and hence allows to generate large values for . This has 
been used to demonstrate optomechanical coupling be- 
tween a toroidal mi crocavity and a nearby Si N nanome- 
chanical resonator (Anetsberger et al 2009). Another 



( Barker and Shneider 2010 Chang et al 2010 Romero- 



Isart et al\ [2010 ). This implementation a llows a direct 
extension to matter- wave interferometry ( Romero-Isart 
et a/.[ |2011) and may enable fundamental tests of quan- 
tum theory in a new macroscopic parameter regime (see 
also Sec. |X.D[ ). The necessary para meter regime for such 
tests is experimentally challenging ( Romero-Isart [ 20111) 
and may even require a space environment ( Kaltenbaek 



et al 



2012). Levitation of micrometer-size ( |Ashkin and 



DziedzicI |197^ 



Li et al 2011) and sub-micrometer size 



( [Gieseler et al 2012) silica spheres has already been 
demonstrated in optical dipole traps in high vacuum. An 
alternative approach could be to combine optical trap- 
ping with a low-frequency mechanical suspension (Cor- 



bitt et al 2007b Ni et al 2012), which has been sug- 



gested to lead to thermal decoupling of s imilar quality 
as pu rely optical trapping (Q • / ~ 10^^) (Chang et c 
20121). 



A prominent feature of such setups, with a nano- 
object inside the standing light wave of a cavity mode, 
is quadratic coupling to position. The optical frequency 
shift may be no longer linear but rather quadratic in the 
mechanical displacement, if the object is placed at a node 
or antinode. This could lead to interesting applications, 
such as QND detection of single phonons, as explained 
in Sec. |VI.B.2[ These setups have also been suggested to 
strongly couple two nano-objects, for example a mechan- 



ical membrane to a sin gle atom (Ham merer et al 
Wallquist erd] |2010D (see Sec jX^^ 



2009b 



F. Microwave resonators 

Analogous to optical cavities, LC circuits form a res- 
onator for electromagnetic radiation in the microwave 
regime, i.e. cjc/27r ~ GHz. The motion of a mechani- 
cal element capacitively coupled to this microwave cavity 
results in a shift of capacitance, and thereby of the LC 
resonance frequency {dC/dx oc duOc/dx). Thus, one ob- 
tains the standard cavity-optomechanical radiation pres- 
sure interaction. The first experiments along this line 
have been performed by Braginsky and co-workers (Bra-| 
ginsky and Manukin[|l967||1977[|Braginsky et a/.[|1970p , 
and later in the context of resonant bar gravitational 
wave detection (|Blair et al 1995[ Cuthbertson et 'al\ 
1996); already back then these works have demonstrated 



both cold damping and optomechanical backaction ef- 
fects such as cooling and parametric amplification. Later, 
in the context of ion-trap physics, cooling of a microme- 



related possibility is to co uple two mechanically vibrat- chanical resonator via an LC circuit was shown (Brown 



ing microdisk resonators (Jiang et al 2009 Lin et al 



2009 Wiederhecker et al 



cavities (Eichenfield et al 



2009D or two photonic cry stal 



2009a Roh et al 2010) via 



their optical near field. 

In order to further suppress mechanical clamping 
losses, it has been suggested to levitate the mechani- 
cal objects either by an additional optical dipole trap 
or in the standing wave trap formed by the cavity field 



et al 2007). With the advent of microfabricated super- 
conducting circuits it has become possible to enter the 
size and frequency regime of nano mechan i cal de vices cou- 
pled to microwave cavities (Regal et al 2008). Typical 



available mechanical frequencies range from some MHz 
to some tens of MHz. In order to resemble a low-entropy 
reservoir of the radiation field, which is of particular im- 
portance for quantum optomechanics (see Sec. ra), the 
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microwave photons need to be kept at cryogenic temper- 
atures. For GHz photons, environment temperatures in 
the mK regime are sufficient, which necessitates opera- 
tion inside a dilution refrigerator. Although the momen- 
tum transfer of microwave photons is several orders of 
magnitude smaller compared to photons at optical fre- 
quencies, the bare optomechanical coupling rates go can 
be made comparable to (or larger than) implementations 



in the optical domain ( Pirkkalainen et al. 2013 Roche- 



leau et al. 2010 Teufel et al. 2011a). The essential idea 



is to have a very small coupling gap and to optimize the 
fraction of the total capacitance that responds to the me- 
chanical motion (see also Fig. [8|. 

A current practical challenge for the microwave 
schemes is the sparse availability of quantum optics tech- 
niques such as the preparation and detection of Fock 
states or of squeezed states of the radiation field. How- 
ever, several recent proof-of-concept experiments have 



demonstrated their availabi l ity in principle (Eichler et al. 



2011||Hofheinz et a/. [ |2009| [Mallet et a/.|[2QTT^ 



As a sidenote, capacitive coupling has also been used 
to couple nanomechanical objects directly to two-level 
quan tum systems, e.g. to a superconducting Cooper-Pair 
box ( LaHaye et al. 2009 ) or to a superconducting phase 
qubit (O'Connell et al. 2010). Note finally that the cou- 
pling need not be capacitive. Recently, it was shown that 
a microwave resonator can also be coupled via dielectric 
gradient for ces to the vibrations of a nanobeam |Faust 
et al. 2012). This makes available a larger range of ma- 



where x|*p^ = Y^fi/(2matom^m) denotes the zero-point 
fluctuations of a single-atom in the trapping potential, 
that is \/N times larger than x^pp of the center-of-mass 
motion of the cloud. Here we have assumed that the ex- 
tent of the cloud is small with respect to the wavelength. 
Sometimes the cloud is actually distributed over several 
lattice sites. Note that the trapping potential could be 
provided by another optical lattice or magnetically. Inci- 
dentally, we note that the same kind of der ivation applies 
for trapped dielectric particles (Sec. IV.E). We also men- 
tion that if the atoms are trapped right at a node or 
antinode, the leading optomechanical coupling is to 
inst ead of x , which leads to different physics (e.g. as in 
Sec.|VI.B.2). 



In another experiment, cavity optomechanics was used 
to cool the motion of a t hermal cloud of Cs atoms tra pped 
inside an optical cavity (Schleier- Smith et al. 2011). Fi- 



nally, density fluctuations in a Bose-Einstein condensate 
of 10^ atoms have been used as the mechanical mode in- 
side a Fabry- Perot cavity (Brennecke et al. 2008). In 



both ultracold cases, due to the strong dispersive atomic 
coupling and the small mass (leading to a large zero-point 
motional amplitude), operation was close to the single- 
photon strong coupling regime, qq/k 1. 

More recently, a setup has been demonstrated that 
couples the motion of a vibrating mirror to the motion 
of atoms trapped in a standing light wav e being reflected 
from that mirror ( Camerer et al.\ 2011), without an op- 
tical cavity. 



terials, which could be beneficial for applications. 



G. Ultracold atoms 



V. BASIC CONSEQUENCES OF THE 
OPTOMECHANICAL INTERACTION 



The ideas of cavity optomechanics have also been im- A. Static phenomena: Optical potential and bistability 



plemented by using clouds of up to 10 atoms. Their 
collective motional dynamics can resemble a single me- 
chanical mode that, for the case of ultracold atoms, is 
already pre-cooled to its quantum ground state of mo- 
tion. In one case, the collective motion of a cloud of ul- 
tracold Rb atoms inside a Fabry Perot cavity was used to 
observe sign atures of shot-no ise radiation pressure fluc- 
tuations (Murch et al. 2008). The dispersive coupling 



of the collective motion of the cloud to an optical cavity 
field results in a position-dependent frequency shift and 
therefore to quantum optomechanical interactions. 

Suppose the single-photon dispersive energy shift of 
a single atom sitting at an antinode of the standing 
light wave pattern is 5E = — fi(^Q*)^/Aat, with the 
atom-cavity vacuum Rabi frequency, and Aat the de- 
tuning between atom and cavity resonance. Then the 
coupling Hamiltonian between the cavity mode and an 
atom cloud of N atoms trapped near position x would 
be NdEa^asm'^{k{x -\- x)). Expanding to lowest order in 
f , this yields a bare optomechanical coupling rate 



We first deal with the simplest case, when the light 
force reacts instantaneously to the mechanical motion. 
This will be relevant for Hi ^ Qrn- Then, the radiation 
pressure force F{x) = hGucavix) depends on the dis- 
placement X via ncav(^)7 the photon number circulating 
inside the optical mode. Such a ID conservative force 
can be derived from a potential (Fig. fl2]): 



F{x) 



dx 



(50) 



For the case of a single, high- finesse optical resonance we 
have ncav(^) = n^^"" /[l^{2{Gx^ A)/ 1^)^], where nj^^^ is 
the maximum number of circulating photons, obtained at 
resonance (proportional to the incoming laser intensity). 
As a result, 



Kad(^) 



-n/^nJ?^t''arctan[2(Gx + A)//^] . (51) 



go 



(kxf^^) sm{2kx)\^, 



Note that for the case of photothermal forces, the dis- 
cussion still applies, only with a different prefactor in 
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F{x) oc ncav(^)- The overall potential for the mechani- 
cal motion also includes the intrinsic harmonic restoring 
potential: 



V{x) 



-X^ + Kad(^) . 



(52) 



The first effect of the radiation force is to shift the equi- 
librium position to xq 7^ 0, with V'{xq) = 0. In addition, 
the effective spring constant is changed to its new value 



(53) 



where the second contribution is called "optical spring". 
In particular for low-frequency mechanical modes, this 
term can be orders of magnitude larger than the in- 
trinsic mechanical spring (Corbitt et al. 2007a). Such 



an approach essentially amounts to a variant of opti- 
cal trapping and can be exploited to diminish the un- 
wanted mechanical dissipation and heating connected 
with the intrinsic mechanical spring being attached to a 
substrate. In the limit of low light intensity, the resulting 
correction for the mechanical frequency (obtained from 
(0) = m,^5{^l^)) is 



8A 



[1 + (2A//t) 
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(54) 



This corresponds to the limit n ^ Vt^n of the dynamical 
case discussed below (note the relation between ncav and 



At larger light intensities V{x) may develop into a 
double-well potential, with two local minima, leading 
to what we may term static bistability. Both of these 
minima correspond to stable equilibrium positions, deter- 
mined by the nonlinear equation F(xi/2) = 
Physically, they represent situations with low/high light 
intensity and low/high restoring force. The bistable be- 
haviour will occur at negative detunings, A < 0, roughly 
when ^ —^m-, such that the original equilibrium 

position becomes unstable. As negative detunings are 
necessary for cooling, the bistability limits the achiev- 
able cooling laser intensities for the case k > 

More quantitatively, analysis of the equation F{x) = 
mQf^Vfix shows that bistability sets in first (at a single 
value of A) when the maximum correction to the spring 
constant, —dF{x)/dx^ obtained at Gx-\-A = — /^/(2V^), 
equals the intrinsic spring constant meff^^. This hap- 
pens at a critical detuning of A = — ^/3/^/2 and at a 
critical laser power determined by 



(55) 



where we employed the relations G = ^o^zpf and x|pp = 
h/{2mei^rtrn)' At higher light intensities, the range of de- 
tunings for which bistability is observed widens. In ex- 
periments, bistability is revealed in hysteresis, e.g. when 



V effective mechanical potential 




bistability 







position X 


^cav 








\ ^ 
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^ circulating intensity 



laser detuning A 

Figure 12 Optomechanical static bistability occurs when the 
laser intensity is sufficiently high to generate two stable lo- 
cal minima in the effective mechanical potential (top). This 
results in bistable behaviour and hysteresis when recording 
the circulating photon number ricav or the transmission as a 
function of laser detuning (bottom). 



recording the transmission or phase shift while sweeping 
the detuning up and down. The first experiments on op- 
tomechanical bistability, with a macroscopic mirror, were 
reported and analyzed already in the 1980s, both in the 



the microwave domain (Gozzini et al, 1985). 



optical (Dorsel et al. 1983 Meystre et al. 1985) and in 



For low- finesse systems, nearby optical resonances may 
also become relevant for the mechanical motion, leading 
to a more complicated effective potential, possibly with 
several local minima (Fig. 12). 



B. Dynamical backaction 

We now turn to dynamical effects, due to the retarded 
nature of the radiation pressure force. 

To derive the dynamics arising from the optomechan- 
ical coupling, one can solve the linearized coupled equa- 
tions of motion for the light and the mechanics, as pre- 
sented in section |III.C[ This is best done in frequency 
space. We employ the classical linearized equations 
Eqs. (48) and ( [49| ) as the basis for our following anal- 
ysis. This is possible since we will be only interested in 
the linear response to an external mechanical force, and 
the averaged linearized quantum equations are identical 
to their classical version (and do not contain the noise 
sources anymore). 

In the absence of optomechanical coupling, the me- 



22 



(D 
O 



o 

ID 
05 




> 

heating/ 
amplification 



red-detuned blue-detuned 



, where we employed the relation G = ^o/^zpf and ob- 
tain hG^ = 2meff^m^^- For now, we just define 



such that the new terms have the structure that is sug- 
gested by the form of the original susceptibility, leading 
to: 



X 



Figure 13 "Thermodynamic" style schematic diagram depict- 
ing the work done by the radiation force during one cycle of 
oscillation. The work is given by the finite area swept in the 
force-displacement diagram, which is due to the retardation 
of the force (finite cavity decay rate) . The work is negative or 
positive, depending on whether one sits on the red-detuned 
or blue-detuned side of the resonance. This then gives rise to 
damping or amplification, respectively. 
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Xxx (^) = ^eff (^m + 2CJ^1]^ {u)-UJ^-iu [Tm + Topt {uj)]) 

(58) 

The real and imaginary parts then yield the frequency- 
dependent mechanical frequency shift Srtm{uj) and op- 
tomechanical damping rate ropt(cj), whose meaning will 
be discussed further in the next section. The ex- 
plicit expressions are (by taking real and imaginary 



parts, 5Qrn{(^) = 
— ImS (cj) /meffCj) : 



Rei;(cj)/2cjmeff, and ropt(c^) = 



ropt(^) =9 



2 

(jj 
2 



(A + cj)2 + A>:V4 (A - cj)2 + A>:V4 



_ (A + cj)2 + /4 (A - a;)2 + tz^ /4 



Figure 14 Schematic optomechanical feedback loop 



chanical oscillator has the susceptibility Xxx oi^) — 
^eff[(^m ~ ~ ^^m^] (see Sec. ILB.S). We wih now 



assume that a weak test force F acts on the mechani- 
cal oscillator in the presence of the optomechanical in- 
teraction. Solving the coupled set of equations, we can 
find the mechanical response to that force. This defines 
the modified mechanical susceptibility, which can be ex- 
pressed in terms of the original susceptibility plus some 
modification^^: 



Xxxi^) = xJoM + ^(^) 



(56) 



The coupled equations (|48|) and (|49|) are solved by ex- 
pressing Sa[uj] in terms of x[uj] and inserting the result 
into the equation for x[ijo]. This yields the modification 
of the linear response to an external force: 



(A + cj) + m/2 (A - uj) 



m/2 
(57) 



These expressions provide an exact solution of the lin- 
earized problem that is also valid in the regime of strong 
coupling, where g > n/2 (Sec. VILB). Also note that the 
effect is linear in the laser drive power, i.e. in the circu- 
lating photon number: g'^ = ^o^cav In the next sections 
we discuss the resulting physical phenomena, i.e. optical 
spring effect and amplification and cooling. 

The frequency dependence of E(cj) will in general yield 
a non-Lorentzian lineshape for the susceptibility that will 
even turn into a d ouble-peak structure at strong cou- 
pling (Sec. VILB.l). However, for sufficiently weak laser 
drive {g <C /^), it is permissible to evaluate dQ^ni^) and 
ropt('^) at the original, unperturbed oscillation frequency 
UJ = Qyn- Then, we just have a shifted and broadened 
mechanical resonance. This picture also explains why 
we need the assumption ^ Teff for this approach to 
hold. A high-Q mechanical oscillator samples the optical 
density of states at cj = il^m, with a small frequency 
linewidth Feff that can be neglected as long as Teff <C hz. 
We will now discuss both quantities, the frequency shift 
and the damping rate, under this assumption. 



1. Optical spring effect 



The notation Ti{oj) reminds us that this modification has the form 
of a "self-energy", in the way that it occurs in the analogous 
Dyson expression for a Green's function of a particle modified 
from the bare Green's function due to interactions. 



We find, with u = (^m, for the frequency shift of the 
oscillator induced by the light field: 
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Figure 15 Optomechanical damping rate and optical spring 
effect vs. detuning and different sideband parameters, with 
theory (black lines) and experimental data (red dots). Left: 
Full mechanical damping rate reff/27r vs. detuning A/k, 
with decreasing values of n/Q^m from top to bottom [n/Q^m — 
3.7, 2.2, 1.4, 0.7). The regime where Feff touches zero is the 
regime of dynamical instability (mechanical lasing). Right: 
Optical spring effect, i.e. l ight-induced mechanical frequenc y 
shift ^r2m/27r. Data from ( Schliesser and Kippenberg| 2010). 
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Figure 16 Optomechanical cooling rate and frequency shift 
in the resolved sideband regime, shown here fo r K/Qm ^ 1- 
(Data courtesy J. Teufel and K. Lehnert, from ( Teufel et al.\ 
2008])) 



In the limit of large cavity decay rate (i.e. the Doppler 
regime, k ^m), this formula yields: 



2A 



This implies that the mechanical oscillator will be spring- 
softenened for a red-detuned laser beam (A < 0), and 
spring-hardened for a blue-detuned laser (A > 0). 

Note that the frequency shift takes a markedly different 
form when entering the resolved sideband regime. Here, 




Figure 17 Scattering picture of cooling and amplication, 
(a) Amplification and heating proceeds by suppressing anti- 
Stokes scattering and enhancing Stokes scattering, via the 
cavity density of states, (b) Cooling proceeds vice versa, 
by suppressing the Stokes process and enhancing ant i- Stokes 
scattering. 



the optical spring effect vanishes at certain detunings and 
the radiation pressure contributes only to cooling or am- 
plification. 



2. Optomechanical damping rate 

Using the same approximation as for the optical spring 
effect, the optomechanical damping rate is given by the 
expression 



/tV4+(A + !^„)2 /t2/4 + (A - i^™)2 

(59) 



This yields the full effective mechanical damping rate: 

Teff = Tm + Topt 

Since Topt can be both positive and negative, it can ei- 
ther increase or decrease the mechanical damping rate, 
i.e. cause extra damping o r an ti-damping. Extra damp- 
ing leads to cooling (Sec. Ill), while anti-damping can 
lead to amplification of thermal fluctuations and finally 
to an instability if the full damping rate becomes nega- 
tive, Teff < (see Sec. |VIIlD . 

This behavior can also be observed experimentally: 
Figure [15] shows the damping rate and optically induced 
frequency shift for different ratios of ^rn/i^- 

The physical origin of the optomechanical damping 
rate can be described in several ways. 

Mechanical picture. Cooling and heating can be un- 
derstood by the following mechanical consideration. As 
the mechanical oscillator is performing its harmonic mo- 
tion, it traces a trajectory in the diagram of radiation 
pressure force vs displacement, as shown in Fig. 13 
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In the limit where k ^ Q^n^ the intracavity field (and 
the associated radiation pressure force) respond instan- 
taneously to the oscillator motion. If, however, we still 
take into account the remaining cavity retardation, this 
leads to the radiation pressure force getting out of phase 
with the mechanical motion. We then can split the force 
F{t) = F + SFi{t) + SFqit) into a component SFi{t) 
in-phase with the motion x{t) (responsible for the opti- 
cal spring effect) and an out-of-phase quadrature term 
SFQ{t) (responsible for cooling or heating). In the dia- 
gram this implies that the oscillator motion traces out 
no longer a line but an area. The sense in which this 
area is encompassed gives the direction of energy flow, 
i.e cooling or amplification. 

Scattering picture. Writing the optomechanical dy- 
namical backaction in the above fashion allows to gain 
physical insight into the origin of cooling and amplifi- 
cation. The expression consists of two terms which are 
essentially the cavity buildup factor ^2 ^^2^4 evaluated 
at the frequencies A = u -\- ftrn and A = u — ftrn- 
These terms describe the strength of the motional side- 
bands of the intracavity field, generated due to the cav- 
ity frequency oscillating because of the motion x{t) — 
xo • sin(l]^t) of the mechanical oscillator. Perturbation 
theory analysis of the classical coupled mode equations 
reveals that the intracavity field consists of sidebands 
( [Kippenberg and Vahala||2QQ7| ): 




Figure 18 Transmission spectrum of an oscillating cavity in 
the resolved sideband regime. The mechanical motion leads 
to the appearance of a series of sidebands spaced by multiples 
of the mechanical frequency, nQm, where n = 0, di, 1, ±2.... In 
optomechanical experiments only the first pair of sidebands 
is relevant. Optomechanical cooling proceeds by pumping 
the lower sideband, which leads to cooling of the mechanical 
oscillator. 



to the absorption spectrum of an oscillating ion ( Sten- 
hStol [19861 ) (Fig.[l8|: 



(A + n(^^)2 + A^74 



(62) 



a(t) ~ ao(^) + <^i(^) 



(60) 



where ao(t) = s • e~*^^^y^/(— + a>:/2) is the unper- 
turbed field and a\(t) contains the Stokes a^it) and anti- 
Stokes aas{t) sidebands: 



Here (3 = denotes the modulation index. For the 

simple case of a weak coherent oscillation (with ampli- 
tude Xq) one obtains to lowest order two sidebands only, 
i.e. J±i(/3)^ = The lower and upper motional side- 
bands that appear in the spectrum are, in analogy to 
trapped-ion physics, related to motional increasing and 
motional de creasing scattering processes (Wineland and 
Itanol 119791) 1^ 



ai{t) = -—ao[t)-{ 



-iQrnt 



JQmt 



i{A + Qm) + 1^/2 -i{A - nrn) + 

(61) 

These two sidebands become asymmetric for nonzero 
laser detuning due to the cavity density of states. Ap- 
plying energy conservation implies that the mode-density 
induced sideband asymmetry extracts or adds power to 
the mechanical oscillator by shifting the frequency of the 
pump photons ^\ i.e. Pl = ■ l/Trt{\as\'^ - Ksl'^) = 
— Fopt • rrief fVt^XQ^ which yields a cooling rate identical 
to the expression derived above (Trt denoting the cavity 
round trip time). Consequently cooling and amplifica- 
tion can be viewed to originate from the imbalance of 
Stokes and anti-Stokes scattering. 

For the case of resolved sidebands (ytrn ^ f^) the cavity 
absorption spectrum 1 — 71{A) develops a series of side- 
bands due to the mechanical oscillator's motion, similar 



Feedback picture. Finally, one can also understand the 
cooling by considering a feedback picture. In this picture 
the mechanical oscillator motion modulates the cavity 
xfield, the latter gives rise to a radiation pressure force, 
which in turn acts back onto the mechanical oscillator. 
Cooling arises again from the phase relationship (retarda- 
tion) between the mechanical motion and the radiation 
pressure force. Note that the feedback is not entirely 
noiseless: quantum noise adds to the intracavity radi- 
ation pressure force, giving rise to a quantum limit of 
cooling treated in section VILA 



It is instructive to consider several limiting cases of the 
cooling rate expression. 

Resolved- sideband regime.- First, in the limit where the 
mechanical frequency is much larger than the cavity de- 



From a quantum mechanical picture each photon that is ab- 
sorbed by the lower sideband removes one quanta of energy from 
the mechanical oscillator (HQrn)- The rate at which photons 
enter the cavity via the lower sideband is given by P^'^/Hujl. 
Consequently the rate at which energy is removed from the os- 



The difference power in the motional sidebands has for this rea- 
son to be multiplied with the pre-factor . 



cillator is: ^ 

expression Fopt = 
classical calculation. 



-^£-(3'^hQm = -ToptE, yielding again the 



'■ ^TT'pQo/'^ agreement with the 
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cay rate {Qrn ^ f^) the cooling rate exhibits pronounced 
maxima and minima near the lower and upper sideband 
(A = The maximum cooling rate is attained on 

the lower sideband (A = —Qrn)' 



- opt \k<^Qm 



9l 



^ (63) 



The cooling rate in the resolved sideband regime can 
also be expressed in a different way, as detailed in 
(Schliesser et al 2008b). In the resolved sideband regime 



(considering an overcoupled single-sided cavity for sim- 
plicity), the relation between intracavity photon number 
and input power is given by ficav = {P/hu) • /^/A^, thus: 



(64) 



As one can see, the cooling rate in the sideband hmit 
does not depend on the cavity Unewidth. 

Doppler regime. - In the unresolved sideband regime 
(Doppler case, i.e. k ^ ^m), the maximum cooling or 
amplification rate is reached for a detuning equal to k/2: 



ropt(A) = g^n, 



-4Ak 



(65) 



ropt(A: 



(k2/4 + A2)2 

/f^m (66) 



Note that in this case the mechanical cooling and am- 
plification rate exhibit a strong dependence on the inverse 
(cubic) cavity decay rate. This shows that the cooling 
rate strongly diminishes in the Doppler regime. 

Note that sub-Doppler cooling in the unresolved side- 
band regime can still be achieved by using pulsed optical 



pumping schemes (Machnes et al 2012 Vanner et al 
20121 [2QlT| [Wang et aL\ \20U^ . 



VI. QUANTUM OPTICAL MEASUREMENTS OF 
MECHANICAL MOTION 

One of the principal advantages of cavity optomechan- 
ical systems is the built-in readout of mechanical motion 
via the light field transmitted through (or reflected from) 
the cavity. In the following, we will discuss several vari- 
ants of optical measurement schemes. We first address 
the measurement of position, where we will find that 
quantum mechanics places fundamental restrictions on 
the overall precision in the regime typically employed in 
experiments (i.e. weak measurements). Then we discuss 
alternative schemes, where there are no such limitations. 
Some supply a measurement of a selected quadrature of 
mechanical motion, i.e. the amplitude of the cos{Qrnt) or 
sm{ftrnt) contribution to x{t). Other schemes measure 
the discrete phonon number. Both approaches are exam- 
ples of so-called quantum non-demolition (QND) mea- 
surements. We close the section by pointing out some 



experimental issues in phase measurements, and by dis- 
cussing feedback cooling based on the possibility of pre- 
cise read-out. 



A. Parametric displacement sensing and the standard 
quantum limit (SQL) 

Measuring the displacement of a mechanical resonator 
via the transmission or reflection phase shift typically 
involves integrating the data over a long time, in order 
to suppress the noise. It is therefore an example of a 
weak, continuous measurement. It turns out that such 
a measurement applied to the coordinate of a harmonic 
oscillator cannot be more precise than what is known as 
the "standard quantum limit". In the following we give a 
brief qualitative and quantitative discussion. More about 



et al 2010a). 
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1. Introduction and qualitative discussion 

The optical cavity represents an interferometer and 
thus allows a direct measurement of the mechanical posi- 
tion via the phase shift of transmitted or reflected light. 
In practice this requires either a homodyne or hetero- 
dyne detector, in which the signal is brought into in- 
terference with a local oscillator that serves as a phase 
reference. The optomechanical cavity frequency pull Gx 
is converted into a phase shift ex Gx/n imparted on 
the photons during their lifetime 1/n inside the cavity 
(assuming slow motion, Vt^ <C /^). If one tries to mea- 
sure this phase shift using N photons passing through 
the cavity (and interfering with a reference beam to read 
out the phase later), then the fundamental uncertainty 
relation between number and phase yields a shot-noise 
limited imprecision of 50 = 1/\/N. It seems that this 
would allow for an arbitrarily precise readout, provided 
one uses a sufficiently large number of photons. Indeed, 
this would be true for an instantaneous readout with a 
very intense flash of light. 

However, in many experiments one rather performs a 
weak measurement: The noisy signal x(t) determined 
from the phase measurement is effectively integrated 
over many oscillation periods to average away the noise 
and get a sufficient signal-to-noise ratio. That this will 
pose a problem can be seen from the general quantum- 
mechanical uncertainty principle which states that it is 
impossible to follow the trajectory x{t) of a particle with 
arbitrary precision (or, to know both position and mo- 
mentum at the same time). It is instructive to see quali- 
tatively how that limitation is enforced in our case. The 
fluctuating radiation pressure force (again, due to the 
photon shot noise) imprints an unavoidable jitter on the 
mechanical motion. Each of the photons imparts a ran- 
dom kick, and their overall effect on the momentum and 
position will grow like ^/N ^ as in a random walk. That 
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Figure 19 Optomechanical systems transduce displacements 
into changes of the optically transmitted (or reflected) phase. 
Upper panel: amplitude response. Lower panel: phase re- 
sponse. 
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Figure 20 Mechanical frequency spectrum for an example 
of an optomechanical system (data from (Schliesser et al 



2008a)). The actual imprecision noise floor (dark data points 
at bottom) is indicated along with the full noise at the stan- 
dard quantum limit. Similar data demonstrating imprecision 
below that at the SQL exists for nanomechanical oscillators 



using a microwave cavity interferometer (Teufel et al.\ [2009 ) 
or nanomechanical oscillators coupled to the near field of an 
optical microresonator. 



be decomposed into quadratures: 



effect is called "backaction noise" and counteracts the 
increase of phase readout precision at large N. Thus, 
there is an optimum at intermediate photon numbers, 
where the sum of the two effects is minimal. 

The quantitative analysis outlined below will be 
phrased in terms of noise spectra, describing the im- 
precision and backaction noise contributions to the over- 
all measurement error. In that context, the appropriate 
question to ask is how large the error is given a certain 
measurement time (which sets the bandwidth over which 
the spectra are to be integrated). In thermal equilibrium, 
the mechanical oscillator's phase and amplitude will fluc- 
tuate on the scale of the damping time l/Tm. Thus, this 
is the longest reasonable measurement time at our dis- 
posal. The outcome of the analysis will be that one can 
determine the trajectory (or rather its two quadratures) 
up to a precision given by the oscillator's zero-point fluc- 
tuations xzpF = (^/2meff^m)~"^^^ during a time l/T^. 
This statement (see Eq. (73)) is known as the "standard 
quantum limit" of displacement detection in a weak mea- 
surement. It is independent of whether the oscillator is 
in its ground state or at high temperatures. 



2. The standard quantum limit 



x{t) = X COs(r^m^) + y Sin(r^m^) 7 



(67) 



where X and Y remain constant during intervals smaller 
than the damping time. Since X = x{0) and Y = 
p(0)/meff^m, the commutator [x^p] = ih carries over to 



X and y, yielding 



X,Y 



ih/mQf[Q^. Heisenberg's 



uncertainty relation thus reads 



AX ■AY> 



2meff^^n 



^ZPF ■ 



(68) 



, and any simultaneous measurement that tries to mea- 
sure both quadratures with equal precision is limited to 
AX = AY = xzPF^^ 

Here we only discuss the most important results. For 
a much more extended recent discussion of the quantum 
limits to weak measurements we refer the reader to ( Clerkl 



et al 2010a). We will consider a single-sided optical 



cavity driven at laser detuning A = (the optimal case) , 
where the number of photons circulating inside the cavity 
is ncav = "^P / (f^J^cav) ' All uoisc spectral densities will 
be symmetrized in frequency (indicated by S). 



Let us first state more formally why there must be a 
standard quantum limit. The oscillator's trajectory can 



Here AX = ({X - (X»2). 
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The noise in the optical phase readout 50 ~ l/yTV 
induces an imprecision in the x measurement, where 
^^imp ^ i<i50/G = n/ {G\/N)^ as discussed above. In- 
serting the photon number N = tN ~ tA^ncav, we find 
((5x^"^P)2 - (A^/ncavG^)^"^. This can be understood 
as the integral of a flat noise spectral density S^^^ ~ 
^/(^cavG^) over the bandwidth set by the measure- 
ment time t. The complete expression for the quantum 
noise limited imprecision noise spectral density reads: 



(^) = 



(69) 



where we also kept the growth of the noise at higher val- 
ues of uj/k. This is due to the fact that the cavity is a 
low-pass filter, suppressing the contribution of motional 
frequencies uj > n io the phase output. On the other 
hand, the phase noise itself is independent of cj, so refer- 
ring it back to the input leads to larger imprecision at 
higher uj. 

At the same time, the backaction noise force has the 



following spectral density (cf. Sec. VILA): 



(70) 



In general, the product of imprecision noise and backac- 
tion forc e noise densities fulfi lls a fundamental inequality, 
see e.g. (Clerk et al. 2010a), a variant of the Heisenberg 



uncertainty relation: 



(71) 



In our particular case, we see that the equality is 
realized, i.e. the cavity displacement detector is as good 
as allowed by quantum mechanics. 

The total noise registered at the detector, expressed in 
terms of x ("referred back to the input"), reads: 

S'^T\^) = S'^icv) + S'^^oj) + SMo^) \x..{oj)\' . (72) 

Here we have added the intrinsic thermal fiuctuation 
spectrum, the imprecision noise, and the effect of the 
backaction noise force on the displacement, calculated 
via the mechanical susceptibility Xxx 

(see Sec.lII.B.Sl). In 



doing so we have assumed a situation where there are no 
cross-correlations between the force noise and the impre- 
cision noise. See (Clerk et a/. ,|2010a) for a more complete 
discussion including the general case. In the following, 
we will denote the sum of the imprecision and backac- 
tion noises as the total added noise, S^^^. Inserting the 
relation (71), and treating Sff as variable (e.g. by tun- 
ing ncav), we can minimize S^"^^. The minimum (at any 
given, fixed frequency) is reached at Sff = ^/{'^IXxxDi 
and this yields 5'^^ (cj) > h\xxx{uj)\ > h\lmXxxH\- By 
using the quantum FDT at T = (see Eq. (16)), we can 



introduce the spectral density of mechanical zero-point 



fiuctuations, S^^ {ou) = h\lmxxx{^)\' We arrive at the 
fundamental inequality 



ZPF 

XX 



(73) 



This is the standard quantum limit (SQL) of weak dis- 
placement detection. The measurement adds at least 
the zero-point noise, on top of the intrinsic fiuctuations. 
Roughly speaking, the effect on the noise looks as if the 
oscillator's energy were increased by h£t^/2^ i.e. half a 
phonon. However, only the backaction contribution re- 
ally corresponds to a physical increase of the oscillator's 
effective temperature, and at the SQL this contribution 
is half of the overall effect, i.e. hVtrn/^-, the other half 
being provided through the imprecision noise. 

If we measure at the mechanical resonance, uj = 
then the added noise of the cavity displacement detector 
is S^t^in^) = S^n^m) = h/{m,^nM' This corre- 
sponds to the zero-point fluctuations, in a measurement 
time t ~ r~^, just as stated in the introduction. Ob- 
viously, if one is limited by laser power, it is better to 
have a high-quality oscillator (small Fm), which boosts 
S'^^^{Qrn) and makes it easier to reach that imprecision 
level. The power to reach the standard quantum limit 
therefore is a natural expression to characterize a trans- 
ducer. It is given by the expression: 



SQL 



1 



(74) 



Both the imprecision noise and the backaction noise 
are shown as a function of laser power (or optomechan- 
ical coupling) in Fig. 21 When referring to "precision 
beyond the standard quantum limit" in this context, one 
wants to emphasize that one can make the imprecision 
noise alone lower than the SQL (which implies the back- 
action noise is already appreciable). This situation has 
been achieved in optomechanical systems for mechanical 
oscillators of nanoscale f Anetsberger et al. 



et a/. I [2009) and microscale dimensions (Schliesser et 



2009 Teufel 



2008arWestphal et al 2012). 



The observed thermal noise at any large temperature 
can also be used to obtain the value of the standard quan- 
tum limit via the relation 



(^) = 2nth(o^) = 2kBT/hiO 

This is a useful expression, since it is independent of the 
calibration of the x measurement. 

In the context of measurements at the SQL, an impor- 
tant step for optomechanical experiments is to observe 
the effects of radiation-pressure shot noise on the me- 
chanical oscillator. T his has now been achieved already 
in cold atom setups ("Brooks et 



2011 Murch et 



2008|), which are conducted routinely at low tempera- 



tures, and where the particularly low effective mass of 
the atomic cloud leads to a very strong single-photon 
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Figure 22 Observation of quantum radiation pressure force 
fluctuations through the energy transferred to a near-ground- 
state mechanical oscillator. Experiments were performed with 
an ultracold atomic gas serving as the mechanical element 
within a Fabry-Perot optical cavity. In (a), the energy trans- 
ferred to the gas was quantified via the rate at which atoms 
were ejected from a finite-depth optical trap. The force fluctu- 
ation spectral density at the mechanical oscillation frequency, 
*S'FF(^m), is thereby obtained at different detunings A be- 
tween the cavity probe and resonance frequencies. In (b), 
from the power difference between the red and blue motional 
sidebands observed in the emission of a resonantly driven op- 
tical cavity, one obtains the heat flux into the mechanical 
system via the cavity probe. The observed heating, given in 
units of mechanical energy quanta per second, matches well 
to that predicted for intracavity shot noise from a coherent 
opti cal field (gray hne). From (a) ( |Murch et al\ |20Q8| ) and 
(b) (Brahms et al. 2012| ) [Courtesy of D. Stamper-Kurn] 



Figure 21 The full measured noise spectrum contains con- 
tributions from the intrinsic fluctuations of the mechanical 
oscillator, but also extra noise due to imprecision in the mea- 
surement (typically flat in frequency) and noise due to the 
backaction heating of the oscillator. Bottom: The value of the 
measured noise spectral density evaluated at the mechanical 
resonance, plotted as a function of the power of the measure- 
ment beam. At lower powers, imprecision noise dominates 
(few photons yield bad phase resolution), while at higher 
powers the backaction noise represents the most important 
contribution. The standard quantum limit minimal noise is 
reached at intermediate powers. 



has been demonstrated for a model situation with delib- 
erately introduced classical noise instead of the quantum 
shot noise of a laser beam (Verlot et al 2009). In an- 



other experiment, it was demonstrated how the radia- 
tion pressure backaction can be employed for amplifying 
an interferometric signal, which can lead to a sensitivity 



lower than the SQL (Verlot et al. 2010). Recently, first 



signatures of a direct observat ion of radiation pr essure 
shot noise have been reported (Purdy et al 2013). 



coupling rate go. This is of advantage, since the ratio of 
quantum backaction to thermal force noise (at A = 0) is 
given by: 



SM^m)/S'M^m) = go^m/i^kBTT^) (P/^cav) {^^W 

(75) 

These atomic cloud experiments have allowed to access 
the radiation pressure shot noise spectrum, e.g. via track- 
ing the heating of the cloud (see Figure 22). 

Current solid-state based devices still exhibit both sig- 
nificantly smaller ratios go/hz and deleterious effects of 
thermal noise, which make the observation of radiation 
pressure shot noise effects a challenging task. One possi- 
ble strategy is to measure the cross-correlations between 
a strong beam exerting radiation pressure force fluctu- 
ations and another beam measuring the resulting dis- 
placement fluctuati ons (([Heidmann et a/.[|1997[ ); see also 
(Borkje et al. 2010 ) for a more recent analysis). This idea 



B. Optical QND measurements 

The weak displacement measurements discussed in the 
last section effectively try to measure non-commuting 
observables simultaneously, namely the two quadrature 
/ 4)components of motion. They are therefore limited funda- 
mentally by the Heisenberg uncertainty principle. How- 
ever, it is also possible to perform measurements of a sin- 
gle observable. This observable can be measured with ar- 
bitrary precision, thereby approaching an idealized pro- 
jection measurement. Repeating the measurement before 
the state had a chance to decay would reproduce the same 
outcome. This is because the system's Hamiltonian com- 
mutes with the observable (neglecting decay). Therefore 
such me asurements are called quantum non - demolition 
(QND) dBraginsky and Khalili| [l996| p95| [Braginsky 



et al. 1980). These have been successfully realized in the 



quantum optical domain (Haroche and Raimond[ 2006 
Leibfried et al. 1996 Lvovsky and Raymer[ 2QQ9|. 
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1. Single quadrature measurements 

It is possible to optically measure only one quadrature 
component of the mechanical motion to arbitrary preci- 
sion ( Clerk et aL \ [200 8 ). This is important, since it can 
be used for a full reconstruction of the mechanical quan- 
tum state, extracting its Wigner density using quantum 
state tomography (see below). In addition, it can serve to 
measure the correlation of mechanical quadratures with 
either the light field quadratures or those of another me- 
chanical object. The resulting correlators can then be 
used to test for entanglement. The fundamentally lim- 
ited precision of a standard displacement measurement 
would not be sufficient for such tests. 

One way of achieving this makes use of a simple prop- 
erty of harmonic motion: Any force applied at time t does 
not affect the position a full period later, at t + 27T/Qm. 
Thus, the unavoidable perturbation of the momentum 
connected with a position measurement does not destroy 
the precision of "stroboscopic" periodic observations at 
times t -\- n27T/Qrn- An equivalent but more practical ap- 
proach is to do a displacement measurement with a laser 
beam (at detuning A = 0) whose intensity is modulated 
at frequency Q^- This reads out only one quadrature 
(say, X^) in the decomposition 
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Figure 23 (a) Schematic setup for an optomechanical single- 
quadrature measurement, with a laser beam amplitude- 
modulated at the mechanical frequency, (b) In optical fre- 
quency space, the beam carries two sidebands, which are then 
modulated by the motion, (c) Tomographical reconstruction 
of the full mechanical Wigner phase space density is possible 
based on measurements of the quadratures at different phases 



x{t) = cos{Qrnt (f) ^ % sin(^]^t + (f) , (76) 

whose phase cp is determined by the phase of the laser 
amplitude modulation. The backaction noise exclu- 



sively affects the other quadrature (Clerk et al 2008). 



All this can be derived in the Hamiltonian formulation 



(Sec. III.B). Suppose the modulated drive yields an in- 
tracavity amplitude a = 2aocos{Qrnt + 9^), in a frame 
rotating at the cavity resonance. Then the standard lin- 
earization gives 



Hint = —2hGao cos{Qrnt-\-(p){Sa-\-Sd^)x ~ —hGao{Sa-\-Sa^ 

(77) 

where we have omitted rapidly oscillating terms and 
adopted the rotating frame for the mechanical resonator 
in the second step. This is a QND Hamiltonian for mea- 
suring X^p. Deviations from this idealized picture are 
small in the resolved sideband limit, i.e. for n/Vtjn <^ 1- 
As a side-effect, the state conditioned on the measure- 
ment result becomes squeezed, since the variance of X^ 
after the measurement approaches zero. First experimen- 



oscillations has advanced. This approach has been pro- 



posed and analyzed in (Braginsky and Khalili 1995 Bra- 



ginsky et al. 1978 Vanner et al. 2011). Because of the 



short pulses, this scheme operates in the non-resolved 
sideband regime n ^ Vtrn- First experiments in this di- 



rection have recently been reported in (Vanner et al. 
2012). It has also been pointed out that properly con- 
trolled pulse sequences can alter the effective optomome- 
chanical interaction, enabling for example sub-Doppler 
cooling rates in the non-resolved si deband regime ( ^Cer-| 



rillo et al. 



2011 Wang et al. 2011). 



/^"^^ 'In both of these approaches, quantum state tomog- 



raphy ( Lvovsky and Raymer , 2009} then would work by 
repeatedly preparing the same mechanical state, measur- 
ing the probability densities of the quadratures at a large 
number of phases and applying the inverse Radon 
transform to obtain the Wigner density. 

Reconstruction of the Wigner density of quantum 
states of vibrational motion has, for example, success- 



fully been achieved for ions (Leibfried et al. 1996) 



tal steps along these lines have been taken in ( jHertzberg 
et a/.||2009D . 



There is an alternative approach for QND detection 
of single quadratures: measuring much faster than the 
oscillation period of the mechanical resonator. Such an 
approach requires very short, intense laser pulses that 
essentially implement an instantaneous projection mea- 
surement of displacement. Picking another quadrature 
to measure then simply involves performing the measure- 
ment at another time, when the phase of the harmonic 



2. Mechanical Fock state detection 

Another crucial observable in the mechanical oscillator 
is the phonon number h = b^b. It is especially important, 
since measuring the discrete Fock states n = 0,l,2,...is 
a very direct proof of the quantum nature of an oscil- 
lator. We note that the first measurements of phonon 
number in a fabricated, mesoscopic mechanical oscilla- 
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tor have recently been performed in a nanoelectrome- 
chanical system, exploiting the strong interaction be- 
tween a piezomechanical vibration and a superconducting 
qubit (O'Connell et al. 2Q1Q| . Here we discuss a route 
towards observing quantum jumps between mechanical 



Fock states in an optomechanical system ( Jayich et al 



2QQ8| [Thompson et al\ |2QQ8 ). 

The idea is that any measurement of instead of x 
will be closely connected to the oscillator's energy, and 
thus the phonon number. This then permits QND de- 
tection of the phonon number in a mechanical resonator 
(for an early detailed analysis of this concept in an an- 
harmonic two-mode nanomechanical system, see (Santa- 



more et 



2004)). In practice, for an optomechanical 



system this requires changing the standard setup to an- 
other variant where the light field couples to the square of 
the displacement. That can be achieved for example by 
placing a thin membrane (or any other dielectric object) 
inside an optical cavity and positioning it near a node (or 
antinode) of the standing light wave that forms one opti- 
cal mode. In that case, the light intensity at the object's 
position, and thereby the object's effect on the optical 
resonance frequency, depends quadratically on the dis- 
placement X. The optomechanical coupling (sometimes 
termed "dispersive" in this case) thus reads 



where ^q^^ 



(78) 



fiecting, an alternative way of deriving Eq. (|78|) is to focus 



on two optical modes (to the left and right of the mem- 
brane). Their coupling, provided by photon tunneling 
through the membrane, then leads to an avoided level 
crossing in the optical spectrum as a function of x. At 
the degeneracy point, both the upper and lower opti- 



cal branch frequencies are extremal (Sankey et al. 2010 



Thompson et al. 2008), leading to Eq. (78) for each of 



them. Quadratic couplings have recently also b een ob- 
served in a cold atom setup ( Purdy et al.\ 2010| ). 

The phase shift observed in refiection from the cavity 
will measure {b-\-b^)'^ cx f ^, but averaged over the photon 
lifetime l/^c. If Hi <^ Qrn^ the time-average of this term 
reduces to 2pb -\- 1. Alternatively, one may keep only 
this contribution in Eq. ( [78| ), invoking the rotating wave 
approximation. Thus, this setup in principle allows a 
QND measurement of the phonon number. Cooling the 
oscillator to near its ground state and then measuring 
the phase shift vs. time, one should be able to observe 
quantum jumps between mechanical Fock states, similar 
to the jumps t hat have been observ ed between photon 

number states (Guerlin et al. 2007). To this end, the 

TT) '. 

phase shift ~ do /i^ induced by a single phonon must 

be resolved during the lifetime (F^nth)"^ of the ground 
state, by refiecting a sufficient number of photons from 
the cavity. Note that one must suppress any deviations 
from the ideal case, such as imperfect placement away 
from the degeneracy point (Jayic h et al.\ |2008[ |Thomp-| 
[son et al.\ |2008j ) or deviations from a perfectly one-sided 



cavity. Such deviations lead to extra noise acting on the 
membrane that hampers the phonon QND measurement. 
A careful analysis ( |Miao et al. 2009) shows that this 
requires go > a^o, where hcq includes any absorption in 
the membrane or cavity, but not the "good" decay chan- 
nel through the entrance mirror. In practice, for setups 
where n ^ this is identical to the single-photon strong 
coupling regime, to be discussed further in Sec. |X.F[ 

Detailed quantum jump trajectory simulations of the 
phonon detection pr ocess have been first presented in 
(Gangat et al.\ 2011), for the idealized case of a single 



optical mod e. More recently, stu dies of the complete two- 
mode setup ( [Ludwig et al. 2012[) have been able to reveal 
the physics discussed above ( Miao et al. 2009 ) . 



The detection of quantum jumps between mechanical 
Fock states is obviously extremely challenging. However, 
one might detect quantum signatures of the mechanical 
oscillator in another way: When driving it strongly into 
a coherent mechanical state of mean phonon number n, 
there will be Poissonian phonon shot noise of size ^/fl. A 
setup like the one described above may then detect this 
noise, which reveals t he granularity of mechanical quanta 
(Clerk et al. 2010b), much in the same way that mea- 



surements of electrical current noise can reveal the charge 
of the charge carriers (|Blanter and Buttiker( [200^ . Fur- 
thermore, higher moments of the phonon shot noise may 
be observed in this way as well. These display distinctly 
nonclassical features - for an in-depth analysis see ( Clerk 



If the membrane is highly re- 2011). 



3. Optical feedback cooling (cold damping) 

The high sensitivity provided by the cavity readout 
of mechanical motion can also be used for directly cool- 
ing the mechanical motion via active feedback. The 
main idea is to obtain the oscillator position by a phase- 
sensitive detection of the cavity output and to use it to 
generate a negative feedback on the oscillator, i.e. a 
force F = —mef^STx proportional to the time deriva- 
tive of the output signal. This increases the damping 
rate of the system by ST without increasing the thermal 
noise (cold damping). The scheme has been suggested in 
(Mancini et al. 1998| ) and has been exp erimentally re- 



alized in several optomechanical devices (Arcizet et al. 



2006bl [Cohadon et al. , 1999; [Kreckner and Bouwmeester 



,2006, |Poggio et al. |2007J with radiation pressure as the 
feedback force. Because the scheme relies on the pre- 
cise readout of the instantaneous oscillator position, the 
ideal configuration comprises both weak coupling and a 
fast cavity decay, i.e. hz ftrn ^ 9 (adiabatic regime). 
The quantum limits of this cooling method have been 



tah et al. 2002 



discussed in (Courty et al. 2001 Genes et al, 2008 



Vi- 



and it has been shown that ground-state 
cooling by cold damping is possibl e. A detailed discus- 
sion can be found for example in (Genes et al. [2008 ). 
The maximum amount of cooling is limited by the im- 
precision of the read-out. An important aspect here is the 
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Figure 24 (a) The "membrane-in-the-middle" setup used to 
generate an optomechanical couphng to . It can be viewed 
as two optical modes to the left and right of the membrane, 
with photon transmission through the membrane coupling 
those modes, (b) Resulting optical spectrum, with an avoided 
crossing between the left and right optical mode as a function 
of displacement x. Coupling to x^ is obtained at the degen- 
eracy point, (c) Future dispersive setups along these lines 
(with greatly improved parameters) may enable detection of 
quantum jumps between mechanical Fock states in the phase 
shift of light reflected from the cavity. An idealized quan- 
tum jump trajectory (without measurement noise) is depicted 
here, starting from the mechanical ground state (0 phonons). 



phenomenon of "noise squashing", where the noise on the 
detector and the noise-driven motion of the mechanical 



oscillator become correlated (see Fig. 25). 

In the context of cold atoms, feedback cooling of a 
single neutral atom has been implemented in a cavity 
QED setup (Koch et a/.||2QlQD- 



VII. OPTOMECHANICAL COOLING 

Optomechanical quantum control requires the mechan- 
ical oscillator to be in or near its quantum ground state. 
Unless the mechanical frequency Qrn is in the GHz range 
(whic h is true on l y for some recent nanomechanical se- 
tups dChan et a/.} I^QHl fEichenfield et a/.^ 2QQ9b} ), even 



dilution refrigerator temperatures of 20mK are not suffi- 
cient to ensure ksT <^ hQ^n- Thus, additional cooling of 
the selected mechanical mode is needed. 



In a previous section (Sec. V.B), we have already dis- 
cussed dynamical backaction effects and the resulting op- 
tomechanical damping rate. It is obvious that this can be 
used for cooling the mechanical motion. The purpose of 
the present section then is to develop the quantum the- 
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Figure 25 Optomechanical feedback cooling of a cantilever. 
Note the 'noise squashing' at strong feedback (reduction of 
the spectrum below the noise floor due to interference between 
the classical noise in the reado ut and the signal caused by the 
mechanical motion) (Poggio et al. 2007). 



ory of cooling, which in particular describes the limits for 
cooling that cannot be obtained from a discussion of the 
damping rate alone. We will focus in this section on in- 
trinsic cavity cooling. For a comparison of this approach 



to feedback cooling (discussed above), see (Genes et al. 
2008). 

A simple classical theory of an oscillator at initial tem- 
perature Tinit subject to extra damping Fopt (propor- 
tional to the laser power) would predict that its tem- 
perature is reduced down to 



^final 



- opt 



(79) 



However, this classical expression ceases to be valid at 
sufficiently low Tfinai, when the fluctuations of the radi- 
ation pressure force due to photon shot noise set a lower 
bound to the achievable temperature. In the following, 
we discuss the full quantum theory that permits to cal- 
culate the quantum limits to cooling and predicts that, 
in many cases, ground state cooling is possible only in 
the resolved sideband regime n <C l^m- 

It should be noted that sideband resolution does not 
intrinsically prohibit optomechanical ground state cool- 
ing. For example, the displacement may not only couple 
to the cavity frequency, c^cav = ^cav(^), but also to the 
decay rate, yielding n = n{x). This mechanism is quite 
plausible both for nano-objects that scatter light out of 
the cavity (depending on where they sit in the standing 
light wave pattern), and for vibrating objects that are 
evanescent ly coupled to some tapered fibre, as well as for 
other geometries. This kind of behaviour has been ob- 



served experimentally (Li et al. 2009c), and it has been 



predicted (Elste et al. 2009[ ) that this can yield novel be- 



haviour for cooling, with the possibility of reaching the 
quantum ground state even when > Cl^. Similar re- 



sults hold for pulsed cooling schemes ( Cerrillo et al. 2011 
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Vanner et al. 2011 Wang et al. 2011) and photothermal 



|1,n+1> 



cooling. 

Another interesting option is to consider a nonlinear 
medium inside a cavity, which may enhance the efficiency 
of cooling (Huang and Agarwal 2009a). 

Finally, one could do away with the cavity by exploit- 
ing other sharp changes of radiation forces as a function 
of wave length. This has b een proposed e.g. for a Bragg 
mirror ( Karrai et al. 2008 ) or for microspheres with their 



narro w internal whispering gallery resonances (Barker 
20101). 



We will focus on radiation pressure cooling in cavity 
setups, which is conceptually the simplest case. Recently 
it has been argued that photothermal forces, which have 
been exploited for optomechanical cooling early on (Hoh- 
berger Metzger and Karrai^ 2004 Metzger et al. 2008a ) , 



could in principle also lead to the quantum ground state 
(Pinard and Dantan 2008 Restrepo a/. 2011). 
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Figure 26 Principle of cavity optomechanical cooling and cor- 
responding transition diagram. Three processes can be distin- 
guished. First a photon can be absorbed by the cavity leaving 
the motional states unchanged (carrier transitions). Second, 
the photon can create a phonon (corresponding to absorption 
on the upper motional sideband). Third, the photon leads to 
the annihilation of a phonon (corresponding to a photon be- 
ing absorbed by the lower motional sideband). Appropriate 
red-detuning of the laser can lead to effective suppression of 
the first two processes and efficient cooling. 



A. Quantum theory of radiation pressure cooling 

In the following, we will work in the weak coupling 
regime, ^ <C where a perturbative picture a pplies. 
The quantum theory of optomechani cal cooling ([Genes 
eTalj [20081 [Marquardt TTal] [2QQ7| [Wilson-Rae et al.\ 
2007|) is related to earlier approaches fo r trapped ions 
(Itano et al. 1992 Neuhauser et al. 1978) and for cavity- 



assisted laser cooling o f atomic and molecular motion 
dHechenblaikner et al.\ [1998| [Vuletic and Chu[ [2000 ). 
The idea is best explained in a Raman-scattering picture. 
Photons impinging at a frequency red-detuned from the 
cavity resonance will, via the optomechanical interaction, 
preferentially scatter upwards in energy in order to en- 
ter the cavity resonance, absorbing a phonon from the 
oscillator in the process. As a consequence, they will 
be reflected blue-shifted by l^m, carrying away a quan- 
tum of mechanical energy. These "anti- Stokes" processes 
happen at a rate A~ (to be calculated below). More pre- 
cisely, the transition rate from phonon state n to n — 1 
will include a bosonic factor n. 



nA' 



(80) 



The "Stokes" process, where photons return red-shifted 
and leave behind an extra phonon, happens at a smaller 
rate ^4+ (with F^^^+i = (n + 1)^^) due to the suppres- 
sion in the final density of available photon states, if the 
laser is red-detuned. Note that we choose to follow the 
notation for these rates, as it is used in the context 
of atomic laser cooling. 

Given these rates (^4+ for upward transitions in the 
mechanical oscillator, A~ for downward transitions), the 
full optomechanical damping rate is the net downward 
rate. 



opt 



A+. 



(81) 



The average phonon number n = X^^o (with the 
Fock state populations P^) changes according to the rates 
Tn^nii, leading to: 



n 



(n+l)(A++0-n(A-+A-J. (82) 



Here we have introduced the extra transition rates due to 
the oscillator's thermal environment, which has a mean 
phonon number nth* = ^thrm and A~^ = (nth + 

l)rm. In the absence of optomechanical effects, these 
would establish equilibrium at n = nth- Now, however, 
the steady-state final phonon number (requiring n = in 
([821)) win be 



(83) 



opt 



Even in the optimal case, i.e. in the absence of any cou- 
pling to a mechanical thermal environment (Fm = 0), 
this leads to a minimal phonon number of 



opt 



A+ 



A+ ' 



(84) 



The rates can be calculated using Fermi's Golden 
Rule, applied to the interaction of the oscillator with 
the fluctuating radiation pressure force: H±nt = —xF^ 
where F = hGa)a according to Eq. (19). In the 



weak coupling limit, all the transition rates can be cal- 
culated once the quantum noise spectrum Sff{^) = 



dte'^UF{t)F{0)) of the force is known (see e.g. 



( [Clerk et al.\ 2010a) for more on quantum noise spec- 
tra) . Re-expressing the result in terms of noise spectra 
yields: 



(85) 
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where we introduced the photon number noise spectrum 
SnnH = /^^(ite^^*((ata)(t)(ata)(0)>. Exploiting 
shifte d photon operators, one can show ( [Marquardt et al. 



2QQ7| ) that the photon number spectrum of a laser-driven 
cavity is 



/^74 + (A + cj)2 



(86) 



Inserting this into Eq. (85), one obtains and from 
there the optomechanical damping rate Fopt = A --A+. 
which coincides with the expression obtaine d usi ng th e 
hnearized equations of motion (Sec. V.B.2, Eq. ( |59[ )). 
We obtain the final minimum phonon number (84) as 



A- 
A+ 



{K/2f + (A - ^mf 
(«/2)2 + (A + 



(87) 

Experimentally, one can still vary the laser detuning A 
to minimize this expression. In the resolved sideband 
regime n <C l^m, this leads to 
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Figure 27 Quantum noise spectrum Sff{oo) for the radia- 
tion pressure force fluctuations acting on the end-mirror of 
a laser-driven optical cavity. The spectrum is asymmetric in 
frequency, due to the quantum nature of the fluctuations. For 
red detuning, the bulk of the spectrum sits at positive frequen- 
cies. This implies (see main text) that the fluctuations induce 
more downward (cooling) transitions than upward transitions. 



40. 



< 1, 



which permits ground state cooling, while in the opposite 
limit (k » VLra) we find 



> 1. 



(89) 



These two results are almost identical to the case of 
atomic laser cooling (in the atomic laser cooling expres- 
sions, a scalar prefactor is different due to the directional 
dependence of the spontaneous emission). 

In the presence of a thermal environment, the final 
phonon number (83) can be written as the result of cou- 
pling to two baths at average occupations n^i^ and nth 
with coupling rates Fopt and Fm, respectively: 



Uf 



Fopt^min H~ FjjjTT/^/j, 



(90) 



opt 



Finite thermal cavity occupancy. This formula is mod- 
ified when the cooling field has thermal occupation. This 
is in particular the case for microwave fields due to their 



th 



low frequency. If the cavity occupation is given by fi^ 
then the final occupation is modified to (in the resolved 
sideband limit): 



rif nth- 



-th 



- eff 



(91) 



This implies that the final phonon number can never be 
below the effective thermal occupation of the drive field 
( [Dobrindt et al\ [2008 ). It should be noted that when 
the radiation field and the mechanical oscillator initially 
have the same bath temperature (as will be the case in 
equilibrium, without extra absorption), the equilibration 
of these two oscillators of frequency cjcav and Vt^n will lead 
to an effective cooling of the lower frequency mechanical 
oscillator, as nth = ^'lo'"*'' > ^cav = 

Equations- of -motion approach. An alternative way to 
calculate the final phonon number is to write down the 
linearized equations of motion for the oscillator and light 
field, eliminate the light field dynamics, and exploit the 
known spectra of the quantum Langevin f orces d riving 
both the mechanics and the cavity (see Sec. III.C)^^. 



The similarity to atomic Doppler cooling is due to the fact that 
the Doppler cooling limit is a result of the Heisenberg uncertainty 
principle. Each photon that leaves the cavity has an energy un- 
certainty given by its photon decay time, i.e. AE = h/ At, where 
in the cavity cooling case At is the inverse cavity decay rate . 
The minimum phonon number can therefore never be smaller 
than this value, i.e. AE = hn, or expressed as an occupancy 

■^min ^/^m- 



^ Note that there can also be experimental situations where the 
mechanical oscillator mediates scattering of photons out of the 
cavity, such as in the case of levitated particles in an optical 
cavity. Denoting the photon scattering rate caused by the me- 
chanical oscllator as "yscat this additional heating causes a mod- 
ification of the quantum limit of cooling. In the sideband limit 
this yields: fif = ^^ptn^ir^+r^nth ^ jscat ^ The analytical form 

•' i opt +i m 1^ 

of the scattering rate depends on the geometry of the mechani- 
cal oscillator under consideration. For the case of a small sphere. 
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We do this by considering the influence of the quantum 
backaction force noise driving the mechanical oscillator. 
Adopting the classical equations of motion for the me- 
chanical oscillator (but keeping the quantum part of the 
symmetrized noise spectrum), we flnd an average energy 
of: 



{Err 



-CO 



duj |Xeff(^)|^ SFF{^)meffUj'^ 



where Sff is the symmetrized version of the spectrum 
introduced above and we assumed the laser to be domi- 
nant, Feff ^ Fm, and also Q^n ^ Teff. This leads to the 
correct quantum result: 



{Err 



1 



A- 



This expression reveals that the energy, consisting of the 
zero point motion and an additional residual term, is 
caused by the quantum fluctuations of the laser fleld. 
The fact that the quantum fluctuations of the laser fleld 
also give rise to the zero point energy can be under- 
stood by noting that the dominant bath that the me- 
chanical oscillator is coupled to is the laser fleld itself 
when Fg// ^ F^. In this sense then, the correct quan- 
tum limit for cooling can be formulated on a semiclassical 
basis, where it is suflicient to consider only the quantum 
fluctuations of the drive fleld. 

Optical output spectrum. The output spectrum of the 
light that emerges from the optomechanical cavity can b e 
calculated using the input-output formalism (Sec. IIFC). 
In close analogy to the fluorescence spectrum of a laser- 
cooled trapped ion, the specrum is given by the expres- 
sion: 

S{uj) ^IqS {uj - ujl) 

+ A-n//side(c^-^m) (92) 

+ (n/ + l) ^side(^ + ^m) (93) 

This spectrum refers to the frequency-resolved photon 
flux. Here 



^side 



(c) 



Teff 



27r 



{uj-ujl) + 



eff 

4 



is the Lorentzian shape of the sidebands, where the width 
of the sidebands is given by the overall effective mechan- 
ical damping rate Feff. The weight of the main peak^^ 



trapped within a Fabry Perot cavity, an analytic solution of 'jsc at 
can be derived ( [Chang et aL[ |2010[ [Romero-Isart et aL||2010| ). 

It is interesting to consider tne total (integrated) power m the 
generated sidebands. The weight of the sideband for the case 
where the final occupancy is given by hf = 1 is given by 

Psideband/Pcarrier ~ ^ and is thus independent of the ini- 
tial occupancy (or temperature). 




(B) 



Lower sideband 




-2 -1 -1 

(o) - CO J/y 
Upper sideband 

A = +Q, 




Figure 28 (a) Sideband thermometry revealing the growing 
asymmetry between photons scattered into the cavity reso- 
nance when increasing the cooling laser power (Figure cour- 
tesy O. Painter, from ( Safavi-Naeini 2012)) (b) Scheme of 
the readout laser for measuring the photons scattered into 
the cavity resonance. 



at the laser frequency is (for a single-sided Fabry Perot 
resonator) 



A2 



1 ^2 



As expected, the spectrum consists (in addition to 
the carrier) of blueshifted (anti-Stokes) photons at = 
ujl -\-^in (second line of Eq. (93)) and redshifted (Stokes) 
photons at cj = ujL — ^m (third line). The sideband asym- 
metry changes as a function of cooling laser power. De- 
tailed balance causes the initially asymmetric sidebands 
(as A~ ^ A+) to become progressively more symmetric, 
with A+(nj + l) = A~nf in the limit where n/ is entirely 
determined by the cooling process (Fopt^min ^ Tm^th)- 
This provides a means to determine the temperature of 
the mechanical oscillator via the spectral weight of the 
two sidebands, a technique that has been widely used 
in the trappe d ion community and is termed sideband 
thermometry (Diedrich et al. 1989). There, the spec- 



tral weightcajibed^ via optical shelv- 

ing! Haroche and Raimond, 2006). A variant of this 



method has recently been de monstrated in an opt ome- 
chanical cooling experiment (Safavi-Naeini 2012). In 



this method the excitation laser is placed on the upper 
(A = -^Vtrn ) and subsequently on the lower sideband 
(A = —^rn) and the rate of scattering into the cav- 
ity TV^'^^ is measured. This yields 7V^'^^(A = Vtrn) = 
^A+(n/ + 1) and 7V^«^(A = -Qrn) = ^^-(n/) pro- 
vided the laser is quantum limited. A treatment in the 
presence of laser noise, which can yield a contribution to 
the asymmetry, is given in Appendix |XIII. C] 

Experiments and practical limitations 

The first experimental attempts to damp the motion 
of a mechanical oscillator with radiation pressure dy- 
namical backaction were carried out by Braginsky and 
coworkers (Braginsky et al. 1970). In these experiments 
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Figure 29 Experimental results for optomechanical laser- 
cooling. We display the initial and final phonon numbers 
(log scale), vs. the sideband resolution parameter Qm/i^- The 
quantum limit for the minimum achievable phonon number 
is plotted as well. MIT: ( |Corbitt "al] |2007a|) , L KB: ( [aT 
cizet et al. 2006a ), Yale: ^Thompson et al.\ 2008 ), Vienna : 



dGroblacher et al.\ l2009b| , MPQ: dSchliesser eToI] |2009|) 
JILA : ( [Teufel et aq|2008t, Cor nell 2010: ( [Rocheleau 
20T0t, Cal tech 2011: ^han et a /., 201 1), EPFL 2011: (|Ri^ 



iere et aL, 2011), Boulder 2 011: (|Teufel et a/.||2011at , MIT 
2011: ( |Schleier-Smith et al] [MTT] ) 



a microwave cavity was employed, and a modification 
of the damping rate of the end-mirror pendulum could 
be observed. Microwave cooling deeply in the resolved- 
sideband regime was moreover implemented in the field 
of gravitational wave detectors, in the form of a high Q 
cryogenic sapphire transducer, where it se rved the role 
to reduce the effective noise temperature (Cuthbertson 
[7fa/.H1996p . 

Optical feedback cooling of a micromechanical mirror 
using the radiation pressur e force was demonstrated in 
1999 (Cohadon et all 1999). Dynamical backaction cool- 



ing in the optical domain was achieved using photother - 
mal forces in 2004 (Hohberger Metzger and Karrai[ 2004), 



and us ing radiation pressu re forces in 2006 
2006a[ [Gigan et al, 2006 



Schliesser et al. 



Arcizet et al. 



2006). 



These early experiments operated in the Doppler 
regime {n > ftm)- As outlined in the previous sec- 
tion, an important requirement for many applications is 
the resolved sideband regime. This was demonstrated 



in the microwave regime for the first time (Cuthbertson 



et al. 



et al. 



1996) and later in the optical domain (Schliesser 



2008b). Since then, a plethora of no vel systems 
{e^dThompson et a/.[|2008| ), ( [Teufel et a/T|[2008} , (Lin| 
et al. 2009P and many more) have been realized in the 



resolved sideband regime (see Fig. 29). All these exper- 
iments were, however, essentially room temperature ex- 
periments. 

To initialize a mechanical oscillator in the ground state 
at thermal equilibrium, the condition kBT/h^rn 1 has 
to be realized. In general, the probability to find the sys- 
tem in the ground state is related to the average occupa- 



tion number n by Pq = 1/(1 + Note that in the field 
of trapped ion experiments, the motional ground state 
can nowadays routinely be prepared with more than 97% 
probability (Diedrich et al 



1989 



Leibfried et al. 



2003). 



TTis 



Reaching the ground state (often referred to as n < 
challenging for low frequency oscillators, as the thermal 
freezeout for a 1 MHz oscillator would equate to 50 jiiK. 
These temperature are far below those attained with a 
dilution refridgerator. Ground state cooling with conven- 
tional cryogenics can therefore be reached only for GHz 
oscillators. Indeed, the first demonstration of quantum 
control at the single phonon level was demonstrated with 
a 6 GHz piezoelectric mechanical oscillator cooled to be- 



low 50 mK (O^Connell et al. 2010). 



It is challenging to reach the quantum ground state 
of micromechanical oscillators at lower frequencies. A 
widely pursued strategy has been to combine cryogenic 
precooling with dynamical backaction laser cooling. The 
precooling thereby allows to reduce the starting mode 
temperature. This technique has over the past years al- 
lowed a substantial reduction of the phonon occupancy. 

While initial cryogenic experiments had demonstrated 
cooling to a level of a few dozen quanta in the optical 
domain fGroblache r et al. 2009b Park and Wang', '2009[ 
Schliesser et al. , 200^, further experimental progress al- 



lowed to reduce the motional energy to a level close to 
the zero point motion in several experiments both in the 
microwave domain (|Rocheleau et al. 2010[ Teufel et al.\ 



2011a) and in the optical domain (Chan et aTy2011[Riv 



iere et al. 2011 ), with ground state probabilities ranging 



from 0.2 to 0.7. 

Specifically, in the microwave domain, increasing the 
coupling strength by improved cavity designs resulted 
in cooling to around 4 quanta (|Rocheleau et al. 2010) 
(Po = 0.2) and to n/ = 0.38 (Pq 0.72) using a super- 
conduc ting resonator coup led to micromechanical drum 
mode ITeufel et al, 2011a), respectively (n/ = 0.38). In 



both experiments the limitation for the occupancy was 
set by the fact that the cavity had a low, yet finite resid- 
ual thermal occupation. This also led to the observation 
of squashing in the mechanical noise spectrum data of 
this experiment for some drive powers. All these experi- 
ments have been performed in dilution refrigerators with 
base temperatures as low as about 25 mK. 

Also in the optical domain, new geometries allowed 
to improve the cooling performance to hf = 1.7 quanta 
(Pq = 0.37) for im proved microtoroidal resonators ( ,Ver-] 
|hagen et al] |2012D and to 0.8 quanta (Pq = 0.54) for a 



nanomechanical mode of a photonic crystal beam (Chan 
et a/.[[2011} . In ( [Riviere et al. |2011J a spoke-supported 



microresonator was precooled to ^ 200 quanta in a He-3 
buffer gas cryostat (ca. 750 mK temperature) and opti- 
cally pumped by a Ti:Sa laser. In this experiment, the 
occupancy was limited by the cavity decay rate and 
in addition by reheating of the sample due to the laser, 
which in the case of glass reduces the Q-factor due to two- 
l evel systems in th e amorphous glass. The experiment by 
(Chan et al. 2011 ) involved a 3.67 GHz mechanical mode 
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and an external cavity diode laser drive. The cooling was 
limited by reheating of the sample to its initial temper- 
ature of ca. 20 K inside a He-4 flowthrough cryostat and 
the fact that the mechanical Q factor exhibited a strong 
temperature dependence, as well as unwanted effects at 
higher drive powers. In all cases, higher fidelity ground 
state preparation - as achieved in ion trapping - will re- 
quire further reduction of the effects of laser heating of 
the sample, further increases of the mechanical quality 
factor or higher optomechanical coupling rates. 

The technique of radiation pressure dynamical backac- 
tion cooling has many similarities to atomic laser cooling 
dPiedrich e^"aH |1989[ [Leibfried et a/.||2003| ). Neverthe- 
less there are substantial differences. For instance, op- 
tomechanical cooling does not proceed by the use of a 
two-level system, but by a cavity that is excited with a 
coherent laser field. From an experimental point of view, 
optomechanical cooling of a mechanical oscillator to the 
ground state is impeded by several challenges. First, in 
contrast to atomic laser cooling where the atoms are well 
isolated, mechanical systems are usually coupled to a 
high temperature bath with correspondingly large mo- 
tional heating rates. In addition, optomechanical laser 
cooling is extraordinarily sensitive to laser phase (i.e. fre- 
quency) noise. In the following two sections the funda- 
mental limits of cooling are reviewed due to both laser 
phase noise and thermorefractive cavity noise. 

A selection of cooling experiments is displayed in 



Fig. 29 There, we plot the initial and final phonon num- 
bers vs. the sideband resolution parameter Q^n/i^ that 
determines the minimum achievable phonon number. 



a. Laser phase noise In the last section, the laser input 
noises were considered to be essentially quantum noises, 
i.e. the laser beam is considered a perfect coherent state. 
However, real laser systems exhibit excess noise, e.g. due 
to relaxation oscillations that can be derived from the 
dynamical equations of laser theory. In these cases, the 
laser frequency noise cannot simply be inferred by the 
laser linewidth and simplified models fail to provide an 
accurate description of the laser's phase and amplitude 
noise. Lasers exhibit generally significant excess noise 
for frequencies below the relaxation oscillation frequency, 
which can differ strongly, from the kHz — MHz range in 
NdiYAG or fiber lasers to several GHz in the case of 



diode lasers ( 


Kippenberg et al. 2011 Vahala et al. 1983 


Wieman anc 


Hollberg 


1991). This noise wih lead to 



radiation pressure fiuctuations that heat the mechanical 



oscillator, as has been discussed in ( Diosi 2008 Rabl 



et al. 2009 Schliesser et al. |2008b ), Moreover, phase 



noise impacts state transfer and entanglement generation 
( Abdi et'al] 20111 Ghobadi et al.\ 2011) and can impact 
sideband asymmetry measurements. For a more detailed 
discussion we refer the reader to Appendix |XIII.C' 



b. Cavity frequency noise In addition to laser phase noise 
- which in principle can be mitigated by properly fil- 
tered laser systems - a more fundamental limit to dy- 
namical backaction laser cooling arises from the fact 
that at finite tempererature, the cavit y frequency will 
exhibit thermodynamical fiuctuations (Gorodetsky and 



Grudinin 2004). These are due to local temperature 
fluctuations that affect the dielectric properties of the 
cavity itself (for microtoroids, microspheres, or photonic 
crystals) or of the mirror. For more details, we refer to 
Appendix pOlLBl 



B. Strong coupling regime 

We now discuss what happens when the laser power P 
is increased further. At first, this will just improve cool- 
ing, since Fopt oc P. However, as we will see in the next 
section, qualitatively new features start to appear when 
Fopt ~ or equivalent ly when g ^ n. This regime is re- 
ferred to as the strong coupling regime , where the driven 
optical mode and the mechanical mode hybridize to form 
two new modes, with a splitting set by 2g. Furthermore, 
even for lower laser drive powers interesting features in 
the transmission spectrum of the cavity appear if it is 
probed weakly in the presence of a strong drive. This will 
be the phenomenon of optomechanically-induced trans- 



parency, discussed in Sec. VII. B. 2 



1. Optomechanical normal-mode splitting 

The strong coupling regime is discussed most easily if 
we assume the non-dissipative part of the Hamiltonian 
to dominate all decay channels, i.e. g ^ hc^ F^. In that 
case, we can just consider the following part of the lin- 
earized Hamiltonian, 



H = -h/\5d)5a + hQ.J)^h - hg{5d) + 5a){h + P) . (94) 

In the most interesting red-detuned regime, where 
A ^ — ^m, we can even start our discussion employ- 
ing the rotating-wave-approximation for the coupling, 
— hg{5d)h^5ao^)^ which is the beam-splitter Hamiltonian 
of Eq. ( 30 ) . The Hamiltonian of these two coupled oscil- 



lators is then easily diagonalized by going over to the two 
eigenmodes. These modes now represent excitations that 
are hybrids between the mechanical oscillations (h) and 
the fiuctuations of the driven cavity mode {5a) around 
the strong coherent amplitude. Their eigenfrequencies 
are: 




(95) 



In particular, right at resonance A = — ^Im, one observes 
an avoided crossing, with a splitting of ujj^ — uj- = 2g 
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Figure 30 Mechanical frequency spectrum (frequency on ver- 
tical axis) as a function of laser detuning, for a strongly cou- 
pled optomechanical system. An avoided crossing, with a 
splitting of size 2g, appears when the negative detuning equals 
the mechanical resonance frequency. This is due to the hy- 
bridization of the mechanical mode (frequency Qm) and the 
driven cavity mode (effective frequency —A). 



between the two excitation branches. At this point, the 
eigenmodes are symmetric and antisymmetric superposi- 
tions of light and mechanics, with new annihilation oper- 
ators {Sa ± \f2. Far from resonance, one recovers the 
two bare frequencies —A and Vt^^ and the excitations 
become again of purely optical and mechanical nature, 
respectively. 

This is the picture appropriate for g ^ where we 
assume that n ^ Fm is the dominant decay channel. 
In the opposite case g <^ hz (assumed in the previous 
sections), the two peaks at uj± merge and the avoided 
crossing cannot be observed. 

In principle, the complete scenario, including the de- 
cay channels and the transition into the strong-coupling 
regime, is described fully by the solution of the linearized 
equations of motion, Eqs. (41) and (43 



As discussed 

in Sec. |V.B[ one can solve these coupled equations ana- 
lytically. In this way, for example, one arrives at an ex- 
act expression for the mechanical susceptibility, Eq. ( 56 ) . 



When plotting this vs. frequency for increasing values of 
^, one observes peaks at uj± that can be clearly resolved 
for g ^ n. The evolution of the mechanical spectrum 
in the strong coupling regime ^ > /t:/4 as a function of 
laser detuning A is displayed in Fig. [30| The same kind 
of analysis also applies to the transmission spectrum of 
the cavity that can be expressed via the same solution. 

However, instead of referring to these rather lengthy 
exact expressions, we can simplify things by considering 
the regime A ^ — ^m, which allows us to perform the 



rotating- wave approximation already employed above. 
We can then write down the linearized equations of mo- 
tion for the mean values: 





(96) 



Solving for the complex eigenvalues of this non-Hermitian 
matrix, we recover the expression for the two branches 
uj± given in Eq. (95) above, except with the replacements 
A A + m/2 and Vt^ ^ Vt^ - With 6 

—A — r^m, we have 



5 .A^ + Fn 




5^i{T^-H)/2 



In particular, at resonance A 



-^m = 0) we find 

^± = ^^m — 9^ ~ (f where we assumed k, Fm 

to slightly simplify this formula. Thus, the eigenfrequen- 
cies change character at the threshold g = a>:/4, where the 
root changes from purely imaginary {g < n/A) to real- 
valued {g > n/A). This corresponds to the transition 
into the strong-coupling regime, with two well-resolved 
peaks. Each of those peaks is of width (FWHM) K./2. 
This is because both of these excitations are half optical, 
and so each of them shares half of the optical decay rate. 
Away from the degeneracy point A = — r^m, the relative 
contributions of the two decay channels F^ and n get re- 
weighted according to the unequal distribution of optical 
and mechanical excitation in the two branches. 

It turns out that cooling becomes less efficient when 
one approaches the strong coupling regime. In fact, the 
exact solution of the linearized equations of motion can 
be employed to derive the appropriate modification to 
the formulas for the final occupany in the case of strong 
optomechanical coupling: 



Uf 



=4») 



2nn 



- opt 



(97) 



Here fi^p is the standard result (90) derived above for 
the weak-coupling regime g <^ and to simplify the 
expression we assumed the resolved-sideband regime n <C 
where nmin = (^/4^m)^ and Fopt = 4^^//^, as well as 
strong cooling Fopt ^ Tm, neglecting terms of still higher 
orders in n/Vt^. 

The peak splitting in the strong coupling regime {g ^ 
Fm) and the resulting modification to cooling (no- 
tably the limitation ari sing from the finite cavi ty decay 
rate) was predicted in ( Marquardt et a/. l|2QQ7|). It was 
analyzed extensively in QDobrindt et oL 



20M Wilson- 



|Rae al\ 2008); see also (Marqu ardt et a/.| |2008D an d 



the supplementary material of ( ,Groblacher et al\ 2009a ) . 
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Figure 31 Experimental observation of strong coupling of 
a mechanical oscillator to a light field ("parametric normal 
mode splitting"). From top to bottom, the light intensity and 
thereby g — ^oVncav is increased. Da ta from the experiment 
reported in (Groblacher et al\ 2009a). . 



More generally, the full expression ( Dobrindt et al. 2008 



Wilson- Rae et al. 2008) for n/ can be given, including 



any possible thermal occupation of the cavity field (which 
has been observed in microwave setups): 



Uf nth- 



nth- 



'^cav ^ 2 
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(98) 



The first experimental 
tomechanical coup l ing {g 



(Groblacher et al. 



observation of strong op- 
> ^5 Tm) was reported in 
2009a|). Subsequently, experiments 



on other setups have been able to achieve significantly 



larger ratios g/tz, see (Teufel et al. 2011bj). In addition. 



the regime of ^ > F^^t h (see below) w as reached as 
well dTeufel et a/^|2011a|b[[Verhagen a /.[ |2012[ ), where 
coherent quantum state transfer between light and me- 
chanics could take place (see Sec. X.E). 
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Figure 32 Experimental observation of the coherent coupling 
regime, i.e. g > {hv^fim^m}, between a mechanical oscilla- 
tor and a microwave cavity mode as well as an optical cavity 
mode. The data shows the splitting in the mechanical dis- 
placement spectrum as a function of increasing drive power 
(left panel, for the ca se of a superconducting drum resonator 
( Teufel et al.\ 2011b), courtesy J. Teufel) and as a function 
of the laser detuning around the lower sideband (right p anel, 
for the case of a sp oke supported microtoroid resonator (Ver- 
hagen et al. 2012| , courtesy T. J. Kippenberg). 



The photon statistics of an optomechanical system in 
the strong coupling regime have been found to display 
interesting antibunching behaviour and photon correla- 
tions (Huang and Agarwall 2010), and the same analysis 



discussed four- wave mixing in a setup driven by a strong 
pump field and another ("Stokes") field. One may also 
consider having a nonlinear medium inside the optical 
cavity. In (Huang and Agarwal 2009c), the effects of an 



optical parametric amplifier cavity on the phenomenon of 
optomechanical normal mode splitting have been studied. 

The spectroscopic signatures of strong coupling indi- 
cate whether the coupling exceeds the cavity decay rate. 
If the coupling rate exceeds also the thermal decoherence 
rate (7 = F^n^) the interaction between the mechanical 
oscillator and the light field becomes quantum coherent, 
i.e. the timescale of the mutual coupling is faster than 
the timescale for one quantum of noise to enter from the 
environment: 



g > {TnTlm,!^} 



This parameter regime of coherent coupling is a precondi- 
tion for many quantum protocols such as quantum state 
transfer between the cavity field and a mechanical mode 
(see section X.E). Experimentally, this regime has been 
reached using a supercondu c ting m embrane coupled to 
an LC circuit ( Teufel et al.\ 2011a), and in the optical 
domai n using a toroida l spoke-supported microresonator 
( iVerhagen et a/.[|2012D (cf. Figur^). 
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Figure 33 (a) Level scheme for optomechanically induced 
transparency. A strong control laser on the red sideband 
drives the optomechanical system, while a weak probe laser 
scans across the cavity resonance and takes the (b) spectrum 
of the driven system, which displays a sharp transparency 
feature (for a side-coupled toroid setup). The linewidth of 
that feature is given by the mechanical linewidth (at weak 
drivin g), much narrower than the cavit y linewidth. Adapted 
from ( [Schliesser and Kippenberg 2010). 



2. Optomechanically induced transparency 



Electromagnetically induced transparency (Fleis- 



chhauer et al, 2005) is a phenomenon which occurs in 



multi-level atoms and manifests itself as a cancellation 
of absorption in the presence of an auxiliary laser field. 
It arises from electronic interference or, in an equivalent 
picture, is due to a dark-state resonance of the excited 
state. This phenomenon has been demonstrated for cold 
atomic ensembles, giving rise to a host of phenomena, 
ranging from optical pulse storage to slowing or advanc- 
ing of light pulses. An analogous phenomenon also oc- 
curs in optomechanical systems, as predicted theoreti- 
cally in ( Agarwal and Huang 2010 Schliesser 2009) and 



analy zed further for optical pulse storage ( Chang et a/._ 
[2011). Optomechanically induced transparency has been 
observed in experiment s(Safavi-Naeini et al. 2011 Weis 
\et al\ [2010 ). Injecting a strong control laser beam into 
the lower (red-detuned) sideband of an optomechanical 
system, the optomechanical interaction causes the cavity 
resonance, as seen by a second, weak probe laser field to 
be rendered transparent. The simultaneous presence of 
a strong control laser (se"*'^^^) and a weak probe laser 
{5s-e~^^p'^) lead to a transmission \tp\^ of the weak probe 
laser given by: 



-^0<^^Xmech(^)Xaa(^) 



(99) 



Here A 



denotes the detuning of the strong 



control field from the cavity resonance cjcav, and de- 
notes the detuning between the control laser and probe 
laser, i.e. = ujp — uJc Moreover, the coupling efficiency 
^ = ^ex/^ has been introduced, and the mechanical sus- 



ceptibility^'^ Xmech(^) = -^(^ - ^m) + as wcll as 

the optical susceptibility Xaa (^) = — ^(^1 + A) + /^/2. 
Note that, with regard to the transmission, we used the 
terminology appropriate for a waveguide-coupled unidi- 
rectional cavity (e.g. whispering gallery mode resonator), 
and we follow the discussion to be found in the supple- 
Plotting the 



2010) 



mentary material of (Weis et al 
expression reveals that when the resonance con- 

dition Q ^ Qrn is met, a transparency window arises. 
When the coupling laser is placed on the lower sideband 
(A = — r^m), the expression for the transmission of the 
probe in the vicinity of 1] = (which corresponds to 
(jjp ~ cjcav) reduces to: 



4nc 



v9o 



where for simplicity we assumed critical coupling, i.e. rj = 
1/2. Evaluating the above expression for Q — Vim yields: 



C 



C 



Here C is the optomechanical cooperativity: C — 
4^o^cav/(^r^). Thus a cooperativity of unity is required 
to change the transmission to 50%. It is noted that the 
assumption underlying these theoretical considerations, 
of two coherent drive fields, requires from an experimen- 
tal point of view two laser fields that are coherent for 
times longer than the effective mechanical damping time. 

The physical origin of the transparency window can be 
understood by realizing that the beat of the probe field 
and the coupling laser induces a time varying radiation 
pressure force. If the beat frequency matches the me- 
chanical oscillation frequency, then the mechanical oscil- 
lator is driven resonantly. The driven oscillator in turn 
creates sidebands on the intracavity field. Considering 
the strong coupling laser only, in the resolved sideband 
limit, the lower motional sideband is far off cavity res- 
onance and can be neglected. In contrast, the upper 
sideband of the coupling laser, created by the mechani- 
cal motion, has precisely the same frequency as the probe 
field and is moreover phase-coherent with the probe field. 
This leads to an interference that yields a cancellation 
of the intracavity field on resonance, giving rise to the 
transparency window. The phenomenon thereby results 
from the destructive interference between reflection am- 
plitudes for photons scattered from the coupling laser and 
photons of the probe field. 

The width of the (Lorentzian) transparency feature in 
the weak coupling regime (where Fg// <C n) is given by 
the total effective mechanical damping rate: 



Here it is convenient to adopt a definition different from x^x of 
Sec.lII B 31 
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Figure 34 Transmission of the probe beam in the presence 
of a strong control field on the lower sideband (left panel, 
optomechanically induced transparency, OMIT) and on the 
upper sideband (right panel, electromagnetically induced ab- 
sorption, EIA). 
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We note that at stronger drive this feature smoothly 
evolv es into t he normal mode splitting discussed above 
(Sec. VII.B.l). The narrow transparency window is con- 



comitant with a rapid variation of the transmission phase 
of the probe beam. This implies that a pulse with a 
bandwidth smaller than Fqmit will experience a group 
delay (without distortion of the pulse) . The change of the 
phase of the transmitted light is given by the expression 
(for A = -l^rn): 



arctan 



Evaluating this expression for zero detuning {ft = fl^) 
leads to a group delay of 



1 



OMIT 



However, the bandwidth of the pulse needs to be smaller 
than the transparency window for undistorted pulse 
propagation, which limits the delay-bandwidth product 
to TgToMiT ^ 2. Using an array of optomechanical sys- 
tems, the delay-bandwidth product can be increased and 
may therefore serve as a way to store an optical wave- 
form in long lived phonons. Such cascaded optomechan- 
ical systems c an be realized using e.g. photonic crystals 
(pTang et a/.][2QTT]). 



3. Optonnechanically induced annplification and 
electromagnetically induced absorption 

If instead we consider the control laser being injected 
on the upper sideband (A = ftrn)^ an additional opti- 



cal signal in the probe beam is amplified (Massel et al. 



2011 Safavi-Naeini et al\ 



2Q11| ). The analogous effect 



in atomic physics i s referred to as electromagnetically in- 
duced absorption (Lezama et a/.( |1999| (EIA). The am- 
plification process can parametrically amplify a small sig- 
nal (provided the resonance condition is met), by virtue 
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Figure 35 Optomechanically induced transparency observed 
in the experiment. Data from a microtoroid setu p (|Weis 
et al\ [2010) (left panel) and in photonic crystals ( Safavi-_ 
Naeini et al\ |201l| (right panel) . Note that in the case of 
a photonic crystal the OMIT signature is visible both in the 
reflected and transmitted signal. 



of the constructive interference of the light scattered 
from the pump (control) to the signal (probe) frequency, 
in direct analogy to the above phenomenon of optome- 
chanically induced transparency. Theoretically the phe- 
nomenon can be described by the same equations as the 
effect of transparency, except for the fact that now the 
mechanical damping is reducing with increasing power 
on the upper sideband. The maximum gain is set by the 
maximum power which can be injected onto the upper 
sideband (A = -^Vtrn ), which is limited by the onset of 
the parametric oscillatory instability, in which the coher- 
ent amplification of mechanical motion from the noise 
occurs. The maximum average gain in this case is g iven 
by by Ga^(A = 0) 4(41]^//^)HM assel et a/.'j|2011[ ). As 
for any non-degenerate parametric amplifier, the ampli- 
fication process has t o add half a quantu m of noise for 
fundamental reasons (Clerk et al. 2010a), and the total 
added noise is given by Uadd = + 1/2 kBT/h^m in 
the presence of thermal fluctuations. In the ideal case the 
relative phase between pump and signal is not important 
for the EIA process. However, any relative phase fluctu- 
ations between pump and signal need to take place on a 
timescale long compared to the inverse effective mechan- 
ical damping rate. 



However, this optical amplification process does not 
lead to a stimulated optical amplification process, as in 
the case of optical Brillouin scattering. The reason is 
that unlike the optical Brillouin scattering case, the op- 
tical dissipation is larger than the mechanical dissipation 
{k ^ )• This implies that the mechanical mode can 
experience exponential buildup, while in a Brillouin laser 
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A. Parametric instability and attractor diagram 




Amplitude A 



Figure 36 Optomechanical parametric instability towards 
"self- induced oscillations" ("mechanical lasing"): displace- 
ment X vs. time t. In a system with a sufficiently strong 
blue-detuned laser drive, the mechanical oscillations can dis- 
play anti-damping (Feff < 0). This leads first to exponential 
growth of any initial fluctuations, which then finally saturates 
due to nonlinear effects, resulting in self-sustained mechanical 
oscillations at a stable amplitude A. 



the opposite is the case 
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In this section, we present the classical theory of the 
optomechanical instability produced by radiation pres- 
sure backaction, following (Ludwig et a/., 2008} 



Mar- 



quardt et al. 2006). Our main goal is to discuss the 



VIII. CLASSICAL NONLINEAR DYNAMICS 



amplitude A of the steady-state mechanical oscillations, 
as a function of system parameters such as laser power, 
detuning, and mechanical damping rate. We will find 
that for a fixed set of parameters, A can in general take 
on multiple stable values, corresponding to several sta- 
ble attractors of this dynamical system. That effect is 
known as dynamical multistability, and in experiments 
it may lead to hysteretic behaviour. Our discussion will 
be directly applicable for the radiation-pressure induced 
parametric instability. However, note that very similar 
physics and analogous formulas apply for the case when 
this instability is induced by photothermal forces. In fact, 
if the thermal decay time rth is long (r^ <C ^^m), one 
can re-use most parts of the discussion below, effectiv ely 
only replacing k by (see ( Marquardt et al. 2006 ) for 
a discussion of the differences)^ 

The threshold of the instability can be obtained in a 
linear analysis by demanding that 



Up to now we have mostly discussed effects that can be 
fully understood within the linearized equations of mo- 
tion, i.e. within the quadratic approximation (28) to the 



optomechanical Hamiltonian. However, the approxima- 
tion itself can be used to predict its breakdown: In the 
blue-detuned regime (A > 0), Fopt becomes negative, de- 
creasing the overall damping rate. At first, this leads to 
heating (instead of cooling), enhancing the oscillator's 
effective temperature. Once the overall damping rate 
Tm+ropt becomes negative, an instability ensues. In that 
case, any tiny initial (e.g. thermal) fiuctuation will at 
first grow exponentially in time. Later, nonlinear effects 
will saturate the growth of the mechanical oscillation am- 
plitude (Fig.|36|. A steady-state regime is reached, with 
oscillations proceeding at a fixed amplitude A. These are 
called self-induced (back-action induced) optomechanical 
oscillations. In fact, they are analogous to lasing action, 
but now in a mechanical system and with the incoming 
laser radiation providing the pump. Therefore, this op- 
tomechanical effect can also be understood as one vari- 
ant of "mechanical lasing" or "phonon lasing". Just as 
for laser s/masers, the optomechanical system displays a 
threshold and linewidth narrowing above threshold, de- 
scribed by a Schawlow-Townes type limit. 



In the EIA scenario a stimulated optical amplification is expected 
(and subsequent optomechanical Brillouin lasing) when the me- 
chanical oscillator is more strongly damped than the optical one, 
i.e when the hierarchy Fm ^ of dissipation is satisfied. This 
ensures exponential amplification of the optical field above the 
lasing threshold and prevents the mechanical mode from building 
up significantly. 



- opt 



0. 



(100) 



This will define the limits of an interval where Fm+Fopt < 
0, i.e. where the system is unstable. This interval widens 
as the laser power is increased. 

A simple argument can now be used to obtain the am- 
plitude A. We start from the ansatz 



x{t) 



Aco^iVt^t) 



(101) 



for the self-induced oscillations. This is good for typ- 
ical experimental parameters, where Fni,Fopt 
such that both damping and optomechanical effects only 
show up after many oscillation periods. We will discuss 
the breakdown of this ansatz in the following section on 
chaotic dynamics. 

From Eq. ( [TqT] ), one can obtain the time-dependence 
of the radiation pressure force F(t) (which will depend 
on A and x). In steady-state, the time-averaged power 
input due to this force, (Fx), must equal the power lost 
due to friction, meffFm ( x^). This can be recast into a 
condition resembling Eq. ( |100[ ) , by defining an amplitude- 
dependent effective optomechanical damping rate: 



- opt 
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(102) 



In the low- amplitude limit A ^ 0, this reduces to the 
standard definition of Fopt used up to now. Then the 
power balance condition is simply 



ropt(A) = 0, 



(103) 
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which is an implicit equation for A. This strategy can 
be used for arbitrary optomechanical systems, also con- 
taining more optical modes or other types of radiation 
forces. 

We need yet another condition, to fix the oscillation 
offset X, which is not identical with the unperturbed oscil- 
lator equilibrium position. The time-averaged radiation 
pressure force deflects the harmonic oscillator: 



(104) 



In general, Eqs. ( |1Q3[ ) and (104) need to be solved si- 
multaneously for the unknowns A and x. However, if 
Fm <C l^m, one can already see the instability in a regime 
where the shift x is small and can be neglected, such that 
only Eq. (103) is relevant. 



Mathematically, the onset of small-amplitude oscilla- 
tions, starting from A = 0, is an example of a Hopf 
bifurcation. In this regime, A ex — /th, where / is 
any system parameter (such as the detuning or the laser 
power), and /th its threshold value. 

An important feature is the dynamical multistability, 
i.e. the existence of several stable solutions for a fixed set 
of external parameters. This is observed for sufficiently 
good mechanical quality factors, when higher-amplitude 
attractors become stable. It leads to hysteresis in ex- 
periments, and might also be used for high-sensitivity 
"latching" measurements. For more on this and the ef- 
fects of noise and slow dynamics of the amplitude, see 
dMarquardt et a/.||2QQ6|. 



In experiments, an optomechanical instability due to 
retarded light forces was first studied in a low-finesse 



We still have to obtain F{t). This can be deduced by 
solvi ng the classical equation for the light field amplitude 
(Sec.[llFCt, 



setup with photothermal forces (Hohberger and Karrai 



2004), where the retardation is due to finite thermal 



conductivity and the theory described here applies with 



appropriate modifications, see (jMarquardt et al. 2006 



Metzger et al. 2008c). Subsequent studies of the pho- 



a = - 2 ~ Q^max) + ^(A + Gx{t))a , (105) 
where we defined amax to be the amplitude reached in- 



side the cavity right at resonance (in terms of Sec. |III.C 
we have am a 



= 2ai ^^JK,Q^/ hv). After i nserting the ansatz 

(101), the solution (Marquardt et al. 2006) can be writ 
ten in a Fourier series a{J) 



tothermal setup observed parts of the attractor diagram, 
confirmed dynamical bistability, and uncovered a new 
regime where more than one mechanical mode gets in- 
volved in the nonlinear dynamics ( Metzger et al. , 2008b ) . 

The parametric instability driven by radiation pressure 
backaction, as discussed here, was first demonstrated in 



coefficients 



^ ^ inCL^t -4-^ a mic r otoro i d setup (Carmon et al. 2005 Kippenberg 



et al. 2005 Rokhsari et al. 2005). The fuh attractor 



an 



Jn{- 



GA \ 



i{Gx^A)/K. ' 



(106) 



where Jn is the Bessel function of the first kind and the 
global phase (f{t) = {GA/Q^) sm{Q^t). Now the force 
F{t) = hG\a{t)f and the time averages (Fx) and (F) 
can be calculated. One has (\a{t)f\ = |Q^n|^ and 



diagram still has to be observed in an experiment. 

A recent experiment has demonstrated mechanical las- 
ing (i.e. coherent oscillations) in a setup where two op- 
tical modes are involved and photon transitions between 
those modes provide the power to feed the mechanical 
oscillations (Grudinin et al.\ |2010 ). The attractor dia- 



gram for the parametric instability in systems involving 
more than one optical mode (including a 'membrane in 
the middle') will display qualitatively new features due to 



a{t)\ xj = Al^mlm These series can be 

efficiently summed numerically, to obtain the explicit de- 
pendence of Eqs. (103) and (104) on A, x and the system 
parameters. 

The result is best discussed graphically. In Fig. [37| 
we show the attractor diagram, i.e. the possible ampli- 
tudes A as a function of any system parameter (in this 
case, the detuning A). This diagram can be generated 
by plotting the value of ropt(^, A) and then showing the 
contour lines Fopt = — Fm that indicate possible attrac- 
tors. Note that stable attractors are only those where 
I Fopt I grows with A (the upper half of each line). In 
the case shown here {n/Vt^n = 0.2), the structure of side- 
bands at A = nl^m both for red and blue detuning shows 
clearly. Remarkably, one can have stable self-induced os- 
cillations even on the red-detuned side (A < 0), but only 
for finite amplitude A > 0. This is consistent with the 
fact that the linearized theory there predicts cooling at 
A^O. 



the effects of optical Landau-Zener dynamics (Wu et al. 
2011). 

Just as in a laser, the phase of the self-induced oscilla- 
tions is arbitrary. Thus, external noise, including thermal 
Langevin forces acting on the mechanical oscillator and 
radiation pressure shot noise, will impart a slow phase 
diffusion. The effect of a force SF on the phase scales 
inversely with the amplitude, S(f oc SF/A, as can be 
seen easily in a phase-space diagram. Thus, the diffu- 
sion constant for the phase scales with 1/A^, diverging 
just above threshold (Vahala^^2008^ ). This is the optome- 
chanical analogue of the Schawlow-Townes result for the 
linewidth of a laser (Schawlow and Townes 1958). More 



precisely, the scenario is closer to the case of a maser, 
since the thermal noise is not negligible in the mechani- 
cal oscillator: 



— (2nth + l). 



(107) 



A full discussion of the linewidth narrowing and phase 
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Figure 37 Attractor diagram for "mechanical lasing" in an 
optomechanical system. We show the damping rate Fopt as 
a function of both mechanical oscillation amplitude A and 
laser detuning A, for a system in the resolved sideband regime 
(k/O^ui = 0.2). Contour lines of Fopt correspond to possible 
attractors. For each line, only the upper part is stable, i.e. 
those points on the curve where Fopt increases when going 
towards higher amplitudes. The offset x has been assumed 
to be negligible. (Plot shows Fopt^eff^m/(2/iG'^amax) vs. 
GA/Qm and A/Qm; contours shown at levels —10~^ /white, 
-10"Vorange, -10"^/black ) 



diffusion in optomechanical os cillations can be found in 
([Rodriguez and Armour 2010), where the effects of pho- 



( [Rod 

ton shot noise are taken into account as well. 

Up to now, we have discussed exclusi vely the classi- 
cal dynamics. In the quantum regime (Ludwig et al 



[2008), the parametric instability threshold is broadened 
due to quantum fluctuations, with strong amplification 
of fluctuations below threshold. The existence of attrac- 
tors with finite amplitude shows up, e.g., in phase-space 
plots of the Wigner density, and it changes the phonon 



(as well as photon) statistics (Qian et al. 2012) 



Interesting new collective physics results if many such 
optomechanical oscillators are coupled together. In par- 
ticular, this can lead to synchronization, which we will 
discuss in some more detail in the section on multimode 



optomechanics (Sec. IX). 



B. Chaotic dynamics 

If the laser input power is increased sufficiently, the 
coupled motion of the light field and the mechanical os- 
cillator becomes chaotic. In that regime, amplitude and 
phase fluctuate in a seemingly random fashion that de- 
pends sensitively on initial conditions, even in the ab- 
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Figure 38 (a) Experimental observation of dynamical mul- 
tistability for optomechanical oscillations. The mechanical 
oscillation amplitude vs. static displacement (or equivalent ly, 
detuning) displays hysteresis upon sweeping in different direc- 
tions. A is the optical wavelength for this setup. Data from 
( Metzger et al.\ 2008c'), where the radiation force was of pho- 
tothermal origin, (b) Experimental results for the linewidth 
narrowing above threshold, as a function of oscillation ampli- 
tude, obtained for a microtoroidal setup in ( ,Rokhsari et al.\ 
20061. 




Power [a.u.] 

Figure 39 The transition towards chaotic motion in an op- 
tomechanical system has been observed in (Carmon et al. 
,2007J (courtesy Tal Carmon). 



sence of noise sources. Technically, the ansatz (101) of 



sinusoidal oscillations breaks down, and the full dissipa- 
tive driven dynamics of four degrees of freedom (x, p and 
the complex light amplitude a) has to be taken into ac- 
count. The chaotic regime is characterized by positive 
Lyapunov exponents, where any tiny deviation from the 
initial trajectory grows exponentially with time. 

Chaotic motion in optomechanical systems has been 
explored relatively little so far, a l thoug h it had been 
observed early on (Carmon et al. [2005 ) and has been 
studied more systema tically in a subsequent experiment 
([Carmon et a/.l|2007[). 



IX. MULTIMODE OPTOMECHANICS 

Up to now we have almost exclusively considered one 
optical mode coupled to one mechanical mode. This is 
the "minimal model" of cavity optomechanics, captured 
by the Hamiltonian (22). Of course, it is clear in prin- 



ciple that every mechanical resonator has a multitude of 
normal modes, and every optical resonator likewise has 
many different modes as well. It is relatively straightfor- 
ward to write down the appropriate extension of ( 22 ) to 
the more general case: 
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k j 

-hY,[9o]iiaiai{b,^b])^... (108) 

Here the various optical {dk) and mechanical (bj) modes 
interact according to the optomechanical coupling con- 
stant tensor [goYj^i = [^o];^^, whose entries depend on the 
details of the optical and vibrational modes and their 
mutual interactions. We left out the laser drive and the 
coupling to the radiation and mechanical environments. 

Before we go to the general case, it should be noted 
that restricting one's attention to the minimal model 
is often justified for many purposes. The incoming 
monochromatic laser drive will select one optical reso- 
nance. With regard to the mechanical motion, all the 
mechanical resonances will show up in the RF spectrum 
obtained from a displacement measurement, but we may 
choose to focus on one of those resonances, as long as 
they are well separated. Likewise, cooling or heating in 
the resolved sideband regime singles out a particular me- 
chanical mode via the choice of laser detuning. In the 
bad cavity limit, mult iple modes can be cooled (or ampli- 



fied) s imultaneously (iBagheri et al. 2011 Metzger et al. 
2008cl). 



In the following, we want to discuss some scenarios and 
features where it becomes crucial to go beyond the min- 
imal model. It is clear that going to structures with two 
or more mechanical or optical modes leads to a wealth of 
different possible schemes (Fig. 40), only a few of which 
have been explored so far. Further examples can be found 
in the next section, on quantum optomechanics (Sec.|X|. 

It was pointed out by Braginsky that scattering of pho- 
tons betwee n two optical modes can lead to a parametric 
instability (Braginsky et al. 2001). This analysis was 



intended primarily for interferometric gravitational wave 
observatories, where the free spectral range may match 
relevant mechanical frequencies. 

We have already mentioned a setup with a membrane 
in the middle of a cavity (Sec 



VI.B.2I), in the context of 



QND phonon detection (Jayich et al] 2008 Miao et al. 



2009 Thompson et al. 2008). This setup can be viewed 



as consisting of two optical modes, where photons can 
tunnel between those modes via transmission through the 
membrane. Having two optical modes is essential for the 
quadratic dependence of optical frequency on displace- 
ment that arises from the avoided crossing in the optical 
mode spectrum. 

Besides the QND scheme, setups with two relevant op- 
tical modes can be realized in a large variety of implemen- 
tations and have many additional interesting features. 
Mechanical oscillations can take the system through the 
avoided crossing, potentially resulting in Landau- Zener- 
Stueckelberg physics and optical Rabi oscillations for 
the optical two-state system, coherently shuffling pho- 
tons between the two branches ( (Heinrich et al.^ |2010| ). 



optical 
mode 



(c) 

optical waveguide 



mechanical 
mode 





acoustic waveguide 



Figure 40 Schematic illustrating the possible variety of mul- 
timode structures for optomechanical circuits, (a) Two me- 
chanical modes coupled to a common optical mode (for en- 
tanglement etc.) (b) Two optical modes coupled to a mechan- 
ical mode (QND phonon detection etc.) (c) ID array, with 
mechanical couplings between the cells, (d) 2D array, with 
coupling via common optical modes, (d) Optical and acous- 
tic waveguides feeding into localized optical and mechanical 
modes (e.g. photon-phonon translator). 



Cooling and squeezing in the presence of quadratic op- 
tomechanical coupling has been analyzed (Nunnenkamp 



et al. |2010). Moreover, displacement sensing in a two- 
mode se tup can be much more efficient tha n the standard 
scheme (Dobrindt and Kippenberg 2010). The optome- 



chanical system of two coupled optical modes can also be 
viewed as a photonic version of the Josephson effect and 
its classical dy namics can give rise to chaos ( |Larson and 
Horsdall[20TT|. 



If the mechanical resonance frequency matches the 
transition between the optical branches, one can imple- 
ment a version of optomechanical phonon lasing ()Gru- 



dinin et al. 



2010) that is directly based on a popula- 



tion inversion between the two optical levels, just like 
in a real laser. The Landau- Zener physics mentioned 
above will significantly modify the dynamics of these self- 
oscillations for large amplitudes (Wu et al. 2011). An 



interesting variant of this situation can be implemented 
for whispering gallery optical and acoustical modes in 
toroids or spheres in a stimulated Brillouin scattering 
scheme (Bahl et al. 2011 Tomes and Carmon[ [2009 ) . In 
the same kind of Brillouin setup, the reverse process has 
been demonstrated as well, i.e. cooling via a photon t ran- 
s ition between the lower and upper optical mode (Bahl| 



et al. 2012 Tomes et al. 2011). Recently it has been 



pointed out theoretically that schemes with two optical 
modes can be exploited to enhance the quantum nonlin- 
earities in optomechanical systems fLudwig et al. 2012 
Stannigel et al. 2012), see Sec. X.F The photon-phonon 



translator discussed in Sec. |X.E| is another example of a 
possible device that employs two optical modes. 

Likewise, there are many schemes where more than 



45 



only a single mechanical mode is relevant. These could 
be the various normal modes of a given mechanical struc- 
ture, or several vibrating objects placed inside a cavity or 
coupling to one optical mode (see Sec. X.F). We mention 
only a few illustrative examples. An array of multiple 
membranes inside a cavity was studied theoretically in 



(iBhattacharya and Meystre 2008 Hartmann and Plenio 
[2OO8; Xuereb et al. , 2012) with regard to mechanical nor- 
mal modes, entanglement and collective interactions. A 
recent experiment demonstrated tripartite optomechan- 
ical mixing between one microwave mode and two me- 



chanical modes (Massel et al. 2012). Effects of multiple 



mechanical modes can also be studied for levitated se- 
tups. It has been proposed that an array of levitated 
dumbbell-shape dielectric objects can undergo an order- 
ing transition. 

The nonlinear dynamics of the self-oscillations ("me- 
chanical lasing"; Sec. VIII) become particularly interest- 
ing when more mechanical modes are involved. It has 



been shown experimentally (Metzger et al. 2008b) that 



collective self-oscillations may arise when several mechan- 
ical modes are excited simultaneously, using a strong 
blue-detuned drive for a "bad cavity" whose linewidth 
encompasses several mechanical normal modes of the 
structure (a cantilever in that case). More recently, 
it was pointed out that optomechanical systems are 
very promising for obs erving synchronization phenomena 
(Heinrich et al. 2011). Assume an array of optomechan- 



ical oscillators, each of which consists of a mechanical 
mode coupled to an optical mode that is driven by a 
blue-detuned laser beam such as to go into the "mechan- 
ical lasing" regime (Sec. VIII). If these "clocks" are cou- 



pled mechanically or optically, their mechanical frequen- 
cies can lock to each other, even if they have been dis- 
tinct initially. Under appropriate conditions, a variant of 
the Kura moto-model can be d erived for optomechanical 
systems (Heinrich et al. 2011), which is a paradigmatic 



model of synchronization physics. For two mechanically 
coupled cells, the equation turns out to be of the form 



them implemented as a double-disk SiN structure. Me- 
chanical spectra showed the onset first of self-oscillations 
and then of synchronization as a function of laser detun- 
ing. 

Two coupled mechanical modes of widely different 
damping rates can give rise to Fano lineshapes in their 
excitation spectrum when they hybridize. This has been 



demonstrated in an optomechanical system (Lin et 



2010), and suggested for information storage and re- 



trieval in long-lived mechanical "dark" states. 

Optomechanical photonic crystal structures present an 
opportunity to design more complex optomechanical cir- 
cuits. These might be ID or 2D array structures, where 
many optical and mechanical modes are arranged in a 
periodic layout and coupled to each other. Alternatively, 
one can think of more intricate circuits, more similar to 
computer chips, where different elements fulfill various 
functionalities (sensing, amplification, general signal pro- 
cessing). 

Optomechanical arrays of this type have been studied 
theoretically only in a few works so far, both with respect 
to their collective classical nonlinear dynamics (Heinrich 
et al.\ 2011), their quantu m many-body dynamics ( Lud- 
wig and Marquardt[ 2012 ) a nd for engineering quantum 
dissipation (Tomadin et al. 2012), as well as with re- 



spect to quantum applications: A suitably engineered ar- 
ray of optical and mechanical modes coupled to a waveg- 
uide can slow down and store light |Chang et al. 2011). 
Photon-phonon entanglement in an optomechanical ar- 



ray of three cells was studied in (Akram and Milburn 



2012). Many nanomechanical modes in an array geome- 



try can be entangled via the light field, using a suitable 
parametric drive to select mode pairs (iSchmidt et c 



2012). These studies pave the way towards future archi- 



tectures for (continuous variable) quantum information 
processing with optomechanical circuits. 



Sip = 5Q — K sm{25ip) , 

where 6Q is the intrinsic frequency difference, 6(p the dif- 
ference of mechanical oscillation phases, and K an effec- 
tive coupling constant that can be related to microscopic 
parameters. If K is large enough, synchronization ensues 
{dip = 0). If, one the other hand, several oscillators cou- 
ple to the same optical mode, the behaviour can become 
of a form that is not described by any Kuramoto-type 
mode l . This has been analyzed recently in detail (Holmes] 
et a/.l 120121). 



Synchronization may be important for metrological 
applications, where several synchronized optomechani- 
cal "clocks" of th is type are expecte d to be more sta- 
ble against noise ( [Tallur et al. |2010 ). Experimentally, 
some signs of synchronization in an optomechanical de- 
vice have been observed recently ( Zhang et a/.[|2012[ ) for 
two optically coupled optomechanical oscillators, each of 



X. QUANTUM OPTOMECHANICS 

Quantum mechanics has already figured at several 
places in our discussion so far, notably in setting the lim- 
its for displacement sensing or cooling. We will now turn 
to discuss potential future optomechanical experiments 
where quantum behaviour will take center stage. We will 
discuss ways to create interesting quantum states both in 
the optical and mechanical system, and to create entan- 
glement between the various subsystems. We will then 
turn to nonlinear quantum effects whose understanding 
requires us to go beyond the linearized optomechanical 
interaction, i.e. beyond the quadratic Hamiltonian of 



Eq. (28). Finally, we will see how these effects are envis- 



aged to form the ingredients of future optomechanically 
aided protocols for quantum information processing. 
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A. Light-assisted coherent manipulation of mechanics 

The light field can be employed in principle to generate 
arbitrary quantum states of the mechanical oscillator. In 
this section, we will restrict ourselves to the action of the 



linearized interaction Hamiltonian of Eq. ( 28 ) , 



{lin) 



(109) 



When injecting Gaussian optical states, as is the case for 
the usual laser drive, this can produce arbitrary mechani- 
cal Gaussian states, i.e. coherent and squeezed states out 
of the ground state (which has to be reached first, either 
via optomechanical or bulk cooling techniques). As long 
as the linearized Hamiltonian is valid, non-Gaussian me- 
chanical states can only be produced from non-Gaussian 
optical states. The parameters that can be varied easily 
are the laser detuning A and the laser input power. Any 
pulse of light will generate a radiation pressure force pulse 
that shifts the oscillator's wave function and thereby per- 
mits to create a coherent state. 

It is only slightly more difficult to create a squeezed 
state. As we saw in Sec. V.B.1[ a far-detuned light 
beam creates an "optical spring effect", i.e. a change 
in the mechanical frequency by SQ^n = 2ncav5'o/^- 
well known, a time-dependent modulation of l^m, i-e. a 
parametric driving of the mechanical oscillator creates 
a squeezed state: see (Mari and Eisert 2009) for an 



analysis in the case of optomechanical systems, where 
one can employ a modulation of the laser power. For 
SQjji{t) = SQcos{2Q^t)^ the resulting effective mechani- 
cal Hamiltonian turns out to be (in a frame rotating at 
^Irn, and in rotating wave approximation) the standard 
squeezing Hamiltonian 



^mech = ^(6^ + (6t)^), 



(110) 



such that b{t) = cosh{6Qt)b{0) - i smh{6Qt)b^ {0) in the 
absence of dissipation. Thus, in this time-dependent 
scheme, squeezing grows exponentially with time. How- 
ever, a realistic analysis needs to take into account the 
initial thermal population, as well as dissipation and de- 
coherence, and distinguish between the steady-state sit- 
uation and the transient case. When several mechanical 
modes are coupled to the same optical mode, a mod- 
ulated laser drive will generate two-mode squeezing or 
beam-splitter interactions between pairs of modes that 
are selected according to their frequency. This can form 
the basis for continuous varia ble quantum state pro cess- 
ing in optomechanical arrays (Schmidt et al. 2012). 



Measurements can also be used to generate inter- 
esting mechanical states in a probabilistic manner, 
i.e. conditioned on the measurement result. This in- 
clude s squee zed states via single-quadrature detection 
(Sec. VI.B.1[ ) or mechanical Fock states via phonon num- 
ber readout (Sec. VLB. 2). Further below we will com- 



(including non-Gaussian states), by nonlinear effects 
(Sec. |X.F" ) or various state transfer protocols (Sec. X.E). 



B. Mechanics-assisted readout and manipulation of light 

The optomechanical interaction can be exploited to de- 
tect and manipulate the quantum state of the light field. 
An example that was suggested early on is the pos- 



sibility of a QND detection (jBraginsky et al. 1980) for 


the light intensity circulating 


inside the cavity (Bragin- 


sky et al\ 1977 


Jacobs et al.\ 


1994 


Pinard et al\ 1995). 



The displacement of the end-mirror, induced by the ra- 
diation pressure force, can serve as a noiseless meter for 
the light intensity. In addition, as the radiation pres- 
sure force is proportional to the photon number, it in- 
creases in discrete steps. If the photons are sufficiently 
long-lived and the interaction strong enough, this may 
even enable QND photon detection, by registering the 
resulting mechanical displacement. In practice, however, 
that regime is extremely challenging to reach, as it would 
require ^o/(^^m) ^ 1 (see the section on quantum non- 



linear dynamics below. Sec. X.F). 



Regarding the manipulation of the light's quantum 
state, one of the most straightforward applications of op- 
tome chainc^^onsist sjn_s 

beam (Fabre et al\ 119941 iMancini and Tombesil 11994). 



In this context, the change of the cavity length due to 
the intensity-dependent radiation force can be compared 
to the effect of a Kerr medium inside a rigid cavity. The 
resulting physical picture depends sensitively on the de- 
tuning and the frequency at which the noise of the light 
beam is analyzed. In the simplest case, one may imagine 
that a temporary ffuctuation in the incoming intensity of 
the light beam induces a change in the cavity length via 
the radiation pressure force. This, in turn, shifts the op- 
tical resonance and thereby affects the circulating and 
outgoing intensities, potentially suppressing the noise. 
It should be noted that for a single-sided optical cavity 
(without internal losses) there can be no change of the 
amplitude noise at zero frequency, since the incoming and 
outgoing intensities have to be equal. However, there can 
be amplitude squeezing at finite frequencies (and phase 
squeezing down to zero frequency). The squeezing ef- 
fect diminishes towards frequencies above the mechan- 
ical resonance. Finite temperatures degrade squeezing, 
as the thermal motion of the mirror imprints extra noise 
on the light beam. At low frequencies, and for detun- 
ings on the order of z^:, the light beam's noise is increased 
(Fabre et al. |1994 ) by a factor nth/Qm, where nth is the 
equilibrium phonon number. 

A recent review of the current efforts towards demon- 
strating radiation pressure shot noise effects and sq ueez- 
ing in optomech anical experiments can be found in ( |Ver 
lot et al. 2011). First experiments in this direction have 



ment on other ways to generate more nonclassical states 



simulated the q uantum fiuctuations by classical inten- 
sity noise ( Marino et a/.[|2010| ). More recently, the very 
strong coupling and low temperatures attainable in re- 
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alizations of cavity optomechanics with atomic clouds 
have allowed to obtain first signa tures of genuine opti- 
cal squeezing (Brooks et al. 2011) at the quantum level. 



C. Optomechanical entanglement 

The optomechanical interaction can be used to engi- 
neer entanglement between the light beam and the me- 
chanical motion, or between several light beams or sev- 
eral mechanical modes. However, a prerequisite for all of 
these approaches is to cool the mechanical oscillator to 
near its ground state and to achieve sufficiently strong 
coupling. 

In the following, we ffist describe entanglement in the 
single-photon strong coupling regime (for more on this 
regime, see Sec. X.F) and then turn to continuous vari- 



able entanglement. 

A particularly simple physical picture applies in an ide- 
alized situation where we imagine starting with a super- 
position of photon states inside the cavity. Then, the 
radiation pressure force assumes different, discrete val- 
ues for each photon number, displacing the mechanical 
harmonic oscillator potential by 2xzpF^o/^m per pho- 
ton. If initially the mechanics was in its ground state, 
it will evolve into a coherent state |<^n(^))vib' oscillating 
around the new, displaced origin that depends on the 
photon number. Thus, we immediately arrive at a non- 
factorizable, i.e. entangled state of the form: 



^)cav® l"n(*)>vib 



(111) 



n=0 



Here n is the photon number, are the arbitrary initial 
amplitudes for the photon field (e.g. corresponding to a 
coherent state), and (pn{t) a phase-shift that can be ob- 
tained by solving the time-dependent Schrodinger equa- 
tion for this problem. The state (111) can be interpreted 



as a "Schrodinger cat" type state, where a "microscopic" 
degree of freedom (the optical cavity mode) is entangled 
with a "macroscopic" (or mesoscopic) degree of freedom, 
the vibrating mirror. This picture has been ffist analyzed 
in dBose et a/.[ |1997| [Mancini et a/.||1997| ). 

Several signatures of entanglement exist. In this exam- 
ple, where the overall state is pure, we can simply trace 
out the mechanical vibrations, arriving at the reduced 
density matrix of the optical field, which will be found 
in a mixed state whenever there is entanglement. How- 
ever, at multiples of the mechanical period, light and me- 
chanics completely disentangle, since the coherent state 
l^n(^))vib ^^^^ h.dive returned back to the origin, indepen- 
dent of photon number. At these times, the photon state 
becomes pure again, even independent of the mirror's ef- 
fective temperature. The optical field's coherence thus 
demonstrates decay and revivals. This can in principle 
be tested in an interferometric optomechanical "which- 
way" experiment, where a photon can take either of two 
paths, one of which contains an optomechanical cavity 




Figure 41 An optomechanical which-path experiment can be 
employed to test for the decoherence of mechanical superpo- 
sition states ( Marsh all et al.\ |2003| ) . A photon would leave 
behind a vibration in the mirror of an optomechanical cavity, 
which destroys the photon's coherence with the other path. 
However, if the photon exits after a full mirror oscillation 
period, it will be fully coherent again, unless mechanical de- 
coherence has occured in the meantime. 



(see Fig. [41[). The revivals can be observed in the inter- 
ference visibility, as a function of the time the photon 
has spent inside the cavity. However, any mechanical de- 
coherence occuring in the meantime will spoil these per- 
fect revivals of the photon field's coherence. This can in 
principle be employed for highly sensitive optical tests of 
sources of mechanical dissipation and decoherence and 



fundamental quantum physics in general (Bose et 



shall et al. 



19991 IBouwmeester et a/.H 1998; Fol man et a/.[|2QQl||Mar^ 



2003). Experiments of this kind could quan- 



tify the decoherence of superpositions of heavy objects 
(i.e. the mirror), and thus potentially shed new l ight 
on the quantum-to-classical transition (see Sec. XI. A). In 



particular, mode ls of gravitationally induced decoherenc e 
might be tested fBose a/.l|1999[ |Marshan et a/.[|2003[ ), 
where the (admittedly tiny) hypothetical extra decoher- 
ence rate becomes potentially observable only for rela- 
tively massive objects. Such experiments will require the 
challenging regime of go / Q^n > 1 (when the displacement 
induced by a single photon may become on the order of 
the mechanical zero-point motion). It has been recently 
suggested that nested inte rferometry allows to reduce this 
requirement significantly ( [Pepper et al 2012). 

On the other hand, entanglement can also exist on 



the level of continuous variables ( Braunstein and van 
Loock, 2005), of the type ffist proposed in the famous 



Einstein- Podolsky- Rosen article (Einstein et al. 1935). 
The analysis of entanglement in this context has to be 
performed by taking into account the dissipative nature 
of the systems involved, and typically relies on solving 
the linearized quantum equations of motion (including, 
if needed, input-output theory for treating the reflected 
or transmitted light beam). In the regime where lin- 
earized equations of motion are valid, Gaussian states 
of the mechanics and the light field will be produced. 
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and their entanglement (of the continuous- variable type) 
can be characterized completely once the correlations be- 
tween the various mechanical and optical quadratures are 
known. A typical measure of entanglement that is com- 
monly applied here is the logarithmic negativity, which 



can be calcula ted easily for Gaussian states (Vidal and 
Werner[ |2QQ2), both pure and mixed. 



The entanglement between mechanical vibrations and 
the optical cavity field (already described above in a pic- 
ture appropriate for strong coupling) has been analyzed 
in more deta i l for the continuous- varia ble case in (Pater- 
nostro et al. |2QQ7 Vitali et al 2007j). Recently, it was 
pointed out that a suitable time-dependent modulation 
of the drive can improve the efficiency of photon-phonon 
entanglement ( Mari and Eisert[ |2Q12p . 

It is also possible to create entanglement of two spa- 
tially separate mirrors. When a strong pump beam 
runs through a nonlinear optical x^^^ medium acting as 
a nondegenerate optical parametric amplifier, two- mode 
squeezing produces entanglement between the quadra- 
tures of pairs of emanating light beams. This entangle- 
ment could then be transferred via the radiation pressure 
force onto two spatially separated mirrors that are part 
of optomechanical cavities on which these light beams 
impinge (Zhang et al. 2003). In this way, optomechan- 



ics would help to create mechanical EPR-type entangle- 
ment at a distance. The verification of EPR entangle- 
ment between macroscopic test masses by sensitive mea- 
surements has been studied in more detail for the con- 



text of_gravita^^ setups (Miao 



etal 2010 Miiller-Ebhardt a/. 2009). In the context 



of optomechanics, we usually consider the system to be 
driven by a coherent laser beam. However, it is natural 
to ask whether special opportunities arise when the light 
that is injected displays quantum features. For example, 
injecting squeezed light can be beneficial for entangling 
nanom echanical resonators via the o ptomechanical inter- 
action (Huang and Agarwal, 2009b). 

In another approach, one could achieve the same goal 
without the optical entanglement created by a nonlin- 
ear medium, and instead perform an optical measure- 
ment after the interaction has taken place, in an en- 
tanglement swapping scheme. If two independent light 
beams interact with separate optomechanical cavities, 
then the beams can afterwards be brought to interfere 
at a beamsplitter, and a suitable Bell state measure- 
ment can then be used to generate entanglement between 
the d istant mechanical re sonators ( Pirandola et al 2006 ) 
(see (iBorkje et al 2011) for a similar proposal). 

Alternatively, the driven optical field inside the cav- 
ity automatically induces an effective interaction between 
several mechanical modes, thus providing yet another 
way to generate mechanical entanglement, without the 
need for any optical nonlinearities or entanglement swap- 
ping schemes. In the case of two mechanical resonators 
(or two normal modes of one resonator) coupling to the 
same driven cavity mode, this may be understood as a 
consequence of the optical spring effect (Sec. V.B.I ). For 



the case of a single mechanical mode 6, eliminating the 
driven cavity mode by second-order perturbation theory 
creates an effective interaction term fi(^Qncav/A)(^+S^)^- 
Proceeding through the same argument for the case of 
two mechanical modes gives rise to an effective mechan- 
ical interaction. It is of the form 



rrint 
^eff 



^#^[(^l+^l) + (^2 + 



(112) 



if we assume for simplicity that both mechanical oscilla- 
tors couple equally strongly to the optical mode. 

In order to successfully entangle different mechanical 
modes, one has to laser-cool those modes, since the me- 
chanical vibrations would be far from their ground state 
for typical bulk temperatures. The optically induced 
steady-state entanglement between two movable mirrors 
under si multaneous laser-cool ing has been studied care- 
fully in ( Pinard et al] 2005), and later for so mewhat 
different setups in ([Hartmann and Plenio 2008 |Muller-] 
Ebhardt et al] 2008 ) . There is an interesting caveat for 



such studies: It is insufficient to apply the Markov ap- 
proximation to describe the dissipative dynamics of the 
mechanical vibrations in this context, even though for 
many other purposes in optomechanics that is a very re- 
liable approach. In fact, it can be shown that in a proper 
treatment (Ludwig et al. 2010) there is an optimum in- 



termediate laser-cooling strength for which entanglement 
is maximized (an effect entirely missed by the common 
Markovian treatments). 



The light-induced interaction given in Eq. (112) can 



form the basis of a general scheme for quantum state 
processing with many nanomechanical modes. In order 
to selectively address pairs of such modes for entangle- 
ment and state transfer, one simply has to modulate the 
coupling strength (i.e. the laser intensity) at sum and 
difference frequencies of those modes. Such a parametric 
scheme only requires one appropriately modulated laser 
input to address whole arrays of mo des, if the proper 
layout is chosen (Schmidt et al. 2012). 

Another possibility consists in exploiting optomechan- 
ics to entangle two light-beams. In these cases, the op- 
tomechanical interaction essentially serves the purpose of 
a x^^^ nonlinear medium. This has been proposed e.g. for 
an optomechanical setup where two degenerate, orthog- 
onally polarized cavity modes are driven strongly and 
their interaction with the moveable mirror creates EPR 



correlations (Mancini et al. 2001) between the quadra- 



ture variables of the beams emanating from these modes. 
Another option is to entangle the two sidebands reffected 
from a vibrating mirror, which works even in the absence 



of a cavity, for a strong short incoming laser pulse (Pi- 



randola et al. 



2003). For this situation, a more detailed 
analysis of the entanglement between mirror vibrations 
and the full light field (infinitely many degrees of free- 
dom) was performed in ( |Miao et al. 2010). Optically 
trapped mirrors in a cavity optomechanics setup can als 



be exploited to entangle light beams (Wipf et al. 2008). 
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D. Quantum hybrid systems 

Optomechanical systems already represent a quantum 
hybrid system, i.e. a coupling between two quantum sys- 
tems of a different physical nature: light and mechanical 
vibrations. In general, hybrid approaches may be useful 
for purposes such as quantum information processing, in 
order to combine the advantages of different physical sys- 
tems in one architecture. Some systems may be strongly 
interacting (good for computation), some are very co- 
herent (good for long-term storage), and yet others are 
easily transported over long distances (good for commu- 
nication). 

In principle, other components may be added easily to 
optomechanical setups. This is because both the light 
field and the mechanical vibrations are quite versatile 
in coupling to a variety of systems, such as cold atoms, 
spins, superconducting and other electronic qubits, etc. 
Consequently, there are already several proposals along 
these lines. 

In a cloud of atoms, the total spin can sometimes be 
treated as a harmonic oscillator, identifying its small fluc- 
tuations around a preferred direction as position and mo- 
mentum quadratures. This picture is useful when dis- 
cussing experiments where the state of the light field 
i s transferred to the atom ic spin state and back again 
(iHammerer et al. 2010a). More recently, it has been 
suggested that light might also be used to couple the 
collective spin of an atom cloud to a nanomechanical os- 
cillator (Hammerer et al. 2009a^. In such a setup, a 



light beam passing first through an optomechanical cav- 
ity and then through an atom cloud carries information 
in its quadratures about the sum and the difference of 
position and momentum variables of the mechanics and 
the spin state. A subsequent QND measurement then is 
able to prepare the two systems in an EPR state, condi- 
tioned on the measurement result. As another example, 
coupling to the internal transitions of atoms can enhance 
cooling (Genes et al. |2009b) . One can also consider more 



involved internal level schemes for the atoms. For exam- 
ple, a vibrating mirror could be coupled strongly to the 
collective spin of a cloud of three-level atoms display- 
ing electromagnetically- induced transparency phenom- 
ena (Genes et al. 2011). In addition, it has been sug- 



gested that coupling the collective spin of an atom cloud 
parametrically to a resonator can lead to the phenomena 



known from optomechanics (Brahms and Stamper-Kurn 
2010| ) (amplification and cooling of the spin, frequency 



shifts, and squeezing of light), with the collective spin 
replacing the mirror motion. 

After the pioneering experiments on optomechanics us- 
ing cold atoms (Sec. IV. G), a number of different pos- 
sibilities have been explored theoretically to couple the 
motion of atoms to other systems. For example, the light 
field inside the cavity may be used to couple the motion 
of a single atom to the vibrations of an end-mirror or 



a me mbrane ( [Hammerer et al.\ |2009b| [Wallquist et al. 



within reach. This could be the basis for exploiting all 
the well-known tools for manipulating and reading out 
the motion of a single trapped atom to gain access to the 
membrane motion. It would also be an interesting system 
to observe the entanglement between a microscopic and 
a macroscopic degree of freedom. Another model system 
that has been studied ( [Hammerer et a/1 2010b) and is 
now being implemented ( [Camerer et al.\ 2011) is a cloud 
of atoms in a standing light wave reflected from a vi- 
brating mirror (without a cavity). This could allow long- 
distance coupling between atoms and mechanical objects 
spaced apart by macroscopic distances. 

A variety of other ideas exist for merging cold atom 
systems with optomechanics. Both the optomechanics of 



Bose- Einstein condensates (jChen et al. 2010[ Jing et al. 



2011 Steinke et al. 



Fermi gases (Kanamoto and Meystre 



2011[) and degenerate cold atom 
20101) have been 



explored theoretically in some depth. Even the Mott- 
insulator to superfluid transition of atoms in an optical 
lattice coupled to a vibrating mirror has been analyzed, 
as an example of a strongly interacting quantum system 



subject to the optomechanical interaction (Chen et al. 



2009 



Larson et al. 



2008) 



The idea of doing optomechanics on trapped atoms has 
found an interesting counterpart in proposals for doing 
optomechanics on levitated dielectric objects (Sec. [IV.E ). 
The promise of this approach lies in a drastically en- 
hanced mechanical quality factor. Some proposals con- 



sidered trapping and cooling a mirror (iBhattacharya and 
Meystre[[20Q7a|b[[Bhattacharya et a/.[l2008nSingh et al. 



2010), which has the advantage (over other objects) that 



scattering of the light into unwanted directions is greatly 
reduced. Alternatively, one can have dielectric spheres or 
other p articles trapped in an optical lattice or by other 
means (Barker', [20101 [Barker and Shneider[ |2010| 'ChangJ 



et 



2010; Romero-Isart et al. 2010). A more detailed 



analysis of fu ndamental applications and protocols can 
be found in (Pender et al.\ |2012| [Romero-Isart | [2011 



Romero-Isart et al. 2011a|b| . If the promise of very long 



mechanical coherence times is fulfilled, then these plat- 
forms could offer the best means to test novel decoher- 



ence mechanisms (see the discussion in Sec. X.C). Some 
early experiments and studies on c avity optomechanics 
with sub- wavelength nan o-objects ( [Favero and Karrai 
Favero et al. 2009| ) (and Sec. IV. E) have already 



2008 



explored the optical coupling and some of the scattering 
mechanisms that may become relevant in this domain, 
even though they did not yet benefit from a suppression 
of mechanical dissipation. 

Connecting the world of superconducting or other 
solid-state qubits to optomechanical systems represents 
an intriguing possibility in the context of quantum infor- 
mation processing. This has become particularly relevant 



since the p ioneering experiment at UCSB (O'Connell 
et al. 20To| that demonstrated strong coupling between 



2010), where the strong coupling regime seems to be 



a superconducting phase qubit and the GHz oscillations 
of a piezoelectric nanoresonator, swapping a single exci- 
tation from the qubit into the resonator. First experi- 
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merits have demonstrated how to manipulate a mechan- 
ical nanoresonator both via the optomechanical inter- 



action and electrically (Lee et al. 12010 Winger et al. 



[2011), which is an important ingredient for a future hy- 
brid platform of solid-state qubits with electrical inter- 
actions coupled to mechanical nanoresonators coupled 
to the light field. Theoretical proposals have already 
pointed out how to use systems of this type to map solid- 
state quantum information into photons and back again 



(Regal and Lehnert 


2011 Safavi-Naeini and Painter 


201 1| Stannigel et al 


I |2011| |2010| |Tian| |2012t |Wang| 


and Clerk 2012). More elementary, the nanomechanical 



structure can serve as an intermediary to generate entan- 



glement between microwave and optical fields (Barzanjeh 
[e/a/.n2011D . 

An equally promising avenue is to merge the fields 
of solid-state quantum optics and cavity optomechanics. 
First experiments deliberately introducing semiconduc- 
tor materials in cavity optomechanical setups now exist. 
These include 2/im-diameter GaAs vibrating disk struc- 
tures with very high op tomechanical coupling strength 
dPing et al\ |2Q11 [ |2010D and coupli ng of the li ght field 
to a semiconductor nanomembrane (lUsami et a/.l 12012). 



The excitonic transitions of quantum dots embedded in 
such materials could couple to the mechanical vibrations 
either directly via deformation potentials or indir ectly 
via the light field (Rundquist and Majumdar 2011). 



E. Quantum protocols 

The previous sections introduced some basic quantum- 
physical features in optomechanical systems, such as pro- 
ducing and reading out nonclassical states of light and 
mechanics. As soon as this can be achieved reliably, one 
may envisage building a toolbox for quantum manipu- 
lation in these systems, and exploit it for purposes of 
quantum communication and quantum information pro- 
cessing. This would follow the pioneering ideas and ef- 
forts in the ion trap community, where it has been sug- 
gested early on that one can exploit the motional degrees 
of freedom to facilitat e quantum gates between the in- 
ternal states of ions ( [Cirac and Zoller 1995 Leibfried 
et al. 2003 ) . Micro- and nanomechanics offers the added 



value that they can be functionalized and hence couple 
to many different physical degrees of freedom. Optome- 
chanical devices therefore offer a fruitful addition to the 
vast array of physical systems that are being explored for 
quantum information processing ( jZoller et a/!|[2QQ5| . 

One of the most prominent protocols is quantum state 
transfer. In the context of optomechanics, this would al- 
low to realiably convert an optical pulse into a mechanical 
excitation (and vice versa). In principle, this is straight- 
forward, since the linearized optomechanical interaction 
describes a coupling between two oscillators (mechani- 
cal and driven cavity mode) , which can be tuned via 
the laser intensity (Sec. III.B). For the red-detuned case, 
A = — r^m, in the resolved-sideband limit, we found a 
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Figure 42 Optical wavelength conversion in an optomechani- 
cal setup: Two modes at different frequencies are illuminated 
by strong red-detuned control beams. When a signal is in- 
jected at the resonance of mode 1, it is converted into a signal 
emanating from mode 2. 



beam-splitter type of interaction: 



(113) 



Thus, the excitations Sa on top of the strong coher- 
ent laser drive can be swapped onto the mechanical 
resonator. Such an operation would be performed in 
a pulsed scheme. The coupling g{t) becomes time- 
dependent via the laser intensity, and it would be 
switched on for just the right amount of time to per- 
form a complete state swap between the two oscillators 
{Sa and b). Two laser pulses are needed for this scheme. 



A red-detuned control pulse at uj^ 



control 



Qrn de- 



termines the time-dependent coupling g{t). A second, 
"signal" pulse serves to excite the ^a-oscillator into some 
target state that then will be written onto the mechan- 
ics. Since Sa oscillates at the frequency —A = Q^n in the 
frame rotating at cc;£°"*^°^ , the signal pulse has to be in- 



I.e. 



jected at a frequency co^l^^^^ = cj£°"*^°^ 
right at the cavity resonance. Ideally, the whole swap- 
ping pulse sequence is shorter than the cavity decay time, 
which however requires a two-mode setup with different 
decay rates for the modes (since the control pulse cannot 
be shorter than I/hc for a single mode). A detailed anal- 
ysisof swapping protocols can be found in (Wang and] 



Clerk||2012). 

There is also a simple classical picture for what hap- 
pens during such a pulse sequence: The superposition of 
signal and control beam leads to a beat-note in the in- 
tensity at ftrn- This translates into a radiation pressure 
force that resonantly excites the mechanical oscillator. 

The optomechanical interaction can also be exploited 
for conversion between different optical wavelengths. 



That scheme is illustrated schematically in Fig. 42 It 
involves two optical resonances, each of which is driven 
by a strong control beam, red-detuned by Qrn- When an 
input signal (e.g. a pulse) is injected at the resonance 
of mode 1, it will be converted to an output signal em- 
anating from the resonance of mode 2. This has been 



demonstrated recently in (Dong et al. 2012 Hill et al. 
^2012^. 

First experimental proof-of-principle demonstrations 
have shown t hese concepts in the high-temperature, clas - 
sical regime (Fiore et al. 2011 Verhagen et al. 2012). 
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A = A = 
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Figure 43 The linearized Hamiltonian of cavity quantum op- 
tomechanics describes three different kinds of interaction be- 
tween the mechanical resonator and the driven cavity mode, 
depending on which laser detuning is chosen. This is the ba- 
sis for elementary quantum protocols, such as storing optical 
information into the mechanical degree of freedom. 



In this context it was pointed out ( Romero-Isart et al 



2011b Verhagen et al. 2012) that the requirement for 



coherent state transfer of this type is g > k,^ Tm^th- The 
next challenge will be to drastically improve th e fidelity 
and demonstrate true quantum state transfer ( Parkins 
and Kimblel |1999| [Ritter et 'al] |2012[ ), e.g. by reco n- 



structing the mechanical quantum state (Sec. VI.B.l). 



Recently, state transfer between a mechanical mode 
and an itinerant microwave coherent state has been re- 
ported (Palomaki et al. 2012). In this experiment, a 



suitably shaped microwave pulse was written onto the 
motional state of a micromechanical membrane and later 
retrieved via quickly switched control beams, hence real- 
izing a coherent mechanical memory for microwave pulses 
in the weak coupling regime. Another protocol that does 
not require the strong c oupling regime is based on quan- 
tum state tel eportation (Hofer et al. 2011 Romero-Isart 
eran|20TTb|. 



If a single photon is sent into the setup and transferred 
to the mechanical resonator, this will prepare the res- 
onator in a Fock state or some nonclassical state in gen- 
eral dAkram et a/. 1 12010[ [Kahili et a/.||2010D. That may 



be the most efficient route towards generating nonclassi- 
cal mechanical states in optomechanical systems as long 
as the single-photon strong coupling regime (Sec. X.F) 
has not been reached. 

Up to now, we have only described how single localized 
phonons (stored inside the mechanical resonator) can be 
converted to photons and back again. An equally inter- 
esting, or perhaps even more useful, scheme would take 
traveling phonons and convert them to photons. Such a 
device has been proposed recently and has been termed 
an opt omechanical "traveling wave phonon-photon trans- 
lator" ( Safavi-Naeini and Painter] [2011 ). Although the 
frequency is shifted by many orders of magnitude, the 
wave function of the outgoing single photon is designed 
to be a faithful replica of the incoming phonon's wave 
function. 

The basic idea is the following: Phonons are travel- 
ing down a phononic waveguide and enter a localized 



phononic mode, where they experience the usual optome- 
chanical interaction. Frequency upconversion is achieved 
by having a high-intensity stream of incoming photons. 
In a Raman-type scattering process, the phonon com- 
bines with one of those photons to form a single photon 
at a slightly different frequency, so energy conservation is 
obeyed. The strong pump beam and the weak stream of 
outgoing converted photons can be efficiently separated 
by being coupled to two different optical modes. Indeed, 
a suitably engineered opto mechanical structure (Safavi- 
Naeini and Painter) |2011|) has two optical modes cou- 
pling to the phonon displacement field in the following 
way (which is conceptually identical to the membrane- 
in-the-middle setup discussed in Sec. VI.B.2): 



H = -hgo{b + b'^){a\a2 + + • • • • 



(114) 



If mode ai is pumped strongly, we can replace the field 
operator by the classical amplitude: ai{t) ^ ae~^^^^ (as- 



sume a real- valued). If 002 
terms to be retained are: 



H = —hgQaiha^e 



' cJi + r^m, then the resonant 



-\-iujLt 



) + ... (115) 



This clearly displays the elementary process of convert- 
ing a single phonon to a single photon (and back). It 
also shows that the conversion rate can be tuned via the 
pump strength a. All the remaining challenge lies in en- 
suring that 100% of the phonons arriving at the device 
are indeed converted into photons. This is essentially 
an impedance-matching problem, since the coupling of 
the phonon mode to the phonon waveguide (ideally set 
by Fm if other losses can be neglected) is usually much 
weaker than the coupling of the optical mode to the pho- 
ton waveguide (set by Hi2)- Without extra fine-tuning, 
most of the phonons would be reflected. This can be 
overcome by a judicious choice of the coupling. Indeed, 
from t he p oint of view of the phonon mode, the cou- 
pling ( |115[ ) to the lifetime-broadened photon mode a2 
gives rise to a Fermi golden rule transition rate 4^o<^^//^2. 
This is the rate at which a given localized phonon would 
be converted and decay into the photonic waveguide, 
and it is identical to the optomechanical cooling rate 
in the sideband-resolved regime. By matching this to 
the coupling to the phononic waveguide, i.e. demanding 
Ag^o? / i<i2 = r^, one creates a situation that is equiv- 
alent to a two-sided cavity with equal mirrors, where 
100% transmission can be achieved on resonance. Conse- 
quently, in the present setup ideally 100% of the phonons 
can be converted if 2gQa = \/k2T^. A detailed anal- 
ysis fSaS\d^^ae^^ 2011) considers the full 



scattering matrix that describes scattering of incoming 
phonons into the photon waveguide (or reflection back 
into the phonon waveguide), and it includes the unwanted 
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effects of extra intrinsic losses and noise^^. 

On a classical level, the device described above takes 
the slow amplitude and phase modulations of the phonon 
field, i.e. of a soundwave traveling down the waveguide, 
and transposes them into the optical domain by shifting 
the carrier frequency from mechanical frequencies (e.g. 
GHz) up to optical frequencies. The fact that the band- 
width is set by the smallest damping rate in the problem, 
which is Fm, can be an advantage if one uses the device as 
a narrow-bandwidth frequency filter. For example, two 
photon-phonon translators in series (the first one oper- 
ated in reverse, i.e. going from photons to phonons) im- 
plement a potentially very narrow filter in the optical 
domain. 

Phonon-photon conversion has been analyzed for 
transferring solid-state quantum information to the opti- 



cal domain (Stannigel et al. 2011 2010 Wang and Clerk 



[2012") and for transferring optical pulses between differ- 
ent wavelengths and pulse shaping (Safavi-Naeini and^ 
Fainter! [2QTT| |Tian[ [2QT2{ |Tian and Wang[ |201Q^ ' 



One example for the more advanced possibilities of op- 
tomechanical quantum protocols is to perform continuous 
variable quantum teleportation. Besides employing this 
for generating entanglement between distant mechanical 
oscillators (as mentioned above), one may also teleport 
an arbitrary input state of the light field onto the me- 
chanics. The generic idea is to start with the entangle- 
ment between the mechanical motion and a light beam, 
and then to let the beam interfere at a beamsplitter with 
another light beam, carrying an arbitrary input state. 
A subsequent measurement in both output ports of the 
beamsplitter (Bell measurement) then yields a classical 
measurement result on the basis of which one manipu- 
lates the mechanical state, leaving it in a final quantum 
state that equals the arbitrary input state. Such a scheme 
has been analyzed for a strong short laser pulse imping- 
ing on a v ibrating mirror in a f r ee-beam setup (withou t 
a cavity) ( [Mancini et al. 2003 Firandola et al. 2003), 
where the refiected optical Stokes and anti-Stokes modes 
at cjl ± get entangled with the vibrations, an d for a 
time-dependent drive of an optomechanical cavity ( Hofer 



et al. 2011 Romero-Isart et al. 2011b). 



F. Nonlinear quantum optomechanics 

The optomechanical interaction 

—hgQo)a{h + b^) 

is cubic in field operators, i.e. the corresponding Heisen- 
berg equations of motion are nonlinear. However, in 
experiments this nonlinearity so far has only played a 



In that reference, k, and 7 refer to half the photonic and phononic 
intensity decay rates 1^ and Fm employed in our notation. 



role in the classical regime of large amplitude oscillations 
(both mechanical and with regard to the light field), see 
Sec. [VIII.A In the quantum regime, we have thus far 
resorted to the "linearized" description, with a quadratic 
interaction Hamiltonian of the type —hgo^^^n^{Sa^ + 
Sd){b + 6^), see Eq. (28). This linearized approach is 



good enough to understand many facets of cavity op- 
tomechanics: displacement detection down to the SQL 
(Sec. VI. A), t he theo ry of optomechanical ground-state 
cooling (Sec. VILA), light /mechanics hybrid ization i n 

> tz (Sec. 



VH.B.l 



VH.B.2 



the strong-coupling regime go\^n^ 
optomechanically induced transparency (Sec. 
optomechanical squeezing of light (Sec. X.B), and almost 
all of the various entanglement and state transfer schemes 
presented in the previous sections. The experimental ad- 
vantage of the linearized interaction is that its strength 
9 = 9o V ^cav can be tuned at will by the incoming laser 
power. In this way, a small value of go (fixed by the 
setup) may be compensated for by a stronger laser drive, 
until technical constraints become important. 

The disadvantage of relying on the linearized inter- 
action is that, by itself, it will always turn Gaussian 
states (of the light field and the mechanics) into Gaussian 
states. These may be squeezed or entangled, but they 
will never have a negative Wigner density, which may 
be required for certain quantum applications. It should 
be noted that there are some ways around this restric- 
tion, by introducing a nonlinearity at some other stage 
of the experiment: For example, one may send in single- 
photon pulses and then transfer these Fock states onto 
the mechanics using the linearized interaction. Another, 
sometimes equivalent, option is to perform single-photon 
detection at the end, thereby creating nonclassical states 
via post-selection of events. These strategies are there- 
fore related to the general schemes that have been ex- 
ploited already for linear optics quantum computation by 
adding sing le-photon sources and photodetectors ( |Knill| 



et a/.||200lD . 



Thus far, we have encountered only two ideas on true 
nonlinear quantum optomechanics: T he optic al QND de- 
tection of the phonon number (Sec. |VI.B.2| ) is such an 
example, and it indeed would prepare (probabilistically) 
Fock states of the mechanical oscillator. Another ex- 
ample, discussed already very early in the literature, is 
the optomechanical "Schrodinger cat" type of entangle- 
ment, where a single photon should ideally be able to 
displace the mechanic al osc illator by about a mechanical 
zero-point width (Sec. X.C). 



These examples require a large value of go, which is a 
challenge. We remind the reader that in a typical (Fabry- 
Perot type) setup the value of go can be estimated as 



90 



XZPF 



1 



L V 2meff^n 



(116) 



Here L is the effective size of the cavity , and meff the 
effective mass. Both can be made small by miniaturiz- 
ing the setup, and consequently record values of go are 
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achieved in micrometer-size devices, like photonic crys- 
tal nanobeams or very small disks and toroids. This is 
even in spite of the fact that miniaturization also drives 
up the mechanical frequency Qrn- In such setups, go cur- 
rent ly take s values on the order of up to some MHz (see 

Sec.IrO:). 



(Sec. VII.B) already for ncav = 1- Alternatively, go/hz 



Next, we should discuss in which sense go can be 
"large". The steady-state displacement produced by a 
single photon on average (ncav = 1) is 



Sx 

XZPF 



Thus, to displace by more than the zero-point width (me- 
chanical ground state width), one needs go > ^m- How- 
ever, if the photon decay rate is large, then one can 
see only the average displacement produced by the pho- 
ton number fluctuating around ncav = 1, and one would 
not resolve the granularity of the photon stream. Such 
a situation should still be well described within the lin- 
earized approximation. To obtain truly nonlinear effects, 
one would rather like to make sure the following pic- 
ture applies. Take any single photon entering the cavity. 
If its lifetime is large enough {Q^ ^ the resolved- 
sideband regime), it will displace the oscillator by the 
amount given in Eq. (117). This then implies ^ ^ as 
a necessary (but not sufficient) condition. 

These considerations directly lead us to consider the 
ratio 



90 



(118) 



If that ratio is larger than one, then the presence of a 
single phonon would shift the optical frequency by more 
than a cavity linewidth. More precisely, the mechani- 
cal displacement produced by a superposition of and 1 
phonons would be measured so efficiently that the pas- 
sage of a single photon through the cavity already de- 
stroys the superposition. We can also look at the light 
field's back-action: In the "bad cavity" limit {n ^ ^m), 
a single passing photon with a lifetime K:~^ gives an av- 
erage momentum kick 6p = 2pzpF9o/i^ to the mechani- 
cal oscillator, which would be larger than its zero-point 
momentum uncertainty pzPF = ^eff^m^zPF if 9o > ^• 
One can make a connection to the Lamb-Dicke param- 
eter used in ion-trap physics, by defining the ratio of 
momentum kick to momentum zero-point fluctuations: 

^Lambd-Dicke = SpXzPv/^ = Sp /{2pzpf) = 9o/l^- 

The ratio go/ Hi has been called the "granularity param- 
2QQ8| ), as ^o//^ > 1 allows to reveal 



eter" (Murch et al 



the granularity of the photon stream (the discreteness 
of individual photons). It has reached values on the or- 
der of and larger than one in experiments with clouds of 
ultracold atoms (see Sec. |IV.G ). 

There are other ways of interpreting ^o/^ > 1 as well: 
The (sideband-resolved) cooling rate ^o^cav/^ would 
be so large that one enters the strong-coupling regime 



can be written as the ratio between the mechanical zero- 
point fluctuations and the width of the optical cavity 
resonance, expressed in terms of a displacement (where 
the optomechanical coupling enters): 



90 



XZPF 
^^cav 



where Sxc^v = i^/G = HiXzPF/90' 

There is an additional interesting aspect about go / Hi. It 



(117) can s erve as a "quantumness" parameter (Ludwig et al. 



2008), with larger values denoting a gradual classical- 



to-quantum crossover. All the parameters of any given 
standard optomechanical setup can be boiled down to 
the following five dimensionless combinations: 



On ^ 

r ' n 

J- m 



9_ 90 

Hi ' K, 



(119) 



Here the first four do not depend on the value of Planck's 
constant. This is obvious for the first three (the sideband- 
resolution ratio, the mechanical quality factor, and the 
laser detuning in units of mechanical frequency). It is 
less obvious for g/tz. However, this can be written as 



Hi 



90 



^opt 



(120) 



Here E^cav = ^^opt^cav is the light energy stored inside 
the cavity, which is connected to the laser driving power. 
In this sense, g/ni serves as a dimensionless classical mea- 
sure of laser power. Onl y the ratio ^o/^ depends on 
as can be seen from Eq. (116): go \fh. 



Thus, one can imagine keeping all four classical ratios 
fixed, and only increasing gol^ by changing parameters 
in the setup. This is then completely equivalent to in- 
creasing Planck's constant, allowing one to resolve more 
and more quantum features as goj^ grows. It should 
be noted that of course even for gol^ ^ 1 one can ob- 
serve quantum effects, but only in the linearized regime. 
Proving the quantumness of these effects produced by lin- 
earized interactions usually requires a quantitative com- 
parison (e.g. with the light field's or the oscillator's zero- 
point fluctuations). In contrast, some features observed 
for larger values of gol^ may even be qualitatively dis- 
tinct from classical predictions. An example noted above 
are the quantum jumps of p honon n umber that could be 
observed for go > /^abs (Sec. VI.B.2). 

The third important ratio involving the coupling is 



This is the ratio between the strength g^ / ft^ of the effec- 
tive photon-photon interaction induced by the mechan- 
ics, and the optical linewidth. When this starts to be 
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larger than one (and K: <C ^^m), then the presence of 
one photon shifts the resonance sufficiently that a sec- 
ond photon cannot enter the cavity. That leads to the 
photon-blockade phenomenon (see below). 

In (Ludwig et al 2008), the regime of ^ 



was considered using both full numerical master equation 
simulations and a quantum Langevin approach. There, 
it was found that for increasing values of go/i^ quan- 
tum fluctuations start to have a pronounced effect on 
the mechanical lasing instability that is observed at blue- 
detuned laser driving. The strongly enhanced suscepti- 
bility of the system just below the threshold amplifies 
the effects of these fluctuations, and the threshold is 
smeared and shifted. In this quantum regime, the co- 
existence of several attractors (known from the classical 
case, see Sec. VIII. A) results in non-Gaussian mechan- 



ical Wigner densities and mechanical states with non- 
Poissonian phonon distributions and large Fano factors 
(Ludwig et al 2008). It was found that in this regime 



for appropriate parameters one can even generate true 
nonclassical mechanical states, with partially negative 
Wigner densities ( jQian et al 2012). These states are 
present in the steady state (under constant drive), so the 
Wigner densities could then be read out according to the 
schemes presented in Sec. 



In (Nunnenkamp et al. 



VI 



2011), master equation sim- 



ulations were further extended to discuss the full range 
of detunings and the excitation spectrum of the cavity. 
Multiple optical sidebands are found, and the mechani- 
cal state of the oscillator is seen to develop non-Gaussian 
states particularly at detunings which drive multi-photon 
transitions. 

The nonlinear quantum optomechanical regime leads 
to very interesting photon correlations. In particular, 
under appropriate conditions one may observe optome- 
chanically induced photon blockade ( Rabl , 2011) . This 
shows up in the photon-photon correlations as strong 
anti-bunching, g^'^\t = 0) < 1, which has been calculated 



and discussed in (Rabl 2011) for the case of weak laser 
driving. The regime of optomechanical photon blockade 
requires sideband resolution {n < ft^) as well as strong 
single-photon coupling with > n and > Q^k,. In a 
more recent work, the analysis of photon correlations in 
this regime was extended to cover the full temporal evo- 
lution of g'^'^^ {t) and the Fano factor, as well as higher mo- 



ments of the photon counting statistics ( Kronwald et al. 



[2013} 



In (Hong et al. 2011), an exact solution was presented 



for the regime where only a single photon is transmitted 
through such a strongly coupled optomechanical setup, 
which is important for the generation of entanglement 
(Sec, pre ). 



The understanding of the strongly nonlinear quantum 
regime is aided b y the following picture, alrea dy discussed 
partiall y in Sec. X.C| and ffist employed in (Bose et al. 



|1997; Mancini et al.\\1997 



ncav5 the mechanical equi 

amount 2xzPF^cav^o/^m 



: For any given photon number 
ibrium position is shifted by an 
One should thus consider the 



mechanical Fock states in this new displaced parabolic 
potential. 

If an additional photon enters the cavity, the potential 
suddenly shifts, but the mechanical wave function at ffist 
remains the same. Thus, the overlaps of the given initial 
wave function and the displaced new Fock states will de- 
termine the strength of possible transitions. These over- 
lap integrals are known as Franck- Condon factors from 
the theory of molecules, where vibrations may be excited 
during electronic transitions. 

We conclude by briefly describing the formal treatment 
via the well-known polaron transformation, as applied 



to an optomechanical system (Groblacher et al. 



Mancini et~al] |1997| [Nunnenkamp et al.\ |2011| |Rabl' 



2009a 



2011). We consider the standard optomechanical Hamil- 



tonian (in the frame rotating with the laser frequency) 



H = — flAa^a+fll^m^^^— %0^^^(^+^^) + ^<^L(^+^^)+^bath, 

where q^l is proportional to the laser amplitude. We 
complete the square to obtain 



- -^^^^j - h{gl/Vtm){o)af 

(121) 

for the ffist three terms in H. This shows two things: 
First, an effective photon-photon interaction is gener- 
ated, viz. the {d)a)'^ term, which is crucial for nonlinear 
effects and quantum gates. 

Second, n photons shift h by (^o/^m)^, i-e. they 
shift the equilibrium oscillator position by 2xzpf(^o/^)^ 
to the right. ^This can be accomplished by a unitary 
U = ex.-p[{b^ — b){go/^)a^a] acting on the wave functions. 
After applying that transformation to the Hamiltonian, 
via H^^"^ = WHU^ we obtain a Hamiltonian that would 
be diagonal in the absence of driving and decay: 



(122) 



Note that in these new coordinates the addition of one 
photon shifts the wave function to the left, whereas the 
center of the mechanical oscillator potential is now forced 
to remain fixed (in contrast to the actual physical situa- 
tion, where the potential shifts and the wave function is 
fixed). The exponential (polaron operator) in the laser 
driving term generates the Franck-Condon overlap fac- 
tors. 

The energy level scheme for the nonlinear quan- 



tum optomechanical regime (Fig. 44) displays equally 
spaced phonon ladders whose offset shifts depending 
on the photon number: ^^(nphonon, ^cav)/^ = —(A + 
^§/r2m^cav)^cav + ^mnphonon. In this schcmc, the driving 
laser only induces energy-conserving transitions, horizon- 
tal in the diagram. Different transitions are activated 
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photons 



Figure 44 The energy level spectrum of an optomechani- 
cal system at zero detuning A = 0: E{nphonon,ncav)/h = 
— (^o/^m)ncav + ^mnphonon- Somc posslblc proccsscs are 
indicated: Blue-detuned laser drive at the first sideband 
(A ^ +^lm) leading to creation of a phonon and subsequent 
escape of the photon out of the cavity (blue), a similar cy- 
cle at red-detuned drive (A ^ — ^^m) leading to cooling by 
annihilation of a phonon (red), photon-blockade prohibiting 
transitions ^ 1 ^ 2 towards two photons (green), and a 
resonant two-photon transition 0^2 (orange) via a virtual 
intermediate state (dashed). Note that the last two, nonlinear 
effects require go/^m > i^- 




AoM 2 3 4 



Figure 45 Single-photon strong coupling regime, (a) Illustra- 
tion of displacement depending on photon number (leading 
to Franck-Condon physics), (b) Transmission sidebands, (c) 
Photon-photon correlator as a function of detuning (for weak 
drive). [(a),(b) courtesy of A. Nunnenkamp, (c) courtesy of 
P. Rabl; g in this figure refers to our ^o] 



upon changing the detuning. Multi-photon transitions 
become possible via virtual non-resonant intermediate 
states (e.g. adding two photons at once). 



XI. FUTURE PERSPECTIVES 

The fast experimental and theoretical advances in cav- 
ity optomechanics during the recent years are constantly 
opening up new avenues with respect to applications and 
tests on the foundations of physics. Here we remark 
briefly on the broader outlook. 



A. Foundational aspects 

The ability to achieve coherent quantum control over 
the center of mass motion of massive mechanical objects 
provides a fresh approach to fundamental tests of quan- 
tum theory in a hitherto unachieved parameter regime. 
Specifically, quantum optomechanics offers a universal 
scheme for experiments in the quantum regime of mas- 
sive mechanical objects from clouds of 10^ atoms or 
nanometer-sized solid-state devices of 10^ atoms and a 
mass of 10~^^ kg, to micromechanical structures of 10^^ 
atoms and 10~^^ kg, to macroscopic centimeter- sized ob- 
jects for gravitational wave detectors comprising more 
than 10^^ atoms and weighing up to several kg. In prin- 
ciple, this offers a range of almost 20 orders of magnitude 
in mass and 6 orders of magnitude in size for macroscopic 
quantum experiments. 

A specific example where quantum optomechanics pro- 
vides a new direction for future experiments is the quan- 
tum measurement problem, which addresses the ques- 
tion why quantum superpositions do not seem tooc- 



cur at the level of macroscopic objects (Leggett 2002). 
Various new theories and phenomena beyond quantum 
theory have been suggested in order to achieve an ir- 
reversible decay of superposition states, i.e. decoher- 



ence, into well-defined classical states (Adler and Bassi 


2009 


Bassi and Ghirardi| 2003| Diosi 19841 |Karolyhazy 


1966 


Penrose 


1996 


). Each of these approaches predicts 



a particular scaling of the decoherence rates with parti- 
cle number or mass, and with the actual distinctness of 
the states involved in the superposition. For sufficiently 
macroscopic systems and sufficiently distinct superposi- 
tion states, these predictions deviate significantly from 
the decoherence rates expected from standard quantum 



theory ( Schlosshauer 2008 Zurek 1991 2007). Current 



matter- wave experiments with molecules may soon start 



2011 



iNimmrichter 



to enter such a regime ( Gerlich et al. 

et al. 20111). Systematic tests of the validity of quan- 



tum theory necessarily also involve tests of such scaling 
laws and the large mass range offered by quantum op- 
tomechanical systems will provide a unique opportunity. 
One way of producing quantum states involving super- 
positions of mechanical states is via optomechanical en- 
tanglement (see Sec. X.C). Probing the decoherence of 
such optomechanical superposition states, for example 
via the interference visibility in a single-photon inter- 
ferometer dBose et~al\ |1999| iBouwmeester et al\ |1998 



cisive tests of specific 'collapse' models (Bassi et al. 2005 



Marshall et al. 2003 [Pepper et al. 2012), may allow de- 
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Kleckner et al. |20Q8 ). A particularly exciting perspec- 
tive is to extend these experiments to the large masses 
that are available in gravitational wave interferometers 
(see (|Chen, 2013) for a recent review) and that were re- 
cently cooled to 200 thermal quant a above the quantum 
ground state (Abbott et al. 2009| ). Another route that 



has been suggested is to analyse the contrast of matter- 
wave interference of levitated nano-objects, where super- 
positions of macroscopically distinct position states ar e 
generated via optomechanics ( Romero-Isart et al 2011). 
It has been shown that such experiments would in princi- 
ple allow to enter a regime in which all non-sta ndard de- 
coherence theories can be systematically tested ( Romero- 



Isart 2011). The demanding experimental boundary 
conditions with respect to temperature and background 
pressure (to minize the effects of standard decoherence) 
might require the added benefit of a space enviroment. 
First studies along this line are currently being performed 



(Kaltenbaek et fl/.[|2012D . 

Another fascinating long-term perspective is the possi- 
bility to make use of the accessible large masses in quan- 
tum optomechanics experiments to explore the scarcely 
studied interface between quantum physics and gravity. 
Some of its aspects are already covered by the decoher- 
ence tests discussed above, as the models of Karolyhazy 



( Karolyhazy 


1966), 


rose (Penrose 


1996 



1966D, Diosi ( |DiQsi[ |1984| |2007| and Pen- 



1996 2000) identify gravity as the domi- 
nant player of their state- vector collapse. A completely 
different approach ha s been taken by a recent proposal 
(Pikovski et al. 2011| ) that suggests that quantum op- 
tomechancis experiments could test directly predictions 
from quantum gravity. Specifically, the availability of 
large masses in combination with quantum optical state 
preparation and readout is shown to be sensitive to pos- 
sible deviations from the quantum commutation relation 
even at the Planck scale. This would open the route to 
table-top quantum optics tests of quantum gravity pre- 
dictions. 

The possibility to interconnect optomechanical devices 
in large scale arrays has already been discuss ed in the 
context of investigating synchronization effects ( Heinrich 



et al. 2011|). The dynamics in such arrays may also 



enable the study of many-body quantum effects, which 
could complement the current efforts i n quantum simula- 
tions, yet in a solid-state architecture (Ludwig and Mar- 
[quardt', '2012 Tomadin et al. [2012") . Another scarcely ex- 
plored direction is to exploit the role of nonlinear mechan- 
ical responses. For example, the double-well potential 



of bistable mechanical resonators ( [Bagheri et al. 2011) 



could be the starting point for macroscopic tunnelling 
experiments. 

Finally, it is interesting to note that also in a broader 
context, the topic of controlled photon-phonon inter- 
action is receiving increasing attention. For example, 
a recent experiment has demonstrated quantum entan- 
glement between optical phonon modes of two separate 
macroscopic solids, specificall y of milimiter-sc ale bulk di- 
amond at room temperature (Lee et al. 2011 ), which has 



been generated from photon-phonon entanglement. Be- 
sides providing an interesting alternative to obtain quan- 
tum effects involving macroscopic objects this is also of 
direct relevance for applications, such as solid state quan- 
tum memories. Other examples include new schemes to 
achieve coherent conversion of bosonic modes, and even 
suggestions to exploit the measurement of optomechan- 
ical recoil e nergies for mesoscopic mass standards ( |Lan| 



et al. 2013). 



B. Applications 

Although the field of cavity optomechanics is only in 
its beginnings, several domains of applications have be- 
come obvious already now. In laser sciences, these in- 
clude tunable optical filters, based on the fact that op- 
tomechanical coupling can lead to extreme tuning of the 
mechanical frequency up to several octaves, as well a s op- 
tomechanical implementations of laser stabilization (Ale- 
gre et all 2010|. In addition, the compatibility of some 



optomechanical devices with silicon photonics (see Sec. 
IV ) enables on-chip optical architectures with added ver- 
satility. For example, exploiting strong optical nonlin- 
earities provided by optomechanical cavities (see Secjvj) 
adds an important and long-sought feature to optical in- 
formation processing. Along the same lines, embedded 
optomechanical cavities have been shown to serve as an 



all-optical memory element (Bagheri et al. 2011 Cole 



and Aspelmeyer 2011), or have been proposed as a new 



technology for single-photon detection (Ludwig et 



2012). In the first case, the binary states of a bistable 
nanomechanical resonator are controlled and monitored 
by an optomechanical cavity; in the latter case, a single 
photon would induce a measureable frequency shift to an 
optomechanical cavity, thereby resulting in a detection 
that is in principle destruction-free and photon-number 
resolving. In the long run, these features may provide 
a new momentum to all-optical information processing 
dCaulfield and Dolevl [20To| . 



For sensing applications, cavity optomechanics pro- 
vides several new aspects: for example, while damping of 
mechanical motion has been used to increase the band- 



width _o£_scajinm since decades (Bruland 

1996 Garbini et a/.[ [l996), cavity optomechanical 



et al 



devices allow both readout and damping of much higher 
mechanical frequencies, hence providing faster sampling 
and scanning rates. At the same time, the high sen- 
sitivity of the optical readout allows new integrated op- 
tomechanical platforms for acceleration sensing. A recent 
demonstration using optomechanical crystals achieved an 
on-chip acceleration resolut ion of 10 jiig/VH zwith. a test 



mass of only a few ng(Krause et al. 2012). Optome- 



chanical cooling in combination with on-chip mechanical 
sensors has recently also been suggested to provi de a re- 
duction in thermal noise for the optical readout (Winger 
\et al.\ |201l| . In turn, the ability to coherently amplify 
mechanical motion provides a route to radiation pres- 
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sure driven coherent oscillators with compact form factor 
and low power consumption. Finally, the combination 
with optomechanical preparation of squeezed mechani- 
cal states (see Sec X.A) could lead to a new mechanical 



sensing technology with unprecedented levels of sensitiv- 
ity due to the reduced position variance of the readout 
device. 

From a quantum information processing perspective 
cavity optomechanics offers a new architecture for co- 
herent light-matter interfaces in a solid-state implemen- 
tation. Mechanical motion can serve as a universal 
transducer to mediate long-range interactions between 
stationary quantum systems (see Sec X.E). The specific 



trait of optomechanical systems is the interconversion 
between stationary qubits and flying (photonic) qubits, 
which constitutes one of the main elements of long- 
distance quantum comm unication and a future quantum 
internet ( |Kimble| |2QQ8D . At the same time, strong op- 
tomechanical coupling in the single-photon regime opens 
up the field of non-Gaussian quantum optomechanics 
wit h a w ealth of quantum operations and protocols (see 
Sec X.F). The phenomenon of optomechanically induced 



transparency enables slowing of light pulses or even their 
storage, hence providing an interesting solid-state imple- 
mentation of a quantum memory (see Sec VII.B.2). 

Eventually, combining cavity optomechanics with 
other transduction mechanisms will allow to exploit the 
full functionality of micro- and nanomechanical devices. 
Such quantum hybrid systems utilize the mechanical mo- 
tion to achieve coupling between otherwise incompatible 
or uncoupled quantum systems (see Sec. X.D). A par- 
ticularly exciting perspective of opto-electromechanical 
hybrid devices is their ability for coherent conversion be- 
tween optical and microwave frequencies. Cavity cooling 
in these hybrid structures could also be applied to certain 
modes of a heat bath in integrated electronic circuits, for 
example to suppre ss unwanted therma lization effects in 
spintronic devices (Usami et al. 2012). Another inter- 



esting direction is to couple individual qubits, for exam- 
ple single atoms or single spins, to optomechanical de- 
vices. In combination with large mechanical frequencies 
such structures could allow mechanically mediated qubit 
interactions without additional laser cooling of the me- 
chanical modes, thereby significantly relaxing the exper- 
imental requirements for information processing in qubit 
registers. 
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XIII. APPENDIX: EXPERIMENTAL CHALLENGES 

A. Influence of classical excess laser phase noise on laser 
cooling 

We briefly consider the role of laser noise in optome- 
chanical cooling. Of particular interest is phase noise, 
described by the phase noise spectral density 5^00(^1) 
(or alternatively described by its equivalent frequency 
noise S^juji^) = S (f^cfyip) / Vfi) . Excess phase (and am- 
plitude) noise can be found in many laser systems due 
to relaxation oscillations, which can even in the case of 



diode lasers extend well into the GHz regime (Kippen- 



berg et a7^|2Qll[ |Safavi-Naeini et alj |2Q12| [Wieman and 



Hollberg 1991). Such excess phase noise has has been ex- 
perimentally observed to heat the mechanical oscillator 
(Schliesser et al] '20Q8b) and has been analyzed theo- 



retically (Diosil, 2008; Rabl et al 



Schliesser et al 



2008b ) . The spectral density of force fluctuations caused 



by this noise when pumping on the lower sideband in the 
resolved sideband regime, is given by (with r] = Kq^/k): 



Sujuj (^m) 



(123) 



By comparing this force noise to an effective ther- 
mal Langevin force of the laser {S^pp{^m) = 
2mef fTm^LJ^^m) near the mechanical resonance an 
equivalent laser noise occupation ul can be derived. The 
final occupancy of the mechanical oscillator in the pres- 
ence of optomechanical sideband cooling is subsequently 
nf = T^-^ — (fith + ^l), where nth denotes the average 

J- m~r^^ opt 

occupancy of the thermal bath. The excess contribution 
of the frequency noise is therefore: 



(124) 



The lowest occupancy that can be attai ned in the pres- 
ence of excess phase noise is given by ( |Rabl et al. 2009) 
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9 ^OJOJ 



(125) 



This expression can also be recast into a condition for the 
amount of phase noise that would lead to unit occupancy, 
i.e. Pq = 0.5. Given the cavity photon number required 
to cool near to the ground state (n/ = 1), we find that 
> 1 if the frequency noise exceeds the level ( |Rabl| 



et al 



2009) 



> 



si 



kBT/hQr, 



(126) 



If this level of noise is present, phase noise will preclude 
ground state cooling. This formula also reveals that to 
mitigate the effect of phase noise it is generally desirable 
to have a large vacuum optomechanical coupling rate and 
a low mechanical damping rate, to increase the level of 
tolerable phase noise. 



B. Influence of cavity frequency noise on laser cooling 

The cavity frequency fiuctuations are driven by ther- 
modynamical temperature fiuctuations. Considering the 
situation were light propagates inside a dielectric cavity 
(microtoroid, microspheres, photonic crystals) or pene- 
trates a mirror surface, the temperature fiuctuations of 
the sampled volume V are determined both by the abso- 
lute temperature and the heat capacity of that volume: 



ji2 



Here C is the specific heat capacity (per mass). In 
general, these temperature fiuctuations will exhibit a 
spectral density that depends on the resonator geome- 
try and boundary conditions for the thermal transport. 
For some cases, such as silica microsphere s, the spec- 
tral densities [cj] are known analytically ( Gorodetsky] 
and Grudinin 2004) and the corresponding frequency 
noise, S'^^ [cjj = K'^yi'^tt['^]i can be evaluated (where 



dT 



ujQ in the case of thermorefractive noise and 
Kth = If^^o in the case of thermoelastic noise). While it 
has been noted that on e can in principle com pensate also 
thermorefractive noise (Kimble et all 2008), the noise is 



of particular relevance to optomechanical cooling as it 
provides a limit to the minimal occupancy. This can be 
understood with the model referred to in the last section, 
ie. the fact that cavity frequency noise will translate into 
radiation pressure force fiuctuations in the presence of a 
strong cooling laser. 

Finite thermal cavity occupancy. A further source of 
cavity noise can arise from the fact that the cooling field 
has thermal occupation. This is in particular the case for 
microwave fields due to their low frequency. If the cavity 



occupation is given by n*^^, then the final occupation is 
modified to (in the resolved sideband limit): 



(127) 



This implies that the final phonon number can never be 
below the effective thermal occupation of the drive field 
(Dobrindt et al\ '2008"). It should be noted that when 
the radiation field and the mechanical oscillator initially 
have the same bath temperature (as will be the case in 
equilibrium, without extra absorption), the equilibration 
of these two oscillators of frequency cjcav and Vtrn will lead 
to an effective cooling of the lower frequency mechanical 
oscillator, as nth = ^^f^^ > ncav " ^^^^^ 



C. Influence of classical excess laser noise on sideband 
thermometry 

Excess noise of the readout laser (characterized by an 
occupancy ul ) compromises the self-calibration of the 
sideband thermometry method (Sec. VILA), as it can 
lead equally to an asymmetry, since in the presence of 
laser noise 7V^^^(A = Qm) = ^A+(n/ + 1 + nz,) and 
7V^^^(A = -1]^)^^^^^ = ^A-(nj-nL). This asymme- 
try can in this case be understood by arising from noise 
squashing and anti-squashing of the classical and quan- 
tum noises alike. This noise induced asymmetry on the 
upper and lower sideband can also be viewed as originat- 
ing from the effects of optomechanically induced trans- 
parency (OMIT) and optomechanically induced amplifi- 
cation (or EIA). 
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